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A USER'S  GUIDE 
TO  THE  OA3660  APL  WORKSPACE 


by 

F.  Russell  Richards 
Naval  Postgraduate  School 
Monterey,  CA  93940 

ABSTRACT 

Instructions  are  given  for  the  use  of  the  APL  Public 
Library  Workspace,  2 OA3660,  which  was  developed  as  an  aid  to 
interactive  exploratory  data  analysis.  The  OA3660  workspace 
is  accessible  to  all  users  of  the  computer  time  sharing  system 
at  the  Naval  Postgraduate  School.  The  workspace  contains 
various  data  analysis  functions,  data,  and  complete  internal 
documentation.  This  report  provides  a primer  on  APL,  documen- 
tation on  each  function  contained  in  the  workspace,  examples 
of  the  use  of  each  function,  and  program  listings. 
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A USER'S  GUIDE 


TO  THE  OA3660  APL  WORKSPACE 
by 

F.  Russell  Richards 
Naval  Postgraduate  School 
Monterey,  CA  93940 

I . INTRODUCTION 

This  report  provides  documentation  for  the  OA3660  APL 
workspace  that  was  developed  during  the  summer  and  fall  quarters 
of  1977  in  conjunction  with  the  offering  of  the  OA3660  course. 

Data  Analysis,  at  the  Naval  Postgraduate  School.  The  OA3660 
workspace  contains  APL  functions,  variables,  and  data  sets  that 
are  useful  for  interactive  exploratory  data  analysis.  The  APL 
functions  include  selected  procedures  from  the  STAT101  public 
library  [6],  and  from  Donald  R.  McNeil's  textbook.  Interactive 
Data  Analysis:  A Practical  Primer  [5]  (with  some  modifications). 

The  OA3660  workspace  also  includes  functions  written  by  myself 
and  students  in  my  class  and  various  utility  functions  extracted 
from  other  public  library  workspaces.  Documentation  for  each 
function  is  contained  in  this  report  and  in  the  OA3660  Public 
Library  Workspace. 
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The  data  sets  include  data  contained  in  examples  and  exercises 
in  the  above  mentioned  book,  data  from  John  Tukey's  Exploratory 
Data  Analysis  [8],  data  generated  from  test  scores  in  the 
OA3660  class,  and  data  extracted  from  a few  other  textbooks. 

The  data  are  cross  referenced  with  the  sources  in  a variable 
named  DATAMAP. 

This  report  is  intended  for  the  user  who  has  some  experi- 
ence with  APL  and  who  is  taking  a course  in  Data  Analysis  or  one 
who  is  already  familiar  with  the  basic  techniques  of  interactive 
exploratory  data  analysis  as  described  in  Tukey  [8],  McNeil  [5], 
and  Mosteller  and  Tukey  [7] . This  report  will  not  attempt  to 
explain  the  analysis  techniques  themselves,  nor  will  it  discuss 
detailed  APL  concepts.  It  will  give  a brief  description  of  the 
construction  of  data  arrays  and  transformations  of  data  since 
those  are  key  elements  in  interactive  exploratory  data  analysis. 

It  also  includes  a brief  discussion  of  basic  APL  commands  such 
as  logging  into  the  system,  loading  workspaces,  etc.  For  further 
details  on  APL  and  the  NPS  time-sharing  system,  CP/CMS,  the 
reader  is  referred  to  the  Naval  Postgraduate  School  Technical 
Note  No.  0141-33,  APL (CMS) -An  Introduction  [1];  the  IBM  report 
#GH20-0906-l,  APL\ 360-QS  and  APL\360-DQS  User’s  Manual  [2]; 

APL-An  Interactive  Programming  Language  [3]  by  Gilman  and  Rose; 
and  APL  Programming  and  Computer  Techniques  [4]  by  Katzan. 
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Each  function  contained  in  the  OA3660  workspace  is 
discussed;  the  function  syntax  is  given;  the  function 
parameters  (if  any)  are  described;  APL  listings  of  the  functions 
are  given;  and  at  least  one  example  is  included  to  illustrate 
the  use  of  each  function.  An  attempt  has  been  made  to  assure 
that  every  function  in  the  workspace  is  completely  debugged. 
However,  should  any  problems  be  experienced,  the  user  is  urged 
to  notify  me  of  the  problems.  Users  are  also  encouraged  to 
submit  to  me  interesting  data  and  APL  functions  that  are  use- 
ful for  interactive  exploratory  data  analysis. 


II.  APL  PRIMER 


This  chapter  is  written  primarily  for  the  user  who  is 
not  familiar  with  the  APL  programming  language.  We  attempt  to 
provide  only  that  material  that  the  user  needs  to  know  about 
APL  to  use  successfully  the  OA3660  APL  Public  Library.  There- 
fore, we  discuss  in  this  chapter  login  procedures,  basic  APL 
workspace  management,  error  recovery  procedures,  creation  and 
storage  of  data,  data  transformations,  function  syntax,  and 
logout  procedures.  A user  conversant  with  APL  should  proceed 
to  Chapter  III. 

A good  understanding  of  APL  would  enhance  the  user's 
facility  for  working  with  data  arrays  and  performing  interactive 
exploratory  data  analysis.  Therefore,  we  encourage  th  user  to 
seek  out  more  detailed  information  on  APL.  References  [1,2, 3, 4] 
are  all  recommended. 

A.  Getting  Started:  Login 

In  order  to  use  the  OA3660  APL  Public  Library,  the 
user  must  have  access  to  a computer  terminal  that  has  APL 
capability  linked  to  the  Naval  Postgraduate  School  IBM  360/67 
computer.  There  are  many  different  types  of  remote  terminals 
available  for  use,*  and  each  has  its  unique  features  of  oper- 
ation. Therefore,  the  user  should  check  out  the  operating 
* 

The  IBM  2741  terminals  with  a special  APL  typing  ball,  the 
Intertec  terminals,  and  some  of  the  CRT  terminals  have  APL 
capability.  In  addition,  there  are  several  school  owned 
portable  terminals  with  APL  capability. 
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instructions  for  each  terminal  to  augment  the  general  instructions 
given  here . 

If  the  terminal  is  wired  directly  to  the  computer,  the 

user  need  only  turn  it  on  to  access  the  NPS  CP/CMS  Time  Sharing 

System.  If  not,  the  user  must  link  to  the  computer  via  an 

acoustic  hookup.  First,  turn  the  terminal  on  to  the  correct 

settings.  Then  dial  the  appropriate  telephone  number  for 

* 

connection  to  CP/CMS.  When  the  shrill  audible  tone  is  received, 
place  the  telephone  receiver  into  the  acoustic  coupler  which 
is  connected  to  the  terminal.  If  all  of  the  terminal  settings 
are  correct,  and  if  the  time  sharing  system  is  in  operation, 
the  user  should  receive  the  message  "CP-67  ONLINE".  At  this 
point,  the  user  must  log  onto  the  system  by  typing  (with  the 
non-APL  character  set) : 

LOGIN  XXXXPYY  (Q) 

where  XXXX  is  the  user's  identification  number  assigned  by 
the  computer  center,  P indicates  a private  user  ID  (type  G 
if  a general  user  without  private  disk  space) , YY  is  the  terminal 
number,  and  (cr)  indicates  a carriage  return.^  If  the  identifi- 
cation number  is  valid,  the  system  will  request  the  user's  password. 
* 

The  telephone  number  and  the  terminal  settings  depend  on  the 
terminal  being  used.  For  ASCII  terminals  (nearly  all  terminals 
at  NPS  other  than  the  IBM  2741  terminals)  dial  either  x2611 
(with  speed  setting  at  110)  or  x3025  (with  speed  setting  at  300). 
For  EBCDIC  terminals  (IBM  2741),  dial  x2701. 

On  some  terminals  the  user  must  depress  simultaneously  the  keys 
CONTROL  and  S in  place  of  a carriage  return  even  though  there 
may  be  a key  so  labelled.  We  will  use  to  indicate  a 

carriage  return  for  all  terminals. 
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The  assigned  password  should  be  typed  followed  by  a carriage 
return : 

PASSWORD  (0) 

If  the  password  is  incorrect,  the  system  will  reject  the  login 
and  ask  the  user  to  start  over.  If  the  password  is  correct, 
the  system  will  request  the  user's  four-digit  project  number 
to  be  followed  by  a four-character  cost  center  code.  The 
project  number  must  be  assigned  by  the  computer  center.  The 
cost  center  code  is  the  user's  section  identifier  or  department 
code . 

PPPPCCCC  (0) 

If  the  project  number  is  acceptable  the  system  will  respond  with 
a ready  message  indicating  that  the  user  is  in  the  CMS  subsystem 
of  the  CP/CMS  Time  Sharing  System.  Otherwise,  the  user  must 
repeat  the  entire  procedure. 

Once  the  user  has  entered  CMS  he  should  switch  on  the 
APL  character  set*  and  type 

APL  (0) 

★ 

Switching  to  the  APL  character  set  may  consist  of  flipping  a 
switch,  changing  a type  ball  or  impact  print  set,  or  issuing 
a sequence  of  program  instructions  depending  on  the  type  of 
terminal . 
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to  enter  the  APL  subsystem.  The  system  will  then  respond  with 
the  message: 


A*P*L\N*P*S 

LIBRARY  DOCUMENTATION  SYSTEM...  ) LOAD  1 LIBDOC  TYPE  DESCRIBE.* 
The  user  is  now  in  APL,  and  all  of  its  powerful  features 
are  available  to  him.  He  may  use  APL  in  the  calculator  mode 
somewhat  as  he  would  use  a hand  held  calculator;  he  may  define 
his  own  functions  to  perform  a sequence  of  operations;  or  he 
may  access  the  public  libraries  which  consist  of  commonly  used 
preprogrammed  functions  that  perform  a variety  of  useful  compu- 
tations . 

In  the  next  section  we  describe  briefly  the  use  of  APL 
in  the  calculator  mode  to  generate,  manipulate  and  store  data. 

In  later  sections  we  describe  the  use  of  APL  to  transform  the 
data,  and  the  use  of  APL  Public  Library  Workspaces.  The  user 
is  directed  to  references  [1,2, 3,4]  for  information  about  writing 
functions  in  APL,  and  for  more  detailed  information  about  primitive 
operations  in  APL. 

★ 

If  any  symbol  other  than  the  right  parenthesis  ) appears  before 
the  word  LOAD,  the  system  is  not  properly  translating  the 
character  set  into  the  required  APL  characters.  If  this  happens, 
the  user  should  type  the  symbol  that  appears  (represented  here 
by  jt)  followed  by  the  word  OFF  as  follows: 

^off  (cr) 

This  will  get  the  user  off  the  system  at  which  time  he  should 
consult  with  the  computer  center  staff  as  to  the  special 
requirements  for  his  terminal.  Do  not  proceed  in  APL  until 
the  character  ")"  appears  before  the  word  LOAD  in  the  message 
typed  when  APL  is  entered. 
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B.  Generating  and  Storing  Data 

Generation  of  data  arrays  is  very  simple  in  APL  because 

one  need  not  be  concerned  with  format.  One  must,  however,  become 

familiar  with  the  APL  syntax.  First  is  the  assignment  symbol 

«-  which  plays  a role  somewhat  like  the  = sign  in  FORTRAN. 

It  means  to  take  the  expression  on  the  right  and  assign  it  to 

the  variable  named  on  the  left.  Suppose,  for  example,  one  wants 

* 

to  generate  a vector  of  data  named  X consisting  of  the  four 
observations:  3.15  12.57  8 6.003.  One  need  only  type: 

X 3.15  "12.57  8 6.003  (0) 

with  the  decimal  typed  if  needed  and  with  one  or  more  spaces 
serving  as  delimiters  for  the  separate  observed  values.  (Real 
numbers  may  be  expressed  in  scientific  format  by  use  of  the  E 
notation.  2.5E5  means  2.5  x 10^  or  250000.)  Extra  spaces  before 
and  after  the  assignment  arrow,  «- , are  not  needed,  but  they  will 
not  hurt  anything.  Additional  values  can  be  added  to  the  vector 
by  using  the  catenate  operator  (the  comma)  as  follows: 

* 

Just  about  any  name  can  be  used  for  variables  as  long  as  the 
first  character  is  alphabetic  and  the  other  characters  are 
alphabetic  or  numeric.  The  character  A is  also  acceptable. 

One  must  take  care  to  strike  the  characters  for  the  APL 
character  set.  Several  of  the  APL  characters  are  identical 
in  appearance  on  the  terminal  keyboards  to  non-APL  characters, 
but  different  results  are  obtained.  Among  the  characters 
that  appear  the  same  are  ) ' ( - + x 4 / : ; *. 
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X ^ X,  5.76  "4  1.47  (CR, 


This  statement  says  to  catenate  the  three  indicated  values  to 
the  end  of  the  old  data  vector  X and  to  call  the  resulting 
vector  X.  The  entire  vector  can  be  viewed  by  typing  X 
followed  by  the  ^r) 

X © 

3.15  "12.57  8 6.003  5.76  "4  1.47 

i 

i 

In  the  above  vector,  the  negative  signs  are  typed  using  the 
negative  symbol  (upper  shift  2) , not  the  minus  sign  (upper 
shift  +) . 

If  one  is  entering  many  data  values  into  a vector  he 
may  require  several  lines  of  input  to  do  so.  Additional  values 
could  be  inserted,  as  above,  using  the  catenate  operator.  How- 
ever, a simpler  way  is  to  type  a comma  followed  by  the  quad 
symbol  □ (upper  shift  L)  and  a carriage  return  at  the  end  of 
a line  if  more  values  are  to  be  continued  on  succeeding  lines. 

The  next  line  will  automatically  begin  with  the  quad  symbol 
prompting  the  user  to  continue  entering  values.  This  is 
illustrated  below: 

DATA  «-  4 6 8 12  6 "2 
□ : 4 2 12  9 "3  "17  9 9 2 "1  4 


'4  0 0 1 22  36  29  18 


, □ © 


If  a value  is  typed  incorrectly  and  discovered  before 

the  (cr)  is  struck,  the  value  can  be  corrected  by  backspacing 

* 

to  the  incorrect  entry  and  hitting  the  line  feed  key.  This 
will  erase  the  incorrect  value  and  everything  to  the  right  of 
the  value.  Then  type  the  correct  values. 

Data  correction  or  modification  after  a line  has  been 
terminated  by  a carriage  return  can  be  accomplished  several  ways 
One  way  is  to  determine  the  index  of  an  incorrect  entry  and 
assign  a new  value  to  that  specific  element  of  the  data  array. 
For  example,  suppose  one  wants  to  change  the  third  element  from 
8 to  4.81  leaving  all  the  other  elements  alone.  This  can  be 
done  by  typing: 

X[  3]  4.81  (0) 

Multiple  corrections  can  be  made  simultaneously  as  follows: 

X[  3,  4,  5]  4.81  5.9  “3.26  (cr) 

The  X vector  would  then  contain: 

3.15  “12.57  4.81  5.9  “3.26  1.47. 

The  index  operator  \ (upper  shift  I)  is  convenient  for 
determining  the  index  of  a given  value  (or  indices  of  a set  of 
values)  in  an  array  of  data.  If  X is  the  data  array  and  B 
is  a set  of  values  in  X,  then  X \ B will  generate  the  set  of 
a 

On  some  terminals  the  line  feed  key  is  labelled  LF. 

On  others,  the  attention  or  break  key  must  be  struck. 
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indices  of  the  values  of  B in  X.  If  a value  appears  in  more 
than  one  place  in  X,  only  the  subscript  of  its  first  appearance 
is  given.  For  example,  if  X +■  3.5  ”12.57  4.81  5.9  ”3.26  ”4 

the  index  of  the  value  5.9  is  found  as  follows: 


X \ 5.9 


4 


Let  B «-  ”12.57  ”3.26. 

given  by: 


The  indices  of  B in  X are 


X i B 

2 5 

This  operation  could  be  useful  for  altering  selected  values  in 
an  array.  To  illustrate,  assume  that  we  want  to  change  the 
values  ”12.57  and  ”3.26  in  X to  80  and  83.  The  operation 
below  will  accomplish  this  change: 

X[X  i ”12.57  ”3.26]  80  83 

X 

3.15  80  4.81  5.9  83  “4  1.47 

Additional  values  can  be  inserted  at  any  position  in  the 
vector,  and  values  can  be  deleted  from  the  vector  by  using  the 
take  "t",  drop  "4-",  and  catenate  operations.  The  take 


operation  with  syntax 


selects  the  first  r elements  from  the  vector  X if  r is 
positive  and  the  last  |r|  elements  if  r is  negative.  The 
drop  operation  with  syntax 


r 4 X 


deletes  the  first  r elements  of  X if  r is  positve  and  the 
last  I r | elements  if  r is  negative.  If  |r|  is  greater  than 
the  number  of  elements  in  X,  r t X will  insert  zeros  to  the 
right  of  the  elements  of  X until  |r|  elements  are  obtained 
and  r 4-  X will  result  in  the  empty  vector.  The  examples  below 
demonstrate  the  take  and  drop  operations. 


3.15  "12.57  4.81  5.9  “3.26  4 1.47 


3 t X 


3.15  “12.57  4.81 


'2  t X 


'4  1.47 


9 t X 

3.15  “12.57  4.81  5.9  “3.26  “4  1.47  3.15  “12.57 


2 4 X 

4.81  5.9  “3.26  “4  1.47 


'4  4 X 


3.15  “12.57  4.81 


8 4 X 


The  take  and  drop  operators  can  be  used  in  conjunction 
with  the  catenate  operator  to  edit  data  vectors.  The  examples 
below  illustrate  some  ways  this  can  be  done. 

Let 

X 1 2 3 4 5 6 7 8 9 10. 

Create  a new  vector  Y consisting  of  the  first  four  elements 
of  X followed  by  the  elements  11  12  and  then  the  last 


three  elements 

of  X; 

i 

.e . 

Y = 1 2 

3 4 

Y <- 

(4 

+ X)  , 

11 

12, 

("3  + X) 

© 

or  Y 

(4 

+ X)  , 

11 

12, 

(7  + X) 

© 

Let  W 0 ~1  ~2  -3  -4  -5.  Create  a new  vector  Y con- 

sisting of  the  first  eight  elements  of  X,  followed  by  the  number 
25,  followed  by  the  middle  four  elements  of  W. 

Y <-  (8  t X),  25,  (1  + (5  t W)  ) (0) 

A few  comments  about  the  above  operation  are  in  order.  APL 
always  operates  from  right  to  left  except  when  parentheses  are 
used  to  override  the  standard  order  of  operation.  Thus,  (5  + W) 
is  the  first  operation  executed.  This  results  in  the  vector 

0 ~1  “2  "3  ~4 . The  next  operation  is  (1+0  ~1  ~2  "3  ~4) 

which  gives  ~1  ~2  ~3  _4 . This  is  catenated  to  25  yielding 

25  -1  ~2  ~3  _4  which  is  then  catenated  to  (8  t X)  giving 

1 2 3 4 5 6 7 8 25  “1  “2  ~3  "4. 
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In  data  analysis  one  often  wants  to  compare  different 
groups  of  data  or  to  relate  one  group  to  another  group.  Many 
of  the  functions  in  OA  3660  require  the  data  arrays  be  matrices 
whose  columns  represent  the  different  groups  and  whose  rows 
represent  the  various  observations  for  each  group.  Let  us  see 
how  a data  matrix  can  be  created.  Suppose  for  example  that  one 
wants  to  create  the  matrix: 


The  easiest  way  to  create  A is  to  create  a vector  consisting 
of  the  eight  elements  and  then  reshape  the  vector  into  the 
desired  4x2  array.  The  reshape  operator  p can  be  used  to 
reshape  a vector  into  any  specified  size.  The  matrix  A is 

created  by  the  operation: 


A <-  (4 , 2)  p 12345678  (CR 


The  general  syntax  for  the  reshape  operation  is  A (N,M)pB 
where  B is  a vector,  or  scalar,  M is  the  number  of  columns  of 
the  array  A,  and  N is  the  number  of  rows  of  A.  The  first  M 
elements  of  B will  be  used  for  row  1,  the  second  M for  row  2, 
etc.  If  there  are  less  than  M‘N  elements  of  B,  the  elements 
will  be  repeated  from  the  beginning  as  many  as  are  needed.  If 
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B contains  more  than  M*N  elements,  only  the  first  M*N  will 
be  used.  Try  a few  examples  to  see  what  reshape  does. 

Remember  that  the  reshape  operator  fills  up  matrices 
row  by  row.  Suppose  that  the  user  entered  the  data  into  the 
vector  by  columns;  e.g.  B «-  1 3 5 7 2 4 6 8.  How  can  the 
4x2  matrix  A be  generated?  (4,2) pB  will  not  work  since 
that  would  give 


However,  if  one  requested  (2,4)pB  one  would  obtain 


13  5 7 
2 4 6 8 


What  we  need  now  is  to  transpose  this  to  obtain  A.  The  transpose 

★ 

operator  is  Si  (upper  shift  0 overstruck  with  upper  shift  / ) . 
This  sequence  of  operations  can  be  combined  as  follows: 


A «-  M2,4)pB  (CR 


The  operator  p is  also  useful  for  determining  the 
dimension  or  shape  of  an  arbitrary  variable.  When  used  for 
this  purpose  it  is  called  the  shape  operator.  The  syntax  is 


p Z (CR 


Type  shift  0 (oh) , backspace,  and  type  upper  shift  /. 
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Table  2.1  illustrates  the  results  of  the  shape  operator  when 
applied  to  selected  variables.  The  first  entry  is  the  scalar  3. 
APL  considers  all  scalars  to  have  no  dimension  (not  dimension  0). 
Thus,  pZ  is  the  empty  vector.  In  the  next  to  last  case,  the 
variable  Z is  defined  to  be  iO  (iota  0).  (This  represents 
a second  use  of  the  index  operator  a where  the  operator  has 
only  a right  argument.  If  K is  any  nonnegative  integer,  iK 
is  the  vector  consisting  of  the  first  K positive  integers. 

For  example,  i3  would  be  the  vector  1 2 3.)  The  use  of 

the  index  operator  in  Table  2.1  is  a rather  special  one.  APL 
interprets  t0  as  a vector  with  0 elements,  hence  the  result 
0 for  pZ.  This  probably  is  puzzling  to  the  reader,  but  it  is 
really  quite  handy.  The  reader  should  simply  commit  to  memory 
the  facts  that  scalars  have  no  dimension  and  the  empty  vector 
has  dimension  0.  In  the  last  entry  of  Table  2.1,  Z is  defined 
as  the  1x3  array  consisting  of  1 2 3.  This  array  appears 
to  be  identical  to  the  vector  Y «-  1 2 3,  but  APL  makes  a 
distinction  between  the  two. 


VARIABLE 

pZ 

Z 3 

- 

Z «-  2 3 5 

3 

Z «-  ( 2 , 2)  p 12  3 4 

2 2 

Z <-  iO 

0 

Z <-  (1, 3)  p 2 3 5 

1 3 

Table  2.1:  THE  SHAPE  OPERATOR 
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The  last  APL  operators  discussed  in  this  section  are 

the  grade  up  and  -ade  down  functions.  These  two  functions, 

* 

denoted  A and  f,  give  the  indices  of  the  elements  in  a 
vector  in  order  of  magnitude:  ascending  order  for  grade  up 
and  descending  order  for  grade  down.  This  is  very  useful  for 
data  analysis  since  much  of  our  work  requires  us  to  sort  the 
data.  These  operations  accomplish  this  easily.  Here  are  some 
examples : 

A <-  9 12  6 4 ~2  3 7 "5  11 

1A 

856437192 
f A 

291734658 

In  the  grade  up  of  A the  first  element,  8,  says  that  the 

eighth  value  of  A should  be  taken  first;  the  second  element, 

5,  says  that  the  fifth  element  of  A should  be  taken  next, 
etc.,  to  sort  A in  increasing  order.  Similarly  for  grade  down 
to  sort  in  decreasing  order.  The  sorted  arrays  can  be  written 
as  follows: 

AUA] 

”5  "2  3 4 6 7 9 11  12 

A[¥A] 

12  11  9 7 6 4 3 “2  “5 

•k 

These  symbols  are  upper  shift  H and  G,  respectively,  over- 
struck with  upper  shift  M. 
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We  have  only  scratched  the  surface  of  APL  features 


useful  for  generating,  modifying  and  storing  data.  The  collection 
of  operations  discussed  here  will,  nevertheless,  enable  a user 
to  do  most  of  the  things  that  he  requires  to  generate,  correct, 
modify  and  store  data.  In  the  next  section  we  discuss  arithmetic 
operations  in  APL  so  that  the  user  can  perform  calculations  and 
transform  data  arrays. 

C.  Arithmetic  Operations 

In  this  section  we  discuss  some  of  the  basic  arithmetic 
operations.  We  restrict  attention  to  those  operations  most  fre- 
quently used  in  data  analysis. 

The  symbols  used  for  the  basic  operations  of  adding,  sub- 
tracting, multiplying,  and  dividing  are  the  standard  ones  used 

* 

in  mathematics:  ( + , -,  x,  t).  However,  in  APL  these  operations 

are  more  powerful  than  their  equivalents  in  most  programming 
languages  since  the  operations  can  be  applied  to  entire  vectors 
or  matrices.  When  used  with  vectors  or  arrays,  the  operations 
are  applied  componentwise.  This  requires  that  the  vectors  or  arrays 
have  the  same  size.  An  exception  is  made  for  the  case  where  one 
of  the  arguments  is  a scalar.  When  this  happens  APL  adds, 

★ 

The  subtraction  sign  is  upper  shift  +.  Do  not  confuse  this 
with  the  negative  sign  which  is  upper  shift  2. 
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subtracts,  multiplies  or  divides  the  scalar  to  each  element  of 
the  vector  or  array.  Of  course,  the  operation  can  be  applied 
to  two  or  more  scalars  as  with  any  other  programming  language. 
The  examples  below  illustrate  the  basic  operations: 


A «-  3 5 7 8 

B •*-  ~1  "2  1 2 


A + B (CR 


2 3 8 10 


A - B (CR 


4 7 6 6 


A x B (CR 


'3  "10  7 16 


A r B (CR 


'3  "2.5  7 4 


2 x A (C 


6 10  14  16 


3 + B (CR. 


2 14  5 


C <-  (2,2)  PA  (C 


C (C 


7 8 


D ( 2 , 2 )p  B (C 


D 


1 2 
1 2 

2 3 

8 10 

1.5  2.5 

3.5  4.0 

27 


C + B 


C v 2 


3x4  + 5 


12 


3x2 


The  last  two  examples  are  included  to  re-emphasize  the 
right- to-left  order  of  operations  in  APL.  Thus,  5 is  added  to 
4 to  give  9 which  is  then  multiplied  by  3 to  yield  27.  Similarly, 
in  the  last  case,  2 is  multiplied  by  3 to  give  6 which  is  then 
divided  into  12  to  yield  2.  Unlike  some  other  languages  there 
is  no  hierarchy  of  operations  other  than  right  to  left.  Of 
course,  parentheses  may  be  used  as  in  algebra  to  change  the 
priority  of  operation. 
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(3x4)  +5 


(12  4 3)  x 2 (CR 


8 
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In  addition  to  the  four  basic  operations,  data  analysis 
frequently  requires  power,  log,  reciprocal,  and  exponential 
transformations.  As  with  the  four  basic  operations  these  can 
be  applied  to  scalars,  vectors,  or  matrices. 

Table  2.2  presents  the  syntax  for  these  operations. 


Trans  formation 

Syntax 

xa 

X * a 

in  X 

®X  1 

1/X 

tX 

X 

e 

*X 

Table  2.2.  SYNTAX  FOR  DATA  TRANSFORMATIONS 


Here  we  give  examples  of  these  transformations: 

X 1 2 3 4 5 © 

X * 0.5  (cr) 

1 1.41421  1.73205  2 2.23607 

X * 2 (CR) 

1 4 9 16  25 

X * ~3  © 

1 0.125  0.037037  0.015625  0.008 


1 The  symbol  for  the  logarithm  is  the  upper  shift  0 (oh)  over- 
struck with  the  asterisk  (upper  shift  P) . 
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0 0.693147  1.098612  1.386294  1.609438 


1 0.5  0.333333  0.25  0.2 

X * * “3  @ 

1 0.793701  0.693361  0.629961  0.584804 

(Recall  that  APL  operates  from  right  to  left  so  that  the  operator 
above  raises  each  element  in  X to  the  -1/3  power.) 

*X  ^r) 

2.71828  7.389056  20.085537  54.598150  148.413159 

The  logarithm  of  a number  N to  an  arbitrary  base  B can  be 

determined  by  typing  B 0 N;  e.g., 

10  0 1000 


3 

Three  other  useful  arithmetic  operations  for  scalars, 
vectors,  and  matrices  are  the  ceiling  [",  floor  L,  and 
absolute  value  | . All  of  these  can  be  used  with  a single  right 
hand  argument  or  with  both  left  and  right  arguments. 

The  examples  below  illustrate  the  use  of  these  operators: 


Monadic  (one  argument) 
r 2 . 75  6 0.08  "3.6 

3 g i "3  (Gives  smallest  integer  £ argument. 

L 2.75  6 0.08  "3.6 

2604  (Gives  largest  integer  £ argument.) 


2.75  6 0.08  "3.6 

(Gives  absolute  value  of  each  element. 


2.75  6 0.0  8 3.6 


3 5 r 2 8 


3 8 (Gives  the  maximum  of  each  component.) 

3 5 I.  2 8 

2 5 (Gives  the  minimum  of  each  component.) 

3 I 7 6 8 "2 

★ 

1021  (Gives  the  remainder  after  dividing  each  element  by  3.) 

Clearly,  the  floor  and  ceiling  operators  are  useful  for 
rounding  values.  One  can  select  any  number  of  significant  digits. 

If  one  wants  to  round  values  to  the  nearest  integer,  one  should 
type  L0.5  + N.  If,  say,  four  significant  digits  are  wanted, 
one  could  type  (L.5  + N x 1E4)  i 1E4.  The  examples  below 
illustrate  the  use  of  |_  for  rounding: 

N 0.0835126  12.51877623  1.33333333  5.25 

L .5  + N 
0 13  1 5 

(L.5  + N x 1E4)  i 1E4 
0.0835  12.5188  1.3333  5.2500 

Another  operator,  the  reduction  operator  /,  allows 
the  operations  discussed  above  to  be  applied  to  all  the  elements 
of  a vector  or  to  the  rows  or  columns  of  a matrix.  The  syntax  is 

* i 

All  remainders  of  K|N  are  expressed  as  positive  integers  in 
the  set  0,  iK-1.  If  a remainder  is  negative,  K is  added  to 
it . 
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f/A  where  f can  be  any  of  the  arithmetic  operations  discussed 
above  and  A is  a vector  or  a matrix.  The  result  of  this  oper- 
ation when  applied  to  a vector  A of  size  n is  the  scalar 
A [1 ] fA[2]f  •••  fA[n].  When  applied  to  a matrix  M having  r 
rows  and  c columns,  the  result  is  the  vector  B of  dimension  r 
where  B[i]  = A[i;l]  fA[i;2]  f •••  fA[i;c]  for  i = 1,2, ...,r 
(column  reduction) . The  matrix  M can  be  reduced  over  its 
rows  by  typing  f/[l]M.  The  examples  below  illustrate  the  reduction 
operation : 

+ /I  3 7 4 

15  (the  sum  of  all  the  elements) 

x/1  3 7 4 

84  (the  product  of  all  the  elements) 

r/1  3 7 4 

7 (the  largest  element) 

M «.  ( 3,  2)  p 3 1 2 6 5 4 

M 

3 1 

2 6 
5 4 

-/M 

2 41  (differences  of  elements  in  rows  1,  2,  and  3) 

L/M 

124  (minimum  values  in  rows  1,  2,  and  3) 


L 
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+/  [ 1 ] M 


10  1 


5 6 


(sums  of  elements  in  columns  1 and  2) 

r/[l]M 

(maximum  values  in  columns  1 and  2) 


Finally,  we  describe  the  inner  product  operation  and 
the  matrix  inverse.  The  inner  product,  like  the  reduction 
operation  can  be  applied  with  any  general  APL  operators.  The 
syntax  is  Af.gB,  where  A and  B are  vectors  or  matrices 
(which  must  satisfy  certain  size  restrictions)  and  f and  g 
are  any  general  APL  operators.  For  vector  arguments,  A and  B 
must  be  the  same  size,  say  n,  and  the  result  is  given  by 
f/A[l]gB[l]  A[2]gB[2]  ...  A[n]gB[n],  For  matrix  arguments 

the  number  of  columns  of  A must  be  the  same  as  the  number  of 
rows  of  B.  The  result  is  a matrix  of  size  n x m where 
pA  = (n,k)  and  pB  = (k,m) . The  (i,j)th  element  of  the  result 
is  f/A [i ;1 ] gB  [1; j ] A [i; 2 ] gB  [2 ; j ] ...  A [i ;k ] gB  [k ; j ] . The 

reader  should  recognize  the  operation  above  as  matrix  multi- 
plication when  f is  + and  g is  x.  The  examples  below 
illustrate  the  inner  product. 

X <-  1 4 3 2 

Y «-  2 3 1 2 

X + . x y 

21  ((1x2)  + (4x3)  + (3xl)  + (2x2)) 

X x . + y 

336  ((1  + 2)  x (4  + 3)  x (3  + 1)  x (2  + 2)) 


A 3 2 p i 6 

A 

1 2 

3 4 

5 6 

B«-2  3p3  5 2 3 1 1 

B 

3 5 2 

3 11 

A + . x B 


7 

7 

4 

(the  matrix  product  of  A and  B) 

21 

19 

10 

33 

31 
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The  domino  operator  i (type  0,  backspace,  *)  is  used  to 
solve  for  the  matrix  inverse  of  a nonsingular  square  matrix 
The  inverse  of  a nonsingular  square  matrix  C is  found  by 
typing  SC.  For  example, 

C-<-22pll0  2 

C 

1 1 
0 2 

sc 

1 -0.5 

0 0.5 

C + . x sc 

1 0 
0 1 

28 
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D.  Workspace  Management 


In  the  previous  sections  we  have  seen  how  data  are  input 
and  variables  are  transformed.  In  this  section  we  describe  how 
the  user  manages  his  APL  workspace  so  that  he  can  save  data  for 
use  from  one  session  to  the  next;  he  can  load  public  library 
workspaces;  and  he  can  use  functions  available  in  the  public 
library  workspaces.  For  information  about  writing  functions  the 
user  should  see  references  [2,3,4]. 

When  a user  types  APL  he  is  put  into  a clear  APL  work- 
space. In  this  clear  workspace  data  can  be  created  and  any  of 
the  operations  described  in  the  earlier  chapters  can  be  performed. 
If  the  user  wishes  to  maintain  the  data  or  results  for  future 
use  he  must  save  the  workspace.  This  is  done  by  typing 

) SAVE  WSNAME 

where  WSNAME  is  an  arbitrary  name  (first  letter  alphabetic, 
eleven  characters  or  less)  that  the  user  selects  for  the  workspace. 
This  private  workspace  can  subsequently  be  loaded  into  the  user's 
active  workspace  by  typing 

) LOAD  WSNAME 


This  causes  the  active  workspace  to  be  cleared  and  a copy  of 
the  named  workspace  to  be  written  into  the  active  workspace. 
Alternatively,  if  the  user  does  not  want  to  clear  out  the  contents 
of  the  active  workspace  to  bring  in  the  named  workspace,  he  can 
type 
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) COPY  WSNAME 


r 


This  will  simply  augment  the  existing  contents  of  the  active 
workspace  with  the  contents  of  workspace  WSNAME.  However,  if 
an  existing  variable  or  function  in  the  active  workspace  has  the 
same  name  as  a function  or  variable  in  workspace  WSNAME  the 
latter  will  replace  the  former. 

Any  modifications,  additions,  deletions,  corrections, 
etc.  that  the  user  makes  to  a copy  of  a workspace  will  affect 
the  active  workspace,  but  will  not  affect  the  'permanent  copy  of 
the  workspace  maintained  on  his  private  files  unless  the  user 
saves  the  so  modified  active  workspace.  This  can  be  accomplished 
by  typing  ) SAVE  WSNAME  (or  simply  ) SAVE  if  workspace  WSNAME 
were  LOADed) . 

The  user  can  determine  the  name  of  the  active  workspace 
at  any  time  by  typing  )WSID  (for  Workspace  IDentif ication) . 
Similarly,  he  can  change  the  name  of  the  active  workspace  by 
typing  )WSID  NEWNAME . A list  of  all  the  user's  private  APL 
workspaces  is  obtained  by  typing  )LIB.  Entire  workspaces  can 
be  permanently  destroyed  by  typing  ) DROP  WSNAME.  (Be  careful 
with  this  one!) 

Within  an  active  workspace,  a user  can  obtain  a list  of 
all  of  the  functions  contained  in  the  workspace  by  typing  ) FNS . 
The  functions  will  be  listed  alphabetically  so  that  all  functions 
from  those  beginning  with  a certain  letter,  LETTER,  onward  can 
be  obtained  by  typing  ) FNS  LETTER.  Similarly,  a list  of 
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variables  in  the  workspace  can  be  obtained  by  typing  )VARS  or 
) VARS  LETTER.  Typeout  can  be  terminated  at  any  point  by 
hitting  the  BREAK  or  the  ATTN  key.  Variables  and/or  functions 
can  be  deleted  from  the  active  workspace  by  typing  ) ERASE  LIST 
where  LIST  is  a single  function  or  variable  or  a list  of  functions 
and  variables  to  be  deleted.  Names  in  the  list  should  be 
separated  by  one  or  more  blanks.  Two  other  system  commands 
allow  the  user  some  control  over  his  APL  environment.  These 
commands  control  the  number  of  digits  printed  out  and  the  width 
of  a typed  line.  They  are  ) DIGITS  N and  ) WIDTH  N,  where  N 
is  the  number  of  digits  to  be  printed  or  the  desired  line  width, 
respectively . 

In  addition  to  his  own  private  APL  workspaces,  the  user 
also  has  access  to  all  of  the  APL  public  libraries  available 
at  NPS . These  libraries  have  numbers  between  1 and  999  and  are 
intended  to  hold  workspaces  of  general  interest.  See  [1]  for 
a list  of  the  public  library  numbers.  The  contents  of  a public 
library  can  be  displayed  using  the  )LIB  command  followed  by 
the  library  number.  For  example. 
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) LIB  1 


A-DISK 

R/0 

IOFNS 

17.25 

7/09 

PLOTFORM 

17.25 

7/09 

NEWS 

15.54 

8/2  7 

WSFNS 

17.25 

7/09 

TEXTEDIT 

17.26 

7/09 

FORMAT 

17.26 

7/09 

CATALOG 

15.17 

7/14 

MAILBOX 

12.50 

7/14 

MULTIPLO 

14.06 

7/18 

FILEFNS 

17.18 

8/14 

The  list  contains  the  workspace  name  and  the  time  and  date  that 
it  was  last  modified. 

A copy  of  a public  library  workspace  can  be  put  into  the 
user's  active  workspace  by  typing  ) LOAD  n WSNAME  or 
) COPY  n WSNAME  where  n is  the  library  number  and  WSNAME  the 
workspace  name.  (Recall  the  differences  in  LOAD  and  COPY;  LOAD 
will  first  clear  the  contents  of  the  active  workspace.)  Selected 
functions  and  variables  from  a workspace  (either  public  or 
private)  can  be  copied  by  typing 

) COPY  n WSNAME  OBJECT 
or 

) COPY  n WSNAME  GROUP 
32 


where  OBJECT  is  a single  variable  or  function  that  the  user  wants 
to  copy  into  his  active  workspace  and  GROUP  is  a group  of  variables 
and/or  functions  that  has  previously  been  defined  for  WSNAME.* 

When  a public  library  workspace  is  brought  into  the 
user's  active  workspace,  documentation  can  usually  be  obtained 
by  typing  the  word  DESCRIBE.  It  is  somewhat  standard  procedure 
to  document  public  library  workspaces  with  a DESCRIBE  variable 
giving  general  information  and  with  each  function  in  the  work- 
space being  documented  by  a "HOW"  variable.  Type  the  name  of 
the  function  followed  by  HOW  with  no  intervening  spaces.  For 
example,  documentation  for  the  function  STEMLEAF  in  public 
library  2 OA3660  is  obtained  by  typing  STEMLEAFHOW . "HOW" 
variables  generally  describe  the  function  syntax,  parameters, 
input/output  requirements,  etc.  Of  course,  the  user  can  determine 
if  a workspace  contains  this  sort  of  internal  documentation 
by  typing  )VARS. 

When  the  user  has  completed  his  work  session  he  can  logout 
by  typing  )OFF.  If  he  wants  to  get  out  of  APL  but  not  logout, 
he  should  type  ) OFF  CMS.  Then,  he  should  switch  back  to  the 
standard  keyboard. 


The  group  structure  for  the  OA3660  workspace  is  discussed  in 
Chapter  III. 
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E.  Error  Recover 


This  section  describes  recovery  procedures  that  the  user 
can  employ  if  he  receives  error  messages  during  execution 
of  functions  contained  in  a workspace.  Usually,  error  messages 
received  when  functions  in  public  library  workspaces  are  executed 
are  the  result  of  improper  syntax  or  problems  with  the  shapes 
of  the  arguments.  As  soon  as  an  error  is  detected,  execution 
of  the  function  is  suspended,  the  number  of  the  line  containing 
the  error  is  typed,  a caret  is  inserted  at  the  position  in  the 
line  that  the  error  occurred,  and  an  explanatory  error  message 
is  printed  out.  When  the  function  is  suspended  the  values  of 
all  variables  determined  up  to  the  point  of  suspension  can  be 
obtained  by  simply  typing  the  variable  name.  All  of  this  makes 
error  discovery  and  correction  quite  simple.  From  the  point  of 
suspension  the  user  can  branch  to  any  line  of  the  function  by 
typing  -*■  n,  where  n is  the  desired  line  number.  If  n is 
omitted  a branch  is  made  outside  of  the  function.  Many  types  of 
calculations  can  be  performed  while  a function  is  suspended, 
including  the  execution  of  other  functions.  However,  the  function 
will  remain  suspended  until  a branch  is  made  as  described  above. 
Suspended  functions  tend  to  clutter  up  the  user's  workspace  so 
that  suspended  functions  should  not  be  left  pending.  A list 
of  all  suspended  functions  can  be  displayed  at  any  time  by 
typing  )SI  (state  indicator).  If  several  functions  are  listed 
as  being  suspended,  the  user  should  type  as  many  branch  arrows 
(-►)  , one  per  line,  as  there  are  asterisks  displayed  in  the  list. 
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Since  the  functions  contained  in  the  public  library 
workspaces  have  been  tested  extensively,  most  of  the  errors 
that  are  encountered  result  from  improper  function  syntax  or  the 
use  of  arguments  that  are  not  conformable  or  of  improper  size. 
The  documentation  contained  in  the  next  chapter  shows  the  proper 
syntax  for  each  function,  and  gives  the  requirements  placed  on 
the  arguments. 
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III.  THE  OA3660  WORKSPACE 


This  chapter  provides  documentation  for  the  public 
library  workspace,  2 OA3660.  General  workspace  documentation 
is  contained  in  the  DESCRIBE  variable  and  in  the  lists  of 
functions,  variables,  and  groups.  Short  writeups  are  given 
for  each  function  in  the  workspace  via  "HOW"  variables,  and 
examples  of  the  use  of  each  function  are  provided.  Data  are 
contained  in  the  workspace  to  provide  easy  illustration  of  the 
functions.  The  data  are  described  in  the  DATAMAP  variable. 

Finally,  the  actual  APL  programs  are  displayed. 

A.  General  Documentation 

The  OA3660  APL  workspace  is  contained  in  public  library  2. 
Therefore,  the  user  must  type  ) LOAD  2 OA3660  or  ) COPY  2 OA3660 
to  create  a copy  of  the  OA3660  workspace  in  his  active  workspace 
area.  The  functions  and  variables  contained  in  OA3660  are 
displayed  by  typing  ) FNS  and  ) VARS  as  shown  below. 

) PNS 

A,JL/  A NOV  A A3h  A3hSk  BOX  PLOT  CiilSQUAkk  COUok-,S 

COMPAhk.  CONDENSE  COnTI t'JObUCY  COhttaLATTO N FILL  <■  \T 

I vfJi  kb  LI  kb  LIT  LS  LI  HE  fab  GPO  LTSti  mSTATS  iUXS’J  ■ 

Oi-Jbd  ON  63  PARTIAL  h LG  HESS  SCAT  SHOW  h US  SPLIT  STATISTICS 

STbXLSAF  SUMS*  TWICE  UTCOND 


)VARS 

AGES  ANDHOW 
CtiICKWTS 
CON  DENSE  HOW 
C RIMES  DATAmAP 
DRAPER  EPSILON 
LAUD  ARRAS 
MORTALITY 


ANOVAHOW 
CHISJUAREHOW 
CONSUMPTION 
DEATHS  uEP 
GRADES  1 GRADES2 
LENGTH  LINE  HOW 
MST  ATSHON 


A 3 R HOI  AZRSRHOw  BOXPLOTaOW 

COAL  CODE  RES SOU  COMP  ARP.  HU -t 

CONTINGENCY  HON  CO  RUE  LATI  Oft  H 04  Cf: 
DEPTH  DESCRIES  DISCOVERIES 


HY DROP  LASTS 
LSLINEHOW 
NDIVX  NOISY 


insects  kshow 
MEDPOLI  SilHOW  ■ iISS 
V G A P N 0 ft  -l  ate  C T S 


NU.-l  NUMSUMHOW  OH  Ed  HON  ONES  HON  OPTION  PtiRTT  A LHOW 

PRECIPITATION  PRESRATTNG  REG RESSHOW  RESIDS  RIVERS 

RSeLECT  SCALE  SC  AT  HON  SHOWRESHOW  SPLITHOW  STATISTIC 

STun LEAF  HOW  TESTX  TESTY  TWICE HOW  USPOP  VOLCANO 

WARPBREAKS  WID  WIDTH  ^INTERCEPT 


i 0 v 


Notice  that  there  is  a variable  called  DESCRIBE.  This  variable 
gives  general  documentation  about  the  workspace. 


fl T H E DESCRIBE  VARIABLE  PROVIDES  GENERAL  OOCJhE  STATION 
fit  OH  THE  ENTIRE  WORKSPACE.  SIMPLY  TYPE  ’ DESCRIBE ' TO 
n OBTAIN  THE  I N FORMAT I OH . 


U ESC RI EE 

OAZ660-DATA  ANALYSIS 

DATE:  JANUARY  197  8 

PROGRAMMER : F.  RUSSELL  RICH  ARDS 


THIS  WORKSPACE  COST  AIRS  FUNCTIONS  AivD  VARIABLES  USEFUL  V 
EXPLORATORY  INTERACTIVE  DATA  ANALYSIS.  FOR  INFORMATION  On  T.ie 
USE  OF  TnE  FUNCTIONS  TYPE  THE  FUNCTION  NAME  FOLLOWED  bY  ' ROW ' . 
FOR  EXAMPLE,  TYPE  ST  EM  LEAF  HOW  FOR  DOCUMENTATION  ON  THE  FUNCTION 
STerLEAF.  SOME  FUNCTIONS  ARE  USED  ONLY  AS  SUBPROGRAMS  FOR  OTHER 
FUNCTIONS  And  WILL  NOT  HAVE  a 'HOW'  VARIABLE.  VCR  :0\ iE  Detailed 
DOCUMENTATION  El, a THE  NFS  TECHNICAL  REPORT,  A J SEN'S  GUIDE  TO  7 
•j  A 3 6 b 0 A PL  WORKSPACE , BY  t . RUSSELL  RICHARDS  OR  DOW  A LD 
XCneIL'L  u 0 0 u , I H T E R AC l J V E DATA  A N A Lie  I D • 

THE  WORKSPACE  ALSO  CONTAINS  VARIOUS  SMALL  DATA  SETS  VUE 
ILLUSTRATION  OF  The  FUNCTIONS.  type  DATAmAP  TO  OBTAIN 
IDENTIFICATION  OF  THE  DATA  AhhAYS. 

THE  JSER  CAN  SAVE  SPACE  IN  HIS  ACTIVE  WORKSPACE  t> i Cot' t T « / 
OnLY  THOSE  functions  AuD  V A RT A b Le S jEEDEd  TO  p»  DFQhs  i i 

. uC  eSS  ARY  TASKS.  FUNCTIONS  A J D VARIABLES  ARE  CON  V v.  iTEHl  IM 
COLLECTED  into  groups  which  are  identified  IN  THv,  'Hot' 

VARIABLES.  V OR  example.  It  THE  USER  WANTS  TO  i E V E RATE  ST-  .LA-’ 
Jn  tsOKPLOT  DISPLAY  S He  CAN  COPY  THE  DISPLAY  GROUP  AS  FOLLOW  S : 

)CDPY  2 CM  3 660  DISPLAY 

THE  GROUP  CONTAINS  JOT  ONLY  THE  FUNCTIONS  BUT  ALSU  ALL 
REJOINED  PARAMETERS  AND  DOCUMENTATION . 
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There  is  also  a variable  named  DATAMAP  which  describes  the 
data  contained  in  OA3660  and  provides  a reference  to  the  data  source. 
This  is  shown  below. 


n'/'ttfc  / ARIA&L&  DAT  Ae-iAP  PROVIDES  T UFOR'.'IATT  U‘J  AttG’Ji  7 /•> 
no  AT  A S STS  COL  T AT  USD  Tit  THE  UOHKSPACb  . TT  GT.pS  C.tn 
ftSOURCE  UP  EACH  DATA  SET  AHO  Tun  Sj  Zs  OP  j DA‘i  4 /■ i . 
nST.\PLY  TYPe  ’ DAT  AiiAP  ’ DAT  A DOCJkS  PC  ATT  Oh  . 


J/qi/hi/lP 


TdTS  YOReSPACE  CONTAINS  Se 

V ORAL  DATA 

S E TS  THAT  C A a OE  ‘J  S e O t ' > 

ILLUSTRATE  THE  FUNCTIONS  Oh  TO  SERVO 

AS  UQ'-tR  t/OkK  PHODLE  tS.  7:.. 

E'hTIrl  C 0 ELECT  TO  it  OE  DATA 

ARRAYS  TS  C 

OUT  A T J ED  I . t 7 ::  E GnOJ  P •'*/:  . 

DATA.  THE  DATA  C AhE  t'hOi 

THE  FOLLOTTlvG  SOURCES: 

( 1 ) Oh  A P E ll , it . h . A 0 D H 

. S.ilTH,  APPLIED  a STRESS!  Oh  /.,*  ALYS. 

■J  T LET  A a D SOUS. 

12)  'iC 0 bT L . U . h . , TO? E 

RACTTVE  u AT 

A ANALYSTS,  A DDT  SO  « *»V  ESLeY  • 

(3)  -ll  LLEn.  T.  ADD  J.L 

. PREUuu, 

PhOoAsT  LT'TY  A.n)  STATISTICS 

POh  ENGINEERS,  PhEnTT 

C E-' n ALL . 

( 4 ) h T C ci  A h DS  , F . k . , C LASSROOk  GRADES. 

( 5 ) TJ  KEY  , -JOd.t,  EXPLORATORY  DATA 

AH ALY ST S , AUDI SG:'t  - v’t  SLR  Y . 

Tttb  TAbLE  bbLO.Y  LTSTS  EACH 

u AT  A ARRAY 

, TttR  UT  .e  vSTOd,  T.i:.  SOU  - C 

(USING  AUTltOh  AtitirE/ T ATT  Olt 

5).  ADD  PAG 

E Or.  CHAPTER  :<U  ib  mS. 

DATA 

no  AT  A 

SGj  .‘i  C o 

A 0 E S 

42 

■<C(R.11) 

C it  TCKrtTS 

14  6 

<C(P. 30) 

CUAL 

49 

TU 

COitSU'-iPT  T ON 

5 5 

hC(P.  100) 

C kT 4 ES 

50  4 

MC (P.132) 

DEATHS 

5 4 

f!C(P.  94  ) 

DT SCOV  r,  RT ES 

100 

■iC(P.  121) 

On  APE  EX 

1 3 4 

DS  ( P . 1 7 8 ) 

DRAPE hY 

3 3 

DS(  o.  l ->8  ) 

GhADESl 

15  5 

RT 

G RAO  SS2 

15  5 

RT 

HY  DROP  LAnT  S 

34 

Tu(Ctt.  3 ) 

T D SECTS 

12  6 

PC(P. 12,36) 

L Alt  DA  RE  AS 

4 8 

,C(P.  9) 

40RTALTTY 

5 4 

T'J(CH.  15) 

PRECT  PTT  AT  T Oil 

69 

HC(P . 3 ) 

phesratthg 

1 14 

HC(P . 126) 

RT  V E R S 

1 4 1 

kC ( P . 1 4 ) 

TbSTK 

10  2 

Ht'(P.  253  ) 

T r,  STY 

10 

Kt  (P.  253) 

'JSPOP 

1 'J  2 

TJ(Crt.8) 

V 0 DC  ADO 

219 

TU 

-v  Ah  P b R R A ri S 

9 6 

SC (P .28) 

Finally,  notice  that  there  are  "HOW"  variables  for  most 
of  the  functions  contained  in  the  workspace.  (The  only  functions  with- 
out "HOW"  variables  are  utility  functions  which  are  used  by  other 
functions,  but  which  are  transparent  to  the  user.)  These  variables 
provide  documentation  on  the  use  of  the  functions.  For  example, 
documentation  on  the  function  SCAT  is  obtained  by  typing 
SCATHOW. 


SCATHOW 

SURTAX'.  SCAT  ARRAY 
PARAMETERS  : 


(1)  WID-  CONTROLS  THE  HORIZONTAL  SIZE 
( DEFAULT =30  CHARACTERS) . 

OF 

THE 

DISPLAY 

(2)  DEP - CONTROLS  THE  VERTICAL  SIZE 
( DEFAULT- IB  LINES). 

OF 

The 

DISPLAY 

(3)  NDIVX,  NDIVY - NUMBER  OF  UNITS  ON 
RESPECTIVELY  ( DEFAULTER, 4) . 

X- 

AND 

Y-AXES, 

GROUP:  RELATIONS,  SMOOTH , COMPARISONS 

DESCRIPTION:  SCAT  PRODUCES  A SCATTER  PLOT  OP  Tub  DATA  CONTAINED 

IN  ARRAY.  THE  ARGUMENT  ARRAY  CAN  BE  A VECTOR  OR  A MATRIX  PITH  AS 
MANY  AS  9 COLUMNS.  IP  A VECTOR  OF  SIZE  N,  THOSE  VALUES  ARE 
PLOTTED  VS.  THE  INTEGERS  1 TO  N;  IF  A MATRIX,  THE  SECOND,  THIRD, 
ETC.  COLUMNS  ARE  PLOTTED  VS.  COLUMN  1 ON  THE  SAME  AXES. 
DIVISIONS  ON  THE  AXES  OF  THE  PLOT  ARE  NOT  EXPLICITLY  PRINT  ED , 
EXCEPT  AT  THE  EXTREMES  OF  THE  PLOT.  THE  USER  CAN  CONTROL  THE 
RESOLUTION  OF  THE  PLOT  BY  MODIFYING  THE  PARAMETERS  WID,  UoP, 
NDIVX  AND  NDIVY . PRINTING  TIME  INCREASES  DRAMATICALLY  V IT n 
RESOLUTION ; THEREFORE,  OTHER  PLOT  PROGRAMS  SHOULD  BE  USED  IF  HIGH 
RESOLUTION  IS  DESIRED.  FOR  A SINGLE  GROUP  OF  DATA,  THE  NUMBER  r. 
(2<fliS9)  HILL  BE  PRINTED  IF  N POINTS  LIE  CLOSE  TOGETHER  ON  THE 
DISPLAY.  FOR  MULTIPLE  PLOTS  ON  THE  SAME  DISPLAY,  THE  LETTER  A 
REPRESENTS  GROUP  1,  B REPRESENTS  GROUP  2,  ETC.  FOR  2 POINTS 
CLOSE  TOGETHER  THE  LETTER  HILL  BE  PRINTED  WITH  AN  UNDERSCORE . 
THE  DISPLAY  CANNOT  HANDLE  3 OR  MORE  POINTS  CLOSE  TOGETHER  IF 
THERE  ARE  MULTIPLE  PLOTS. 


The  format  of  SCATHOW  is  followed  for  every  function.  The 
function  syntax,  the  user  controlled  parameters,  the  group  (or 
groups)  containing  the  functions,  the  subroutines  used  by  the 
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function,  and  a brief  description  of  the  use  of  the  function 
are  displayed. 

The  functions  and  variables  are  grouped  into  APL  groups 
for  ease  of  handling  and  to  conserve  space.  The  group  structure 
is  described  below. 

)GRPS 

COMP  A hi  SO  NS  DATA  ul  SPLAY  DOC  ESS  eh  rj  ALS  Gkh 

RELATIONS  SMOOTH  STATS  TaQN  AY 


)GRP  COUP  ART  SOU S 
CONDENSE  CON  DENSE HON 

SC  AT ri ON  LINE  LINEtiON  NDIVX 

COMPARE  COUP  AREHON  FILL 

AND  HON  LSLI.NE  LS  LI  HE  RON 


HUM  MISS  UTCOHO  t'MT  iC  AT 

NDIVY  HID  DEP  DEPT:  .GAP 

HUMS tJM  tiU\SUHdON  hShl-.CT  A •; D 


)GRP  DATA 

AGES  C ill  CENTS  COAL  CONSUMPTION  CRISES  u€  AT 

DISCOVERIES  DRAPERX  DRAPERY  GRADES  1 GRADES  2 HYDRO  PLANTS 

INSECTS  LANDAREAS  MORTALITY  PRECIPITATION  PRESRaTTNG 

RI V EhS  TtSTX  TESTY  USPOP  NARPbREARS  VOLCANO 


)GRP  uIFPlAY 

F u t. i<i  Le  A F F T Ed  L E A E tt C N 

bCXPLCTitCN  LENGTH  RILL 

CR 


)GRP  DOC 

AnDHON  ANOVAiiON  A3  R RON  A3  HSR  HUN  BJXPLOTnON 

CHIShUARRHON  CODERESHON  COHP  AhtriQ  V CONDENSE  lO* 

CONTINGENCY  HON  CORRELATION  HON  DAT  ASAP  DESCRIBE  KSiON 

LI N EU ON  LSLlNEdON  'AEDPG  LI SHtiUH  ytSTATRUN  EJ  iS  J i;iON 

ON  EH  RUN  UN  S3  H ON  PARTI  ALhON  REGRESSION  SC  AT  HON  SnG*  f.S  i *«  v 

SPLIT  HUN  ST  ATT  STICSHON  ST  Ed  LEAF  HON  TNJCStiOY 


NIU'Ih  SCALE  LIT  bCXPL'"T 

hLI'jFC  ! NUAFL'UnC* 
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)jtit  ESSE’iTI ALS 
A 3/i 


A u Li  Ad  00  A 

CUMP  Ah  t.  C 0 U c .7 6'  E 
OEPTti  aFSILOJ  FILL 
L J Li  ,7  £ kE  UPO LI S n 
dUilSU-i  U N Ed 
E it  Orth  t.  hi  o F LIT 
C.U.  hi  L TATS 


A 3 h S L bO A F L 0 T 

CONTI itu  F i ICY 
F VT  In  P-JT 

!‘tl  b S ii  LIZ  ii 

0 ii/3  3 OF  TIG  V FlF'GHESS 

ST EM Lk At  S J '-iS U 

FAh'TlAL  S'l  AXIS'! IC S 


C ' tIS K J Ah  a 
C 'j  h ne,  LA  T I ii  7 
as  Li 

/ i/.--  7 i < Oh  • ■.  i 
RESIOS  i>CAL:. 
Td ICt.  J TCi,  i 0 
R SELECT  LI  ’iTzi 


a 1 1 / •> 

7 r /i  L 
■v  I 
a FT 


L-  i 


)GRP  JFTT 

tiS  K.S  HU  0 C HISU'J  AhE 


C HIS^J  A hF.H 00 


) GkP  k F L AT I C ti S 
Lint.  LInF.ttCO  LSLIt.r. 

Li  LI  lt  EdC.t 

L F a !•  h S S 

f. 

/i  iiOtiC'O  S LF.CT  SC  Ah 

SC  AT  d CO  no 

UF.P 

Ul V 

E INTERCEPT  Fill 

CR 

i’ 


)GPP  SMOOT ti 
SPLIT  OHS  3 OiVift’ 

AlkHOO  A2kS  hHOO 
iU/IZY  OIO  OFF 


A 3 /i  4 3 PS  ft  TO  I C E SCAT 

ONE2UOO  OUKHUOO  TOICSdOO 


SPLIT uO / 
Sf/S^Y  ,Z 


)SAP  5T4:T5 

CGtihF.LATZCd  C h k'i-i‘7 

C C RnFLATICHHC « iST  A I S HCo 


hSTATS  Fa n II AL 
PARI I Aid  CO 


STATISTICS 

U "l  A'~  7 C 'r*e 

t i. 


)SftP  TOGO  AY 

ME  OF 0 LI  Sd  I4EDP0  LISHHOO 

RES IDS  EPSILON  .WRMEFFECTS 
COOERESdOO  AND  AHOdOO 

AUOZAdOO  OPTION  INPUT 


LINE  LTHEHOO  LSLIOE  LSLTfiSaOO 
SdOO RES  CODERES  kSLLECT  SdOOkES.iT 
CONTINGENCY  CONTTiPiE'tCY  100  A 

SU  !'iS  G 


ill  A 


'iOZ  A 
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The  user  can  conserve  space  in  his  active  work  area  by 
selectively  loading  or  copying  only  the  required  functions  or 
groups.  We  have  attempted  to  anticipate  the  type  of  analyses 
that  the  user  will  attempt,  and  the  group  structure  has  been 
selected  with  the  objective  to  combine  functions  and  parameters 
which  naturally  go  hand-in-hand.  A specific  group  can  be 
selected  by  typing 

) COPY  2 OA3660  GROUPNAME 

where  GROUPNAME  is  any  one  of  the  groups  listed  above.  Two  of 
the  groups,  ESSENTIALS  and  DOC,  are  especially  useful.  The 
ESSENTIALS  group  contains  all  of  the  functions  and  default 
values  of  all  of  the  required  parameters.  Documentation  and 
data  are  not  included.  The  DOC  group  contains  all  of  the 

documentation  (DESCRIBE,  DATAMAP , and  "HOW"  variables). 

The  user  may  want  to  load  the  entire  OA3660  workspace,  check 
on  some  of  the  documentation,  and  then  erase  the  group  DOC 
to  make  room  for  user  generated  data  and  other  functions  or 
variables.  The  DOC  group  can  be  erased  by  typing  ) ERASE  DOC. 
Similarly,  any  group  can  be  erased  by  typing  ) ERASE  GROUPNAME. 
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The  following  sections  give  brief  descriptions  (the 


writeups  contained  in  the  "HOW"  variables)  and  examples  of  the 
functions  contained  in  the  OA3660  library.  Each  function  descrip- 
tion gives  the  syntax,  a list  of  the  parameters,  a list  of  the 
groups  which  contain  the  function,  a list  of  subroutines  (other 
functions)  used  by  the  function,  and  a narrative  about  the  use 
of  the  function.  The  narrative  tells  what  the  function  does, 
describes  any  restrictions  on  the  function  arguments,  and 
indicates  how  the  function  parameters  affect  the  output.  The 
order  of  presentation  of  the  functions  follows  what  appears  to 
me  to  be  a natural  sequence  of  exposure  to  tools  of  data 
analysis  . 

The  HOW  variables  are  listed  alphabetically  in  Appendix  A, 
and  the  APL  program  listings  are  given  in  Appendix  B. 
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)GRp  DISPLAY 

SIExLEAF  ST EuL P AF  HOG  'll  Ilk  FT  ALL  L 

cCY.PLCI.iCS  Li.  nUln  FILL  xU  iFU.i  GUkSU 

rti 


Tr  hi/  AL  ) 


2 ) A v»  # v *.  j.  1 1.  t-\  i 

1 'i  i?  3 ( 0 h.  S A'J LI = 1 / 

4 SIR  Pi  A 1<D  LEAF 


it  i r r L h. 


txt.l-ihL  hi n>.  S 

SYuTAkt  ST Ek LEAF  VECTOR 
G ECU  Pi  DISPLAY 
P An  A iF.T t.  RS  : 

(1)  SI  LT  tt-  CCS? ROLF  Ini  mILTu  (~  H Ah  All  t.RS  Pei. 
THE  DISPLAY  ( Df E AU LI =7 0) . 

(2)  SCALE-  VARIES  Ini.  DEPTH  if 'IE  t Lt.'rRZAL) 
DISPLAY  lit  Util'IS  CP  1,  2,  i?  A 0.5  ri-.v  f .-1  f r.2r  A 
F ELECT  As  IkTiuRR  Ft.  OP,  1 TO  3 ( Dh.S  Ai'LT=  1 ) . 

SU  bPhCGR  AS : LI? 

LESC  RIPTICL  i FTEMLEAF  GENERATES  A FT  Rk  ALL  LEAP  LIFPL, 
LECTOR  OF  CPSFRhi/  ATICsS . THE  FUR  Cl  ICS  AUTOMATIC  ALLY  SC, 
DATA  USING  n SC  ALIiiG  ROUT I RE  dASlu  CL  I HE  kAti'JP.  Ai/D  SIDE 
LAI  A VECTOR . bY  VARYING  I HE  FC  ALE:  PARAkETEh  THE  USER  CA. 
THE  SIZE  OF  THE  FTi.k,  BUI  IT  SILL  ACT  sF.CEFFAP.ILY  L E Si 
Tab:  AkCUHI  SPECIFIED,  TRY  VARIOUS  VALUES  LIKE  3,2,  A.D 
SnlCtt  CHOICE  YIELDS  THE  BEST  DISPLAY.  IP  A LEAP  ■ "CpT. 

yi a L y characters  to  lp.  prihieu  On  a dI>/p  , she  li  . .*• 

Thl'SCATED  ALL  TiiE  xUkbEh  CP  TRUNCATED  CHARACTERS  SILL  of 
AT  THE  EHD  OF  TtiE  LIRE.  If  THE  FI  Eli  IsTF.RVAL  IS  ISC 
PC. sEh  OF  lEk , THE  LEAVES  i-iAY  CCiilAIti  IhF  CdARACItEF  4, 
In  AL LI I ICR  IC  THE  DIGITS  0,1,...,  9.  Is  THIS  Ca. 


’ U Ah  ACT E RS  A 


REPRESENT  Tiit.  sULLEhE 


SCALE 
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l i ^ 2 ^ j 3 b b 7 9 
3 3 j 1 2 3 3 4 4 b b 6 7 7 9 
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L 1 2 4 4 4 5 G 
3 3 4 7 3 


3 3 3 7 9 
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0 1 1 1 1 j 3 3 3 3 4 4 4 6 7 8 8 9 9 9 9 9 /’  4 & o^CC  J £ tf  jujjj.y; 
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3 1 1 *.  2 3 3 4 4 5 6 rf  9h4  bf  o Ur.  Gu  • 
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A -Til/',  i h.tr.  :iU.:L>r.  h , 
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1 U bp  liCj  ixAi-i 
j b IP  tiZ  P!  10 
Ob’  Uni  A » 
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PPL  AY 


5 0 0 n Ah  API  t til ) . 


t x 


i 
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NU 14 

SYNTAX t 
PA'RAMb  T Eh 


G i\  C U R 


c 


SU  UPROGRA 
DESCRIPT  I 
In  ALL  CA 
PRCMC  es 
DLL E 1 1 On 
SPREADS, 
USbtUL  SC 


SUMnON 
NUllSUGl 

i 

.■US  Sr, 

A tit.  AY 
C.-iP  A A ISO  it  S 
,4  ; Eh? 

Ci M NU  AS  U 14 


V ECTOR 

NUMBER 
( DEF  AU . 


■V U '4 S U A An  RAY 


US flo  TO  1,1  DU 
/I  = ”99999). 


■SING  VALUE  I. 


on : 


OPERATES  Ok  EITHER  A 


V EOT 0 R 0 A A 

EES,  THE  DATA  ARE  TREAT  Eb  AS  A SI  EG  Lb  okTOtl . 
A NUMERIC AL  SUMMARY  'A RICH  GIVES  THE  SAMPLE  SI 
OF  14 I ESI  nG  VALUES),  THE  EIGHTHS,  THE  f.KTh 
AND  THE  MIDPOINTS  IN  TABULAR  FORM.  THIS  SU 
r TESTING  THE  SYMMETRY  OF  a oATA  SET,  AND  TO 


The  effectiveness  OF  DATA  TRANSFORM  A:, 


i S IN  P n 0 DU 0 1 n u 


A n t a / o U 

■ t L it. 

Li  r ( At  J 

E-.bS , u 
>j  4 A R Y I 
EV ALU  A 
SYMMETRY 


\ t 


NUuSUtu  VOLCAnO 


NUMERIC  AL  SU  rh-i  A rY 


/ a i it L t.  t 1 L E = 

219 

lilUPTS 

LOG/ S/i;  In 

MEDIAN 

L'PQ/  8/ MAX. 

S Pl.t  AO S 

65.00  | 

1 

65.00  | 

1 

1 

66.50  | 

37.00  | 

1 

96.00  | 

59. 00  ! 

72.50  | 

2U.00  | 

! 

121.00  i 

97.00  ! 

100.50  | 

2,00  | 

1 

199.00  | 

2 j 7 . 0 0 ; 

fl  THE  ST E u Lr.  AF , THE  bOXPLOT , AnD  THE  uU-tEHIC AL  Sl'.i.;  AtY  noo 
n SUGGEST  THAT  THE  VOLCANO  DISTRI rUTIOn  IS  POSITIVELY 
n SKEGeD.  TrY  a LOG  TRANSFORM  TO  Err.  IF  THAT  FILL  -.An; 
n Tub  UISTRIDUTION  ,uOht.  •EARLY  SY  .ElE'T  KJC . 


hUmSUu  ® V 0 LC An 0 


NUnEnlCAL  SUnHARY 


SA.;P oe  Slob. 

= 219 

MI  DPT S 

LC(a/  / o / - ‘ii  v 

•lb  ul An 

U PI,  / 8/.’iWA' 

SPREADS 

4,17  | 
4.09  | 

i i 

3.61  | 

4.17  | 

1 

] 

4.56  | 

[ 

: . 9 5 ! 

3.99  | 

3.13  | 

1 

4.80  1 

1.62 

2,99 

1 0.69  | 

1 

5.29  i 

4,60 
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ft  THERE  IS  NOW  A DECREASING  TREND  IN  THE  MI DPT S COLUMN. 
ft  THAT  SUGGESTS  THAT  THE  LOG  IS  TOO  EXTREME  A TRANSFORMATION . 
ft  TRY  SQUARE  ROOT. 


NUMSUM  VOLCANO* . 5 


NUMERICAL  SUMMARY 
SAMPLE  SIZE  = 219 


MIDPTS 

LOQ/8/MIN 

MEDIAN 

UPQ/8/MAX 

S PR  EADS 

8.06  | 

1 

8.06  | 

1 

7.94  | 

6.08  | 

1 

9.80  | 

3.72 

7.95  | 

4.90  | 

1 

11.00  | 

6.10 

7.76  | 

1.41  | 

1 

14.11  | 

12.69 

ft  THAT  LOOKS  A LOT  BETTER.  LET  US  OBTAIN  THE  BOXPLOT 
ft  AND  STEMLEAF  PLOT  OF  THE  SQUARE  ROOT  TRANSFORMED  DATA. 


BOXPLOT  VOLCANO*. 5 

1.4 

14.1 


o O X X X X X X X X X X X X X X X X X X | * I X X X X X X X X X X X X X X X X X X OOC 


LENGTH*-  5 0 

BOXPLOT  VOLCANO*. 5 

1.4  14.1 

I I 

o O x x x x x x x x x x x x x I * I x x x x x x X x X X X x x oo 


ft  THE  BOXPLOT  LOOKS  A LOT  MORE  SYMMETRIC. 
ft  STEMLEAF  PLOTi 


NOW  FOR  THE ■ 


5 2 . L F h l?  / CL~ A:.1'?*  . 5 


1 


01  | u 

0 2 i 2 4 4 4 6 

0 3 I 0 0 2 2 2 3 5 

04 | 001445566777399 

0 5 [ 00111223344456677889999 

Ob! 000122 2 3 444 5 66666 6 668993 

0 7 | 0 0 0 0 0 1 1 1 1 2 2 2 3 3 3 4 4 5 8 5 5 5 6 -j  6 5 7 7 7 7 8 c 

08  10013111122  2 2234444  4 44  5 856  7'?  7 738c, 

0 9 1 Dll 111132  3 3 455556  6677778869 

1 0 ( 0 3 011112  2 <.344555666  7 889 
iij  0 J 0 1 1 1 2 2 4 5 6 7 7 y 8 
12155578 

1 3 | 146899 
1 4 j G 3 
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)GRP  COMPARISONS 


CONDENSE 

CONDENSEHOW 

nun: 

MISS  UTCCnd 

EM? 

SC  AT  HOW  LINE 

LIN  F.IIOW  It  DIVA 

NDIVY 

WID  DEP 

DEPTH 

COMPARE  COMPARE HOW  FILL 

ANDHOW  LS LINE  LSLINEHOW 

RUMS CM 

NVMSUhHOW 

REFLECT 

AH  DHOW 

SYNTAX\  X ADD  Y 

GROUPS : RELATIONS , TWOWAY,  COMPARISONS 

DESCRIPTION t THE  FUNCTION  'AND'  IS  USED  TO  CREATE  A <EW  DATA 
ARRAY  CONSISTING  OF  X AUGMENTED  BY  Y AS  ADDITIONAL  COLUMNS,  X 
AND  Y CAN  BE  SCALARS , VECTORS,  OR  MATRICES  BUT  CANNOT  tOT  > >. 

SCALARS,  IF  AN  ARGUMENT  IS  A SCALAR,  A COLUMN  IS  GENERATE!  ~»Ch 
ELEMENT  OF  WHICH  IS  THE  SCALAR.  IF  THE  ARGUMENTS  ARE  VECTORS  C: 
MATRICES,  THEY  MUST  he  conformable. 


a THE  FUNCTION  } AND ' IS  USEFUL  FOrDATA  ENTRY  AND  iA 
a TION . 


r 


'LA-' 


3 

u 

5 


3 

u 

5 

50 
50 
5 0 


50 
5 0 
5 0 


X *■  3 4 5 
Y- 9 10  1 2 
X AND  Y 
9 

10 
1 2 

Z*-3  2 p i 6 
X AnD  Z 
1 2 
3 4 

5 6 

50  AND  Z 
1 2 
3 4 

5 6 

5 0 AND  X 

3 

4 

5 


a EXCEPT  When  ONE  ARGUMENT  IS  a SCALAR,  THE  ARGU.-.r  NTS  Of 
a AND  mUST  bE  CONFORMABLE  IN  SIZE. 

W+- 4 6 8 1 0 
X AN u W 

ARGUMENTS  OF  AND  ARE  not  CONFORM AcLF. . 

W A N D Z 

AhGU:lENTS  OF  AND  ARE  NOT  CONFORMABLE . 
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COMPAREHOW 

SYNTAX i COMPARE  MATRIX 
PARAMETERS i 

(1)  DEPTH - VERTICAL  HEIGHT  OF  DISPLAY  {DEFAULT- 2 0 
Ll/VffS ) . 

(2)  NUMBER  USED  TO  INDICATE  MISSING  VALUES 
{DEFAULT  ="99999) 

(3)  N GAP  * NUMBER  OF  HORIZONTAL  SPACES  BETWEEN  THE 
BOXPLOT  DISPLAYS  {DEFAULT* 3). 

GROUP \ COMPARISONS 
SUBPROGRAM t FILL 

DESCRIPTION t rOMPARE  OPERATES  ON  AN  N BY  K MATRIX  TO  PRODUCE  K 
VERTICAL  BOX  PLOTS  PLACED  NEXT  TO  EACH  OTHER  TO  ALLOW  VISUAL 
COMPARISON  OF  THE  CENTERS,  SPREADS,  AND  OUTLIERS  OF  THE  BATCHES 
{COLUMNS)  OF  THE  MATRIX.  THE  USER  MUST  FILL  UP  THE  MATRIX  SO 
THAT  THERE  ARE  N OBSF.RV ATIONS  FOR  EACH  BATCH.  THE  PARAMETER, 

MISS,  SHOULD  BE  USED  TO  FILL  IN  MISSING  VALUES.  THE  USER  SHOULD 
ASSURE  THAT  MISS  IS  DIFFERENT  FROM  ALL  VALID  ENTRIES  IN  THE  DATA 
ARRAY . 


DEPTH 


20 

MISS 

"99999 

NGAP 


3 


a WE  WILL  DEMONSTRATE  COMPARE  WITH  THE  CHICKWT  DATA. 
a SINCE  WE  USED  0 TO  INDICATE  MISSING  VALUES  IN  THAT 
a DATA,  WE  MUST  FIRST  ALTER  THE  PARAMETER  MISS. 


MISS*-  0 
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COMPARE  CHICKWTS 

423 

108 


o 


x 


X 

1 


02 

1_  ~ T 

★ | x 

T~ 


1 

I I 
1*1 
I I 

1 1 
T 


O 


x 


1 

* 


T 


x 


a THE  DIAGRAM  REVEALS  RATHER  LARGE  DIFFERENCES  IN  THE. 
a MIDSPREADS  DF  THE  CHICKWT  DATA.  WE  SHOULD  TRANSFORM 
a THE  DATA  TO  TRY  TO  ACHIEVE  H0M0GF.NIF.TY  IN  THE  SPREADS. 
a FIRST  LETS  GET  A NUMERICAL  SUMMARY  OF  THE  DIFFERENT  COLUMNS. 


C0NDENSEH0W 

SYNTAX x CONDENSE  MATRIX  OR  R-CONDENSE  MATRIX 
PARAMETERS  x 

(1)  NL'Mm  CONTROLS  WHAT  STATISTICS  ARE  INCLUDED  IN  THr 
SUMMARY.  NUMcL  1,2,5,71  {DEFAULT*  2). 

NUM* 1 GIVES  MEDIANS  OF  EACH  COLUMN  OF  THE  MATRIX 
HUM* 2 GIVES  MEDIANS  AND  INTERQUARTILE  RANGES 
NUM* 5 GIVES  MIN .QUART I LES .MAX , AND  SIZE 
NUM* 7 GIVES  MIN .EIGHTHS .MAX  AND  SIZE 

(2)  MISS r NUMBER  USED  TO  CODE  MISSING  V ALU E S 

(Z?£F/J£,L7  = ~99999). 

GROUP t COMPARISONS 
SUBPROGRAMS x UTCOND.  FMT 

DESCRIPTION i CONDENSE  GENERATES  SUMMARY  STATISTICS  FOR  EACH 
GROUP  OF  DATA  REPRESENTED  BY  THE  COLUMNS  OF  THE  ARGUMENT  MATRIX. 
SINCE  THE  ARGUMENT  IS  A MATRIX.  THE  USER  MUST  FILL  UP  THE  nlSSIuU 
VALUES  IN  THE  MATRIX  WHENEVER  THE  GROUPS  ( COLUMNS ) HAVE  DIFFERENT 
NUMBERS  OF  OBSERVATIONS.  THE  PARAMETER,  MISS,  SHOULD  BE  USED  TO 
FILL  THE  MATRIX  SINCE  THE  FUNCTION  WILL  RECOGNIZE  THOSF  VALUES  AO 
MISSING  DATA  AND  WILL  IGNORE  THEM  IN  ALL  CALCULATIONS.  THE 

RESULT  FROM  CONDENSE  CAN  BE  USED  AS  THE  ARGUMENT  OF  OTnr.R 
FUNCTIONS  SUCH  AS  SCAT  OR  LINE  IF  AND  ONLY  IF  NUM*2. 
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2 


HUM 


MISS 


CONDENSE 

CHICKWTS 

151.5 

43 

221 

83.5 

248 

72 

328 

33 

263 

83 

342 

102 

A THE  ABOVE  DATA  ARE  THE  MEDIANS  AND  MIDSPREADS  FOR  EACH 
a OF  THE  6 BATCHES  ( COLUMNS ) OF  CHICKWT  DATA . ZERO  IAS 
a ENTERED  INTO  THE  PARAMETER  MISS  EARLIER . THE  MIDSPREADS 
a SHOW  QUITE  A DISPERSION , Fit3M  -45  SMALL  AS  151.5  TO  AS  LA  ROE 
A AS  342. 


A WE  CAN  GET  MORE  INFORMATION  ABOUT  THE  INDIVIDUAL  COLUMNS 
a OF  CHICKWTS  BY  CHANGING  HUM  AS  FOLLOWS i 


NUM+-5 

CONDENSE  CHICKWTS 


MIN 

LOQ 

MEDIAN 

UPQ 

MAX 

SIZE 

108.00 

136.00 

151.50 

179.00 

227.00 

10.00 

141.00 

175.00 

221.00 

258.50 

309.00 

12.00 

158.00 

199.00 

248.00 

271.00 

329.00 

14.00 

226.00 

307.50 

328.00 

340.50 

U23.00 

12.00 

153.00 

242.00 

263.00 

325.00 

380.00 

11.00 

216.00 

271.50 

342.00 

373.50 

404.00 

12.00 

A I N 

SEEKING  OUT 

AN  APPRO 

PRIATE  TRANSF 

ORE  AT  10 

N TO  EVEN 

a THE 

MIDSPREADS , 

IT  IS  0 

Fr;EN  USEFUL  T 

0 PLOT 

THE  LOG  0 

a MIDSPREADS  VS.  THE  LOG  OF  THE  MEDIANS. 

a IN  MANY  CASES  WE  NEED  TO  LOOK  AT  A SrATTEk  PLOT  OF  OUR 
a DATA,  OR  AT  SCATTER  PLOTS  OF  FUNCTIONS  OF  OUR  DATA. 
a LET  US  INTRODUCE  A FUNCTION  SCAT  THAT  WILL  PROVIDE  THE 
a NEEDED  SCATTER  PLOT. 
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SCATHOW 

SYNTAX i SCAT  ARRAY 
PARAMETERS t 


(1)  WID*  CONTROLS  THE  HORIZO 
{DEFAULT- 30  CHARACTERS) . 

NT  A L SIZE 

OF 

THE 

DISPLAY 

(2)  DEP*  CONTROLS  THE  VERTICAL  SIZE 
(DEFAULT*!*  LINES). 

OF 

the 

DISPLAY 

(3)  NDIVX,  NDIVY - NUMBER  OF 
RESPECTIVELY  ( DEF AU LT=H ,4) . 

UNITS  ON 

X - 

AND 

Y-AXES, 

GROUP x RELATIONS,  SMOOTH,  COMPARISONS 

DESCRIPTION  t SCAT  PRODUCES  A SCATTER  PLOT  OF  THE  DATA  CONTAINED 
IN  ARRAY.  THE  ARGUMENT  ARRAY  CAN  BE  A VECTOR  OR  A MATRIX  WITH  AS 
MANY  AS  9 COLUMNS.  IF  A VECTOR  OF  SIZE  N,  THOSE  VALUES  ARE 
PLOTTED  VS.  THE  INTEGERS  1 TO  N i IF  A MATRIX,  THE  SECOND,  THIRD, 
ETC.  COLUMNS  ARE  PLOTTED  VS.  COLUMN  1 ON  THE  SAME  AXES. 
DIVISIONS  ON  THE  AXES  OF  THE  PLOT  ARE  NOT  EXPLICITLY  PRINTED, 
EXCEPT  AT  THE  EXTREMES  OF  THE  PLOT.  THE  USER  CAN  CONTROL  THE 
RESOLUTION  OF  THE  PLOT  BY  MODIFYING  THE  PARAMETERS  WID,  DEP. 
NDIVX  AND  NDIVY.  PRINTING  TIME  INCREASES  DRAMATICALLY  WITH 
RESOLUTION ; THEREFORE,  OTHER  PLOT  PROGRAMS  SHOULD  BE  USED  IF  HIGH 
RESOLUTION  IS  DESIRED.  FOR  A SINGLE  GROUP  OF  DATA,  THE  NUMBER  N 
(2siVS9)  WILL  BE  PRINTED  IF  N POINTS  LIE  CLOSE  TOGETHER  ON  THE 
DISPLAY.  FOR  MULTIPLE  PLOTS  ON  THE  SAME  DISPLAY,  THE  LETTER  A 
REPRESENTS  GROUP  1,  B REPRESENTS  GROUP  2,  ETC.  FOR  2 POINTS 
CLOSE  TOGETHER  THE  LETTER  WILL  BE  PRINTED  WITH  AN  UNDERSCORE . 
THE  DISPLAY  CANNOT  HANDLE  3 OR  MORE  POINTS  CLOSE  TOGETHER  IF 
THERE  ARE  MULTIPLE  PLOTS. 


30 

15 

u 

4 


WID 

DEP 

NDIVX 

NDIVY 


X+-\ld 

Y 1*2  3 5 6 8 10  11  13  15  18 
72^20  20  17  24  21  26  31  30  37  33 
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a BACK  TO  THE  PROBLEM  OF  NONHOMCGENEOUS  MIDSPREADS  IN  THE 
a CHICKWT  DATA.  LET  L'S  LOOK  AT  A SCATTER  PLOT  OF  THE 
a LOG  MIDSPREADS  V.  LOG  MEDIANS. 

A FIRST  WE  RESET  PARAMETER  HUM  TO  2,  SO  IT  WILL  O'l/F 
A MEDIANS  AND  MIDSPREADS . 


NUM+  2 


SCAT  • CONDENSE  CHICKVTS 
RANGE  OF  Xi  5 6 

RANGE  OF  Yt  3.4  4.8 


O 

O 


a THE  SCATTER  PLOT  SUGGESTS  A RELATIONSHIP  BETWEEN  i-.IDSPRR 
A AND  MEDIANS.  PERHAPS  SOME  ROOT  OF  THE  DATA  WOULD  MAKE 
a THE  MIDSPREADS  MORE  HOMOGENEOUS . TRY  SQUARE  ROOT. 


CONDENSE  CHICKWTS* . 5 


12. 30368569 
14.86363273 
15.74801575 
18.11001266 
16.21727474 
18.4912651 


1.717184371 

2.851055496 

2.355341653 

0.9195158937 

2.471407191 

2.852064545 


A THAT  IS  BETTER  ( RATIO  OF  MAX  MIDSPREAD  TO  MIN  MIDSPhE AD 
A IS  ABOUT  1.5.  MORE  WORK  NEEDS  TO  BE  DONE  THOUGH. 


• i 


)GRP  RELATIONS 

LINE  LINEHOU  LSLINE  LSLINEHOW 
ANDHOW  RSELF.CT  SCAT  SCATHCW  WID 
A INTERCEPT  FMT  CR 


REGRESS  REGRESSHO'-J 
DEP  NDIVX  NDIVY 


LINF.HON 

SYNTAX t LINE  ARRAY  CR  Z+LINE  ARRAY 

PARAMETER t 

RSELECT"  SELECTS  THE  OUTPUT 
RSELF.CT  =1  GIVES  ABSCISSA  VS. 
RSELECT* 1 GIVES  FITTED  VALUES 
GROUP i RELATIONS 


MATRIX  : 
RESIDUALS 


F RESIDUAL.' 
( DEFAULT-1 ] 


FITTED  VALUES 


LINE  FITS  A STRAIGHT  LINE 


DESCRIPTION i LINE  FITS  A STRAIGHT  LINE  TO  A SET  OF  ( X,Y ) PCI: 
BY  DIVIDING  THE  POINTS  INTO  3 REGIONS  AND  USInG  THF  MEDIANS 
THE  X AND  Y VALUES  IN  THE  OUTER  REGIONS  TO  DETERMINE  THE  SLOPE. 
THE  INTERCEPT  IS  THE  MEDIAN  OF  THF.  DIFFERENCES  YtSLOPE*X.  L J 
ALSO  COMPUTES  THE  RESIDUALS  AND  GIVES  AN  N BY  2 MATRIX  CONTAIN. 
EITHER  THE  X*VALUES  VS.  RESIDUALS  OR  THE  FITTF.u  VALUES  l 
RESIDUALS  DEPENDING  ON  THE  PARAMETER,  RSELECT.  IF  THE  USER  Di 


SHE  SYNTAX  t 


COLUMN  MATRIX 
X "VARIABLE 


NOT  WANT  RESIDUALS  TYPED  OUT,  HE  SHOULD  USE  THE  SYNTAX t 

Z*-LINE  ARRAY 

THE  ARGUMENT  ARRAY  CAN  BE  EITHER  A VECTOR  OR  A 2 "COLUMN  MATRIX 
IF  THE  ARGUMENT  IS  A VECTOR  OF  SIZE  N,  THE  Xr .VARIABLE 
CONSTRUCTED  TO  BE  THE  FIRST  N POSITIVE  INTEGERS.  IF  THE  ARGUN 
IS  A 2n COLUMN  MATRIX,  THE  FIRST  COLUMN  IS  TAKEN  TO  BE  THE  SET 
X "VALUES  AND  THE  SECOND  COLUMN  THE  SET  OF  Y "VALUES.  LINE  C A\ 
USED  IN  CONJUNCTION  NITH  THF.  FUNCTION  SCAT  TO  PRODUCE  THE  SLO 
THE  INTERCEPT,  AND  A PLOT  OF  THE  RESIDUALS.  TO  DO  THIS,  ENTER 

SCAT  LINE  ARRAY 


16 


X+2  34356767598  10  11  12  14 
PX 


ARRAYS  AND  Y 
LINE  ARRAY 

SLOPE:  1.375  Y •INTERCEPT : 0.0625 


2 

0.1875 

3 

"0.1875 

4 

"1.5625 

3 

0.8125 

5 

1.0625 

6 

"2.3125 

7 

"5.6875 

6 

"0.3125 

7 

"0.6875 

5 

3.0625 

9 

"0.4375 

8 

3 . 9 3 7 5 

10 

2.1875 

11 

"1.1875 

12 

1.4375 

14 

o . 6875 

RSELECT+Q 
LINE  ARRAY 

SLOPE:  1.375  Y *1 NTERCEPT : 0.0625 


2.8125 

0.1875 

4.1875 

"0.1875 

5.5625 

"1.5625 

4.1875 

0.8125 

6.9375 

1.0625 

8.3125 

"2.3125 

9.6875 

"5.6875 

8 .3125 

"0.3125 

9.6875 

"0.6875 

6. 9375 

3.0625 

12.4375 

" 0 . 4 3 7 5 

11.0625 

3.9375 

13. 8125 

2 . 1 8 7 5 

15.1875 

"l. 1875 

16.5625 

1.4375 

19.3125 

0.6875 

RSELECT+ 1 

SC/IT  LINE  ARRAY 

SLOPE:  1.375  Y *1  NT ERCEPT : 0.0625 
RANGE  OF  X:  2 14 

RANGE  OF  Y:  "6  4 


O 
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Y 

3 4 4 5 8 6 4 8 9 10  12  15  16  14 

PY 

16 

LINk  Y 

5 LULL : 1.090909091  Y ■»!  IITkihCEPT  : 0.6818181818 


1 

1 . 17n  7n  7?  7n 

2 

1.136363636 

3 

0.04545454545 

4 

"0. 04545454545 

5 

1.863636364 

6 

" 1 . 2 2 t 2 7 2 17  2 17 

~4. 318181818 

8 

"1.409090909 

9 

"1.5 

10 

"l. 590909091 

11 

"0.6818131818 

12 

1.22n2'7212n 

13 

1.136363636 

14 

"1.954545455 

1 5 

0. 9 5454  54545 

16 

1 .863636364 

*«-i  16 

ARkAY-X  An D Y 
LI  ilk  ARRAY 


SLOPl:  1.090909091  Y -I  I1TF  RC  kPT : 0.6818131913 


1 

i. 

1 . 2 2 " 2 " 2 2 

2 

1.136363636 

3 

0.04545454545 

4 

"0. 04545454545 

5 

1 .8  63636364 

6 

"1 . 22'72',2'72'’ 

7 

"4.318131818 

8 

"l  . 409090909 

9 

"1.5 

1 0 

"l . 590909091 

1 1 

"0.6818181818 

12 

1 . 2272"2'72'’ 

13 

1.136363636 

14 

"1.954545455 

1 5 

0. 9545454545 

16 

1 .8  63636  364 

SCAT  Link  Y 

SLOPE  : 1.090909091  Y -I v TkhCkPT : 0.6818181818 

RAUGb  Ok  X:  0 20 

RAUGk  OF  Y:  ~ 6 2 


O 


o 


o 


o 


o 


o 


o o o 


o 


o 


o 
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RFF.LECT-*-! 

Z-LINE  Y 2 

FLO  PE  t 1.857142957  Y -INTERCEPT  s 17.42857143 


LFLINEHCN 

FI *2 AX  : LFLInE  ARRAY  OR  Z+LFLIuE  ARRAY 

GROUP  i RELATIONS 

DR'FO RIPTICN  i IF II RE  LET  F.Ri-iI  N EF  FRF.  LF AFT  SQUARES  SOLUTION, 

Y = YI+FLQPE*X . 17//?  ARGUMENT  ARRAY  OAR  BF  A VECTOR  Oh  A 2 -00 

MATRIX.  IF  THE  ARGUMENT  IF  A VECTOR,  THE  X-V A LUFF  ARE  TAKE*  TO 
BE  THE  INTEGER  F 1 rC*  iV  ►/FF/rF  A/  = p/IFF/SY . I/*’  THE  ARGUTE.!  IF  A 

2-CCLUMt  i/  MATRIX , COLUMN  1 CONSISTS  OF  THE  X-V  ALU  FF  A HD  CC'LU',.u  2 
THE  Y "VALUES.  REFILL'  ALF  ARE  DETERMINED  AND  FT  OREL  IN  A ,AT;:Ia 
NhCSE  FIRFT  COLUMN  IF  THE  FFT  OF  X -VALUE?  ALL  InOFF  F E C 0 n D CCLVW 
IF  THE  FFT  OF  REFILL' ALF.  IF  THE  USER  nANTF  TO  SUPPRESS  Pul,.  OUT 
OF  THE  REF  I DU  ALF,  Hi  MU  FT  ASSIGN  THE  RF.FULTF  Ot  LFLInE  TO  /. 
VARIABLE',  I.E. , TYPE  Z*-LFLINE  ARRAY . 


SCAT  LFLINE  X AND  Y2 

FLOPEi  1.933333333  Y -INTERCEPT  : 16.8  6 666667 
u'AnGE  C F X t 0 10 

n /3  :v  j F CF  Y:  ~ 4 a 


nIF  raft'  05FA  'HANTS  TO  h F GREFF  A REF  PON  FF  VARIABLE  On 
n / 5F?  OF  C ARRIEhF , HE  CAN  UFE  THE  PROGRAM  RE  GREFF . 
n T h IF  LEAF T FQUAREF  MULTIPLE  REGRESS I ON  PROGRAM  IF 
n /FAY  POWERFUL  AND  FLEXIBLE. 
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REGRESSHOW 

SYNTAX i Z'-r  REGRESS  X 
PARAMETER i 

A INTERCEPT*  DETERMINES  WHETHER  CR  NOT  AN  INTERCEPT  TER’l 
IS  TO  BE  INCLUDED.  NINTERCEPT-1  GIVES  AN  INTERCEPT 
TERM , AND  ^INTERCEPTED  GIVES  NC  INTERCEPT.  {DEFAULT  IS 
1.) 

GROUP \ RELATIONS 
SUBPROGRAMS i FMT  AND  SCAT 

DESCRIPTION i REGRESS  DOES  A MULTIPLE  REGRESSION  ANALYSIS 
RELATING  THE  DEPENDENT  VARIABLE  Y TO  A SET  OF  CARRIERS  X.  THE 
LEFT  ARGUMENT  Y IS  A VECTOR  OF  SIZE  N.  THE  RIGHT  ARGUMENT  X IS 
AN  N BY  K MATRIX  CONSISTING  OF  N OBSERVATIONS  ON  EACH  OF  K 
VARIABLES  OR  A VECTOR  OF  SIZE  N IF  K=l.  OUTPUT  CONSISTS  OF  Ah 
ANCVA  TABLE , RnSQU  ARE , STD.  ERROR,  REGRESSION  COEFFICIENTS  {THE 
FIRST  COEFFICIENT  IS  THE  CONSTANT  TERM  IF  LINT ERC F.PT -1 . ) , 
TvSTATISTICS , V ARI ANCEnCCV  ARI ANCE  MATRIX,  DURB IN -W AT SON 

STATISTIC,  AND  A VECTOR  OF  PREDICTED  Y VALUES  AND  RESIDUALS. 
THERE  IS  An  OPTION  THAT  ALLOWS  THE  USER  TO  INPUT  A VECTOR  OF  X 
VALUES  AND  USE  THE  REGRESSION  EQUATION  TO  FORECAST  Y VALUES.  THE 
USER  CAN  ALSO  OBTAIN  A SCATTER  PLOT  OF  THE  RESIDUALS.  NhEN 
EXECUTION  TERMINATES , THE  PREDICTED  Y VALUES  AND  THE  RESIDUALS 
RESIDE  IN  THE  N BY  2 MATRIX  Z. 

NINTERCF.PT 

1 

p DRAPERX 

13  4 

*LETS  REGRESS  a RESPONSE  VARIABLE  DRAPERY  ON  THE  FOUR 
fl  CARRIERS  IN  THE  MATRIX  DRAPERX. 


61 


Z+DR APERY  REGRESS  DRAPERX 


ANCVA 


SOURCE 

DF 

SUM  SQUARES 

MEAN  SQUARE 

F "RATIO 

REGRESSION 

4 

2 . 667  9E3 

6.6697E2 

1 .1148? 2 

RESIDUAL 

8 

4.786  4? 1 

5.9830E0 

TOTAL 

12 

2 . 7 1 5 8?  3 

R SQUARE i 0.9823756203 
STD  ERROR t 2.446007963 

COEFFICIENTS  T STATISTICS 
62.4054  0.8906 

1.5511  2.0827 

0.5102  0.7049 

_ 0 . 1 0 1 9 0.135 

“0.1441  “0.2032 

DC  YOU  WANT  A PRINTOUT  OF  THE  VARIANCE "COVARIANCE  MATRIX? 

Y 

V ARIANCE*COV ARIANCE  MATRIX i 

4. 9099fc’t-3  “5.0507E+1  “5.0603E+1  “ 5 . 1 6 6 0?*1  4.9  597''’*! 

"5.0507F+1  5 . 5 4 6 8?“ 1 5.1266?“l  5.5425?  1 5.0529 ?“l 

“5.0603E+1  5.126  6?~1  5.2387?“l  5.2570?“l  5.1213?“l 

“5.1660E+1  5 . 5 4 2 5?“l  5.2570?“l  5.6959?“l  5.1688?“l 

“4.9597?*1  5.05  29?“l  5.1213?“l  5.1688?“l  5.0  27  5?*! 

DU  RBI N wW  AT  SON i 2.052596933 
DC  YOU  WANT  TO  FORECAST  A VALUE  FOR  Y? 

20  30  15  40 

DC  YOU  WANT  TO  SCAT  RESIDUALS  VS.  PREDICTED  Y? 

Y 

RANGE  OF  Xx  70  120 

RANGE  OF  Yi  “4  4 
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SPLIT 

) GRP  SMOOTH 

ONE  3 CNF.H 

A3R 

A 3 RSR 

TWICE 

A3RH0W 

A3RSRH0W 

CNE3H0W 

CNEHHOW 

TWICE HOW 

NDIVY 

WID  DEP 

SPLITHO'W 
SCATH ON  frWY* 


0NE3HCW 

SYNTAX t CUE 3 VECTOR 
GROUP t SMOOTH 

DESCRIPTION i CNE3  SMOOTHS  A VECTOR  OF  DATA  USING 
RUNNING  MEDIANS  OF  3 WITH  TU KEY'S  ENDrPOINT  RULE . 

U,Y ) PAIRS , THE  USER  SHOULD  USE  AS  THE  ARGUMENT 
VALUES  ORDERED  ACCORDING  TO  THE  MAGNITUDE  OF  THE  X 

VECTCR+YlkXl 

IT  IS  ASSUMED  THAT  THE  X VALUES  ARE  EQUISPACED. 
SMOOTHING  8Y  MEDIANS  OF  3 CAN  BE  ACCOMPLISHED  WITH  A3R. 
CONSISTS  OF  THE  SMOOTHED  SEQUENCE.  A PLOT  OF  THE 
SEQUENCE  IS  GIVEN  BY  TYPING i 

SCAT  ONE 3 VECTOR 


ONE  PA 
WHEN  SMC 
VECTOR 
VALUES j 


55  OF 
0 THING 
THE  Y 


RE 


PEA T E D 
OUTPUT 
COT  RED 


X 

3237432  15  4536 

ONE  3 X 

233443345453 


SPLIT  HOW 

SYNTAX i SPLIT  A3R  VECTOR 
GROUP i SMOOTH 

DESCRIPTION  i SPLIT  DOES  CNF.  PASS  AT  DIVIDING  MESAS  ( PAIRS  OF 
ADJACENT  POINTS  WITH  A COMMON  VALUE  WHICH  IS  A LOCAL  MAX  OR  MIN) 
USING  Tl’ KEY'S  END-POINT  RULE.  IT  IS  USED  IN  CON 'UNCTION  WITH 
A3R. 


SPLIT  ONE  3 X 

233334345453 


A3RHCW 

SYNTAX i A 3R  VECTOR 
GROUP i SMOOTH 
SUBPROGRAM < ONE  3 

DESCRIPTION i A3R  DOES  REPEATED  SMOOTHINGS  OF  RUNNING  MEDIANS  OF  3 
UNTIL  THERE  ARE  NO  CHANGES  IN  THE  SMOOTHED  SEQUENCE  FROM  ONE 
ITERATION  TO  THE  NEXT.  SEE  CNE3HOW . 
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A3R  X 

333443344444 


A 3RSRH0W 

SYNTAX i A3RSR  VECTOR 
GROUP i SMOOTH 

SUBPROGRAMS  i /«3tf,  CiVf3,  SPLIT 

DESCRIPTION  i 3hSR  DOES  REPEATED  SMOOTHINGS  BY  RUNNING  MEDIANS 
OF  3 FOLLOWED  BY  SPLITTING  MESAS  AND  REPEATING  UNTIL  CONVERGENCE . 
SEE  A3R  AN L SPLIT. 


A3RSR  A 

333334444444 


ONER  HOW 

SYNTAX i CNEH  VECTOR 
GROUP « ' SMOOTH 

DESCRIPTION  I C/Vfrf  4 SINGLE  HANNING  OF  A SEQUENCE  OF  DATA. 

IT  IS  USED  IN  CON •UNCTION  WITH  THE  OTHER  SMOOTHING  FUNCTIONS  AS  A 
FINAL  TOUCH wl'P  TO  A SEQUENCE  OF  SMOOTHED  DATA.  THE  IiTH  RESPONSE 
IS  REPLACED  BY  0 . 2 5*7 [1*1 ] + 0.5*7[I1  ♦ 0. 25*7[J+1] . SEE  A3R  AND 
A3RSR. 


ONF.H  X 

3 2.5  3.75  5.25  4.5  3 5.5  9 7 4.25  4.25  6 


TWICEHOW 

SYNTAX i TWICE  VECTOR 
GROUP i SMOOTH 

SUBPROGRAMS i CNE3,  SPLIT,  A3RSR 

DESCRIPTION i TWICE  SMOOTHS  A SEQUENCE  OF  DATA,  THEN  SMOOTHS  THE 
RESIDUALS  AND  ADDS  THE  SMOOTHED  RESIDUALS  BACK  TO  THE  SMOOTHED 
DATA  TO  OBTAIN  THE  FINAL  SMOOTH.  THE  SMOOTHING  IS  DONE  USIuG 
A3RSR. 
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£W  ICE  X 

33333444444U 


ft.VCV  Lb! 2 US  LOCK  A2  SOME  PLOTS  212  u c.?,.  / ! ^RF. 
n DATA  2 H AT  SiiORS  TtiF.  EFFb.OTS  OF  SvCOT rlllvG . 
flV/t'  V ILL  USE  Tub!  PRES  RATING  DATA. 


2IU*-R  3 
U E P ■*- 2 0 

.'T/5?  PRESkATZhG 
itA.iGr.  Ob'  X i 3 12  0 

RAuGF.  OF  Y\  0 100 


2 ° 2 ° 2 ° o 2 2 ° ° 0 

OO  O 0 / <L  ° 2 < 

O 2°  3 ° 2 0 3 ooo  OO  o o o 


O/OOO  O / o 


O O /> 


PIT  IS  HARD  TO  TEll  il'Ou  t'hOk  THIS. 
p HELPS  ARY. 


t l‘  /"  c err  r r o * . /■>  / • — : 7 
r 1 ».  *.  7 J i.  r . ..  f > * .i  i / 


/5  3/(5a  PREShAT  IhG 
r.AiiGr  Ob'  At  0 1 2 0 

t, of  y i 2 o so 


<;  3 3 ° o 22  ° ° 2 3 ° 


2^3  .4  3° 


o 2 3 3 0 
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$ 

& 


*7  HAT  IF  A LOT  FliOOTHER . THE  PFAKF  An  L)  V ALLEY  F ARE 
a CLEARER  HOW.  AF  A FINAL  TOUT  HUP,  LET  F UFE  Tv  ICE.  ; 
a F14CCTH  THE  REFIDUALF  AN 0 ALL  LACK  TO  THE  F.J'CTHFl, 
a Ail  L THEN  UFE  ON  EH  AF  A FINAL  PCLIFH . 


U L 


FCAT 

ONER 

TWICE 

P REF RATING 

n An  G E 

OF  X. 

: 0 

123 

k An  Crt- 

! 

of  y 

i 2 0 

93 

1 

0 2 

2 

I ° 

23 

2 2 

1 

233 

O O 

2 0 o o 2 3 ° 

| O 

o o 

3 o 2 ° 2 ° 2 

I 

2 o 

O O 

1 2 22 

3 o 

O 

o 2 

CM 

o 

o 

0 

o 

2 

° ° ° 3 2 2 

! 

o o o 

233 

1 

2 2 

O 

0 

! 

o 

O 

1 

o 

0 

1 

[ t r*  *"•  ** 

0 

2 ° 
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) GRP  TWOWAY 

MEDPOLISH  MEDPOLISHHCW 

RESIDS  EPSILON  NCRMEFFECTS 
CCDERESHCN  AND  AN  DHOW 

ANOVAHOW  OPTION  INPUT 


LINE  LINEHOW  LSLINE  LS LINE HOW 
SHOWRES  CODERES  REEL  EOT  SHOWRESriOW 
OONTINGENOY  CONTINGENT!  HOW  A NOV  A 

SUMSQ 


MEDPOLISHHCW 

SYNTAX  i MEDPOLISH  MATRIX  OR  Z+MEDPCLISH  MATRIX 

PARAMETERS t 

(1)  RESIDS - CONTROLS  OUTPUT  OF  RESIDUALS.  IF  1, 

RESIDUALS  ARE  PRODUCED  IN  A TWO-WAY  TABLE ; IF  2, 

COMPARISON  VALUES , ( RE ° . *0  E ) AND  RESIDUALS  ARL 

PRODUCED  IN  A TWO  COLUMN  MATRIX  ( DEFAULT- 1 ) . 

(2)  EPSILON • PROPORTION  BY  WHICH  THE  SUM  OF  THE 
ABSOLUTE  VALUES  OF  THE  RESIDUALS  MUST  BE  REDUCED  AT 
EACH  ITERATION  TO  CONTINUE  POLISHING  {DEFAULT'D .01) . 

(3)  N CRMEFFECT S i DETERMINES  WHETHER  OR  NOT  NORMALIZED 
EFFECTS  ARE  OUTPUT.  THE  DEFAULT  OF  0 SUPPRESSES  OUTPUT 
OF  NORMALIZED  VALUES.  NORMEFFECTS-1  CAUSES  NORMALIZED 
VALUES  TO  BE  PRINTED. 

GROUP i TWOWAY 

DESCRIPTION i MEDPOLISH  ITERATIVELY  SWEEPS  OUT  MEDIANS  FRO  . THE 
ROWS  AND  COLUMNS  OF  A TWO-WAY  TABLE  TO  YIELD  THE  MODEL: 

CBS  = MEDIAN  + ROW  EFFECT  * COLUMN  EFFECT  ♦ RESIDUAL 
IT  YIELDS  RESIDUALS  FOR  TESTING  THE  ADEQUACY  OF  THE  MODEL  (SFr. 
SHOWRES  AND  CODERES).  THE  FUNCTION  CONTINUES  POLISHING  UNTIL  I Hr. 
STOPPING  RULE  CONTROLLED  BY  EPSILON  IS  ACTIVATED.  THE  OVERALL 
MEDIAN.  THE  ROW*  AND  COLUMN -EFFECTS , AND  THE  SUM  OF  ABSOLUT r 
VALUES  OF  THE  RESIDUALS  ARE  GIVEN.  THE  USER  CAN  SUPPRESS 
RESIDUAL  PRINTOUT  BY  ASSIGNING  THE  FUNCTION  TO  A VARIABLE,  I.E., 
BY  WRITING i 

Z*-  MEDPOLISH  MATRIX 


1 

0.31 

0 


RESIDS 

EPSILON 

NORMEFFECTS 
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DEATHS 


11.7 

8.7 

15.4 

8.4 

18.1 

11.7 

24.3 

13.6 

26.9 

20.3 

37 

19.3 

41 

30.9 

54.6 

35.1 

66 

54.  3 

71.1 

50 

Z+MEDPOLISH  DEATHS 

315.4 

40.3 

38.25 

37.575 

TYPICAL  VALl'Ex  25.6 

ROW  EFFECTSt  ”15.3625  ”10.1  “l.4  12.6375  32.875 

COLUMN  EFFECTSt  2.7  ”3.9  12.625  ”3.1375 


ftTHE  RESIDUALS  ARE  STORED  IN  Z.  THE  FIRST  VALUES 
n PRINTED  OUT  ARE  THE  SUMS  OF  THE  ABSOLUTE  VALUES 
a OF  THE  RESIDUALS.  ITERATION  CEASES  WHEN  TWO 
ft  .CONSECUTIVE  VALUES  ARE  WITHIN  EPSILON . 


Z 

1.2375  2.3625  "7.4625  1.3 

0.1  0.1  "3.825  1.2375 

0 0 0.175  "1.7625 

0.0625  "3.4375  3.7375  0 

4.825  "0.275  0 "5.3375 


♦/■'■/  I 2 

37.2375 


SHOW  RESHCW 

SYNTAX  i SHCWRES  MATRIX  OR  SHOWRES  MEDPCLISH  MATRIX 

SUBPROGRAM t AND 
GROUP i TWOWAY 

DESCRIPTION i SHOWRES  AIDS  IN  THE  ANALYSIS  OF  RESIDUALS  FROM 
MEDIAN  POLISHING  A TWOmWAY  TABLE.  THE  USER  SHOULD  SET  RESIDS= 1 
IN  MEDPOLISH  IF  SHOWRES  IS  TO  BE  USED.  THE  OUTPUT  OF  SHOWRES  IS 
A MATRIX  CONSISTING  OF  THE  SYMBOLS  Co°*X  DEPENDING  ON  THE  SIZES 
OF  THE  RESIDUALS.  LET  R REPRESENT  THE  RESIDUAL , M THE  MIDSPREAD 
OF  ALL  RESIDUALS,  L THE  LOWER  QUARTILE , AND  U THE  UPPER  QUARTILE. 
THE  CODES  ARE  PRINTED  AS  FOLLOWS i 
0 i R<Ln 1 . 5 M 


0| 
0 i 

X I 


Lm 1 • 5 M £R  £L  * M 


Xt  R>U+1.SM 
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StiCNkSS  Z 

o C ° 
o o ° 


C 0 DF.RF.S  RON 

SYNTAX  t CODERF.S  MATRIX.  Oh  CODESkF.S  ME  DSC  LIS  a 

SU  bPRCGRArn  Abb 
URDU  Pi  TNCNAY 

DESr hJPTICii  i CODER ES  PkC DUCES  A DIAGNOSTIC  ARRAY 
Tar  nESI DUALS  FROM  MEDIAN  POLISHING  CF  A TNC  ~N  A Y 
SHOULD  SET  RESIDS-i  In  hEDPOLlSH  WREN  CCDEhr:S 
OUTPUT  MATRIX  CONSISTS  CF  I tit.  SYMBOLS  ~°f  DF.PFNDI 
OF  Tak  RESIDUALS . 

VALUE  bt.  LON  LONER  QUAeTILE 
o t V ALU  E b t.  T N E F.  N QUART  I L E S 
ft  VALUE  LARGER  THAN  UPPER  QUART ILLS 


E A 

TA.bLh 


■nr  Ulr: 


CODER E S 


a LExt  St. 


NR  AT  RAP  PE.  NS  NnEn  t\ESIuS  = . 


Rr.  S m us*~  / 


t,*-Slt.U POLISH  ut  AT kS 

315.4 
4 0.3 
3 8.15 
37.575 

TYr'iCAL  VALUE  1 2 5.6 
h ON  t F r EOT S t 15.3625  10,1 


12.6375  32.875 


COLUMN  EFFECTS t 2.7  3.9  12.625  3.1375 
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SCAT  Z 

RANGE  OF  Xi  "lO  20 
RANGE  OF  Y i ~8  6 


0 I 


O 


O 


fl  rtfffffff  5ff£W5  TC  iff  /)  RELATIONSHIP  BETWEEN  THE  RESIDUALS 
n AND  THE  COMPARISON  VALUES,  ( ffff • . *CE ) iTV . 


A+-LINE  Z 

SLOPE i 0.6216761995  Y-INTERCEPT i '0.1348011389 


ANCVAHCW 
SYNTAXi  ANOVA 
PARAMETER i 

OPTION-  SELECTS  THE  TYPE  OF  ANOVA  PERFORMED.  OPTION* 1 
CI/ff5  /I  ONE-WAY  ANOVA  AND  OPTION =2  GiTffff  /l  2V<T-V/jy 
MCW/  {TREATMENTS  WITH  BLOCKING).  ( DEFAULT =1) 

GROUPt  TWOWAY 

SUBPROGRAMS i INPUT,  SUMSQ,  FMT 

DESCRIPTION \ ANOVA  DOES  A ONE-WAY  OR  A TWO-WAY  ANALYSIS  OF 
VARIANCE  DEPENDING  ON  THE  VALUE  OF  THE  OPTION  PARAMETER.  ANOVA 
WILL  INTERACT  WITH  THE  USER  TO  OBTAIN  THE  REQUIRED  INFORMATION 
AND  DATA.  WHEN  ENTERING  DATA,  SEPARATE  THE  DATA  POINTS  WITH  AT 
LEAST  ONE  BLANK  AND  USE  NO  OTHER  DELIMITERS . OUTPUT  CONSISTS  OF 
AN  ANOVA  TABLE,  ESTIMATES  OF  THE  OVERALL  MEAN,  AND  ESTIMATES  OF 
THE  TREATMENT  {AND, IF  OPTION =2,  BLOCK)  EFFECTS. 
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I 


OPTION 


An  OVA 

ENT  Eh  NUMBER  Of  TREATMENTS . 


ENTER  y ECTOR  OF  NUMBER  OF  CBS.  FOR  TREATMENTS  1 70  u 
J ! 


3 3 3 3 

E it  T F.  n 

3 OBSERVATIONS  FOR 

TREATMENT 

1 

| 5 

45  46  51 

F.  n TER 

3 OBSERVATIONS  FOE 

7 RE  AT  ME  NT 

2 

42  44  50 

ERIE  h 

3 0 BS  ER\A  AT  ICRS  FOR 

treatment 

3 

u « 

36  41  43 

E n 7 r.  R 

3 CbSERV ATIOnS  FOR 

treatment 

4 

49  47  54 


SOURCE 

Tt\f.  ATmF.nT  3 313.9  2 36.97 

ERROR  3 354.00  19.25 

I CT  AL  11  ^'o4.92 

A ~S  tU  A l\R  - 0.419 

Cf/r.RALL  REAR  - 4 6.08 

TREATMENT  EFFECTS  1.25  "0.75  _4.42  3.9^ 


A n 0 ’/ A TABLE 

c c 


36.97  3.32 

19.25 


a LET  l'S  TAKE  THE  SAME  DATA  AND  DC  A TWO  JAY  ANCVA 


OPTION*-  2 


ANCVA 

ENTER  NUMBER  OF  TREATMENTS. 

□ i 

4 

ENTER  IE£ZQ.&  OF  NUMBER  OF  OBS.  FOR  TREATMENTS  1 TO  4 

□ i 

3 3 3 3 

ENTER  3 OBSERV ATIONS  FOR  TREATMENT  1 

□ i 

45  46  51 

ENTER  3 OBSERV  ATIONS  FOR  TREATMENT  2 

Ul 

42  44  50 

ENTER  3 OBSERVATIONS  FOR  TREATMENT  3 

□ « 

36  41  48 

ENTER  3 OBSERVATIONS  FOR  TREATMENT  4 

□ « 

40  47  54 

ANCVA  TABLE 

SOURCE  DF  SS  MS  F 

TREATMENT  3 60.67  20.22  5.56 

BLOCKS  2 204.17  102.08  28.05 

ERROR  6 21.83  3.64 

TOTAL  11  286.67 

RmSQUARE  = 0.924 

OVERALL  ME  An  = 45.33 

TREATMENT  EFFECTS  2.00  0.00  "3.67  1.67 

BLOCK  EFFECTS  "4.58  "0.83  5.42 


CCtiTIliGEmCYhOW 
SYNTAX  i CONTINGENCY  MATRIX 

GhCUPi  TV  OH  AY  , 

DESCRIPTION  \ CONTINGENCY  TAKES  A T VC -WAY  TABLE  AliJJ  _ 

TEST  CF  INDEPENDENCE.  OF  ROWS  AND  COLUMNS . * ue-  •'  itl-c  n 

ST  AT  I ST  IP  A.il)  ITS  DEGREES  OF  FREEDOM  ARE  OUTPUl  . 


CONTINGENCY  DEATHS 

CHI-SQUARE  V ALU  E-  2.920832627  DF-  12 


fl 

n 

n 

a 


COMPARE  THE  ABOVE  VALUE  WITH  THOSE  IN  r ,il -SQUARE 
TABLES  FOR  12  Dt  TO  DETERMINE  IF  THE  HYPOTHESIS 
STATING  THAT  THE  ROWS  AND  COLUMNS  ARE  In  DEPENDENT 
CAu  BE  REJECTED  AT  THE  SELECTED  SIGNIFICANCE  LEVEL . 
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) GRP  STATS 

CORRELATION  CR  FMT 

CQRRELATICNHC N M ST  AT SHOW 


MET  ATS  PARTIAL  STATISTICS 
P ART  I ALLOW  STATISTICS  ti  C i 


ST  AT 1ST IC SHOW 

SYNTAX  i STATISTICS  l /ECTOR 
GROUP i STATS 
SUBPROGRAM i fAjr 

DESCRIPTION  i STATISTICS  DETERMINES  THE  ME  AD,  VARIANCE, 
DEVIATION,  COEFFICIENT  OF  VARIATION,  LO.VER  AND  UPPER 
:-i EDI  An,  TRIieAN,  AIDMF.An,  MIDRANGE,  RAnGE,  mEAN 

deviation,  interquartile  range,  s ken  ness,  and  kurtcsis 

DATA  IN  VECTOR . 


ST  An D A n d 
QUART ILE E , 
A DSC  LUTE 
E C n TEE 


STATISTICS  VOLCANO 
ME  AN t 7 0.246575  34 
VARlAuCEi  1850.562775 
STD.  DEV.i  43.01816796 
COEFF.  OF  VARIATION i 0.6123881165 
LONER  QUART ILE \ 37 
UPPE n QUAhTILEi  96 
ME  01  Alii  6 5 
T RIME  Am  6 5.75 
MI  DAE  AN i 64.47747748 
R A:iG E i 19  7 
Ml D R AN G t.  t 100.5 

MEAN  ABSOLUTE  DEVIATION:  33.68493151 

Interquartile  range t 59 

COEFF.  OF  SKEWNESS:  0.8325767059 
COEFF.  OF  KURTCSIS : 0.4281831164 


CCRRELATICNHON 
SYNTAX  i R+CCRRELATICti  W 
GROUP : STATS 

SUBPROGRAM:  FMT 

DESCRIPTION  I CORRELATION  DETERMINES  THE  SIMPLE  PE  ARSON  PRODUCT 
MOMENT  CORRELATIONS  BETWEEN  EACH  PAIR  OF  VARIABLES  RF.PRE SENTE  D bY 
THE  0 COLUMNS  OF  W.  THE  OUTPUT  IS  A C bY  C CORRELATION  MAThIX. 
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CORRELATION 

GkADb.F  1 

1.00 

2.00 

3.00 

4.00 

5.00 

0 0 1.00 

0.34 

0.5  0 

0.^5 

0.64 

0 0 0.34 

1 .00 

"0.07 

0.54 

0.68 

00  0.30 

"0.07 

1 .00 

0.52 

0.64 

00  0.25 

0.54 

0.52 

1.0  0 

0.8  7 

00  0.G4 

0.53 

0.6  4 

0.87 

I . J 0 

p 

1/7  THE  CORRELATION  OUTPUT , THE  TCP 

ROW 

A.'ib  Tub.  b it.  FT 

p 

COLUMN  ARE  LAHELF  FOh  Tat.  CCRRFLATI 

OIF. 

FOR  EX  A.-.  PL F 

p 

THE  CORRELATION  BETWEEN  THE  FFCOWu 

A'u  U 

FOURTH  COLU 

p 

OF  GhAub.Fl  IF  0.54. 

pGRAUEFl 


FYNTAXi  MF  I 
GROUP t FT  AT 
FUuFh  OGHAM : 
it  RFC  hi  FT  IF  in  i 


LEV  I AT  I'' 
CCLV  -mF 

lit  Ait  :t 


I CNF, 
OF 
BY 


b i'jT 

,\FT  AT  F UETERnlHEF  CCRRF.LATIO  aF  , 
LOWER  itUARTILEF,  ME ul AiiF , ALL/  UPPER  „ 
at.  An(iL"4EuT  tiATRIX.  W . TtiF.  CORRF.LATK 


]U  ART  I L 
1 u F A tit 


" j ' A n L'  /i . L 
o’  t C it  f i r 
L 1 F PL/i  l r L 


I MATRIX , WtiF.RE  :<i  IF  THE  HUM  HER  Cb  COLUkWF 


Tab.  OUTPUT  FhC 


iUCti  LIKE  T:i A j 


?f  COKRb  LA?  [C 


1 1 FT  AT  F (jRAuEFi 


1.00 

2.0  0 

3.00 

4.00 

5 . 0 0 

1.00 

1.00 

0.34 

0.50 

0.25 

0.6  4 

2.00 

0.34 

1.00 

"0.07 

0.54 

0.68 

3.0  0 

0.50 

"0.07 

1 .00 

0.52 

0.64 

4.00 

0.25 

0.54 

0.52 

1.00 

0.37 

5.00 

0.64 

0.68 

0.64 

0.87 

1 .00 

i’JEAuF 

0.86 

0.31 

0.88 

0.81 

0.8  4 

FTU.utV 

0.00 

0.11 

0.10 

0 . 1 0 

0.0  7 

L.  uni LF 

0.8^ 

0.69 

0.92 

0.74 

0.7  8 

. ;t.  ul  AnF 

0.85 

0.83 

0.88 

0.79 

0.83 

V.  uR'TLF 

0.0  4 

0.92 

0.93 

0.91 

0.9  0 
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PARTI ALHCW 
SYNTAX t PARTIAL  W 
GROUP’.  FIATS 

SUBPROGRAMS  : CCRRR  LAT ICN  , PUT 
DESCRIPTION  : PARTIAL  DETERMINES 
SPECIFIED  ORDER  FOR  THE  COLUMNS  CF  THE 
ASKED  FOR  TdE  COLUMNS  TO  BE  PART I ALLE D 

must  be  Integers  between  i and  c , wu 
COLUMNS  OF  IV. 


PARTIAL  CORRELAT 


MATRIX  W. 
no,  2 RE 
r is 


CUT 

:EE 


PARTIAL  GrADESI 


EXT  Eh  rCLUi-ui 

5 CF 

tc 

BE  PART 

I ALLE d 

J : 

A 

i 

1.00 

2.00 

3.00 

4.00 

M 

« 

o 

o 

1.00 

"0.29 

0.50 

0.64 

2.03 

0.29 

1.00 

0.4  7 

0.49 

3.0  0 

0.50 

0.47 

1.00 

0.95 

4.00 

0.64 

0.49 

0.9  5 

1.00 

1 UT 


T U r. 
USER 
TnE 


CNR 
USEE 
'S  hr. 

ju  P V 


i f H it  i 

TILL  u r. 
SPC'.Sr.S 
t.  r.  R 0 1 


* of  CArEFUL  WITH  INTERPRETING  THE  LABELS  IN  TdE  OUTPUT 
n OF  PARTIAL.  ABOVE , LABEL  1 DOES  NOT  MEAN  COLUMN  ONE  , 

« BUT  TdE  FIRST  COLUMN  REMAINING  ASTER  TilE  SELECTED  VARIABLES 
a HAVE  BEEN  PARTI ALF.D  CUT.  SIMILARLY,  2 REFERS  TO  Tnv. 
a SECOND  REMAINING  COLUMN. 


PARTIAL  GRADES! 

ENTER  COLUMNS  CF  MATRIX  TO  BE  P ART I ALLE  D CUT. 
Li  I 


2 3 


1.00 

2.00 

3.00 

1.00 

1.00 

"0.45 

0.23 

2.00 

"0.45 

1.00 

0.77 

3.00 

0.23 

0.77 

1.00 

a THE  ABOVE  EXAMPLE  SHOWS  THAT  mCRE  ThAN  CNF.  VARIABLE 
a CAN  BE  PART  I ALF.D  CUT.  LABEL  2 ABOVE  REFERS  TO  THE 
a SECOND  COLUMN  REMAINING  AFTER  THE  SELECTED  VARIABLES 
a HAVE  BEEN  PArTIALED  CUT . IN  TdlS  EXAMPLE  2 REFERS  TO 
a COLUMN  4. 
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)GRP  GFIT 

KS  KSHC'W  ChISQUARE  CHISQUAREdO W 


C'HISQU  ARF.tiOW 
SYNTAX:  T CHISQUARE  X 


GROUP  t GFIT 


DESCRIPTION  i CHISQUARE 
PROBABILITY  HASP  FUNCTION 
EMPIRICAL  FREQUENCY  Fl'NCTI 
WELL  THE  PROBABILITIES  7 
GCODNESSvCF- F Z'Z  S 7 A 7 ZF'ZZC 
DEGREES  OF  FREEDOM, 


COMPARES  A THEORETICAL  DISC  RET'.' 

T HAZZNG  i\l  PCSSZbLE  VALUES  WITH  A / 
01 1/  X OVER  THE  SAME  N VALUES  70  TEST  aOW 
FZ7  THE  DATA  X . THE  C til -SQUARE 

ZS  PRINTED  OUT  WITH  THE  DUMBER  OF 


a THE  THEORETZC AL  PROBABILITY  MASS  FUNCTION  2 AND  THE 
a EMPIRICAL  FREQUENCY  FUNCTION  X MUST  oE  OF  THE  SAmF 
a SIZE;  I.E. , p?=pX.  THE  ACTUAL  OBSERVED  VALUES  ARE 
a IMPLICIT  IN  THE  POSITION  OF  THE  VECTORS. 


T *-  .1  .2  .3  .2  .1 
X*-  5 20  38  32  1 2 


A 

A 

A 

A 

A 

A 


THE  VALUES  IN  T REPRESENT  PROBABILITIES . THE  VALUES  IN 
X REPRESENT  OBSERVED  FREQUENCIES . THE  ACTUAL  ObSFRVATZC 
STATED  EXPLICITLY , BUT  IT  IS  UNDERSTOOD  ThAT 


THE  POSITIONS  OF  X AND 
EXAMPLE,  IN  THE  CASE  AD 
TO  THE  VALUE  "l,  THE  SECO 


REFER  TO  THE  SAME 
VF,  TaE  FIRST  FOSI 
VD  0,  ETC. 


VALUES. 

'ION  HA Y 


1 L fi 

■FER 


7 CHISQUARE  X 

CHI-SQUARE  VALUE  = 9.62087227U  DF  = 4 


a TO  DETERMINE  IF  ONE  SHOULD  hF*hC?  THE  HYPOTHESIS  THAT 
a THE  DATA  WERE  DRAWN  FROM  A POPULATION  WITH  PROBABILITY 
a MASS  FUNCTION  7,  COM? ARE  THE  VALUE  ABOVE  WITH  THOSE 

a in  cm -square  tables  with  h de  . 
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KShOW 

SYNTAX i F XS 
GROUP t OF IT 
DESCRIPTION i 
OF  THE  EQUALIj 
DISTRIBUTIONS. 
0 bS  EnV AT  IONS . 
SIZE.  OUTPUT 


IS  DOES  A ROLHOGORGV  rSMIRNOV  GOODR-SS  OF  hi? 

■{  ct'  TWO  DISCRETE  EMPIRICAL  PRC  BAD 
TUF  ARGUMENTS  F ADD  G ARE.  F.A.Cti  I /ECO  Or, 
THEY  NEED  NOT  be.  SORTED , NOR  MUST  TREY  tsE  THE 
"F  KS  ARE  TdF.  VALUE  St 

MAX  F(X)fG(X)  AND  MIN  F(X)iG(X) 


TEST 
T LIT  Y 
S 0 E 

sale: 


Ft-  3 2 5 9 19  16  3 5 


71920126388977772137 


2102013  367  4^11386  3317 


E Ei S G 

.i^A'  HX'fGiX)  =3.1274  5 098  3 4 
FU)rG(X)  ="  0.31  372  5493  2 


A TO  DETERMINE  IP  THE  TWO  SAMPLES  Oh  uAT A WERE  TAKe.N  Ft,  ■ 
a THE  SAME  POPULATION  {OR  POPULATIONS  WITH  ThESA.-lF  CUF , 
n COMPARE  The.  LAX  AND  MIN  VALUE'S  WITH  THE  DEL 7 : V A\  lc  , - 
a IN  TAoLRS  FOR  TdF.  TWO  SAMPLE'  KS  TFST. 
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AN  DHOW 

SYNTAX’.  X AND  Y 

GROUPS:  RELATIONS,  TWOWAY , COMPARISONS 

DESCRIPTION:  THE  FUNCTION  'AND'  IS  USED  TO  CREATE  A NEW  DATA 

ARRAY  CONSISTING  OF  X AUGMENTED  BY  Y AS  ADDITIONAL  COLUMNS.  X 
AND  Y CAN  BE  SCALARS , VECTORS , OR  MATRICES  BUT  CANNOT  BOTH  BE 
SCALARS.  IF  AN  ARGUMENT  IS  A SCALAR,  A COLUMN  IS  GENERATED  EACH 
ELEMENT  OF  WHICH  IS  THE  SCALAR.  IF  THE  ARGUMENTS  ARE  VECTORS  OR 
MATRICES,  THEY  MUST  BE  CONFORMABLE. 


ANOVAHOW 
SYNTAX:  ANOVA 

PARAMETER: 

OPTION - SELECTS  THE  TYPE  OF  ANOVA  PERFORMED.  OPTION-1 
GIVES  A ONE-WAY  ANOVA  AND  OPTION-7  GIVES  A TWO-WAY 
ANOVA  ( TREATMENTS  WITH  BLOCKING).  {DEFAULTED 
GROUP:  TWOWAY 

SUBPROGRAMS:  INPUT,  SUMSQ,  FMT 

DESCRIPTION : ANOVA  DOES  A ONE-WAY  OR  A TWO-WAY  ANALYSIS  OF 

VARIANCE  DEPENDING  ON  THE  VALUE  OF  THE  OPTION  PARAMETER.  ANOVA 
WILL  INTERACT  WITH  THE  USER  TO  OBTAIN  THE  REQUIRED  INFORMATION 
AND  DATA.  WHEN  ENTERING  DATA,  SEPARATE  THE  DATA  POINTS  WITH  AT 
LEAST  ONE  BLANK  AND  USE  NO  OTHER  DELIMITERS . OUTPUT  CONSISTS  OF 
AN  ANOVA  TABLE , ESTIMATES  OF  THE  OVERALL  MEAN,  AND  ESTIMATES  OF 
THE  TREATMENT  (.AND, IF  OPTION-7,  BLOCK)  EFFECTS. 


A3RH0W 

SYNTAX:  A3R  VECTOR 

GROUP:  SMOOTH 

SUBPROGRAM:  ONE  3 

DESCRIPTION : A3R  DOES  REPEATED  SMOOTHINGS  OF  RUNNING  MEDIANS  OF  3 
UNTIL  THERE  ARE  NO  CHANGES  IN  THE  SMOOTHED  SEQUENCE  FROM  ONE 
ITERATION  TO  THE  NEXT.  SEE  0NE3H0W. 


A3RSRH0W 

SYNTAX:  A3RSR  VECTOR 

GROUP:  SMOOTH 

SUBPROGRAMS:  A3R,  ONE  3,  SPLIT 

DESCRIPTION : A3RSR  DOES  REPEATED  SMOOTHINGS  BY  RUNNING  MEDIANS 

OF  3 FOLLOWED  BY  SPLITTING  MESAS  AND  REPEATING  UNTIL  CONVERGENCE . 
SEE  A3R  AND  SPLIT. 
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BOXPLOT  HOW 

SYNTAX : BOXPLOT  VECTOR 

PARAMETER : 

LENGTH-  CONTROLS  THE  HORIZONTAL  SIZE  OF  THE  DISPLAY 
( DEFAULT^ 50  CHARACTERS) . 

GROUP : DISPLAY 

SUBPROGRAM : FILL 

DESCRIPTION : BOXPLOT  GENERATES  A BOX  PLOT  DISPLAY  FOR  A VECTOR 

OF  DATA.  A RECTANGULAR  BOX  WITH  ENDS  CORRESPONDING  TO  LONER  AND 
UPPER  QUARTILES  IS  PRESENTED  WITH  THE  MEDIAN  MARRED  WITH  AN 
ASTERISK.  WHISKERS  ARE  DRAWN  ON  EACH  SIDE  OF  THE  BOX  WITH 
CROSSES  MARKING  THE  LOWEST  AND  HIGHEST  DATA  VALUES  WITHIN  AN 
INTERQUARTILE  DISTANCE  OF  THE  QUART! LES.  DATA  VALUES  OUTSIDE  THE 
CROSSES  ( OUTLIERS ) ARE  MARKED  WITH  CIRCLES  AND  THOSE  MURE  THAN 
1.5  INTERQUARTILE  DISTANCES  GET  HEAVY  CIRCLES. 


CHISQU AREHOW 
SYNTAX : T CHISQU ARE  X 

GROUP : GFIT 

DESCRIPTION:  CHISQU ARE  COMPARES  A THEORETICAL  DISCRETE 

PROBABI LITY  MASS  FUNCTION  T HAVING  N POSSIBLE  VALUES  WITH  A<i 
EMPIRICAL  FREQUENCY  FUNCTION  X OVER  THE  SAME  N VALUES  TO  TEST  HOW 
WELL  THE  PROBABI LITIES  T FIT  THE  DATA  X.  THE  CHI-SQUARE 
GOODNESS-OF-FIT  STATISTIC  IS  PRINTED  OUT  WITH  THE  NUMBER  OF 
DEGREES  OF  FREEDOM. 


CODERESHOW 

SYNTAX:  CODERES  MATRIX  OK  CODESRES  MEDPOLISH  MATRIX 

SUBPROGRAM:  AND 

GROUP:  TWOWAY 

DESCRIPTION:  CODERES  PRODUCES  A DIAGNOSTIC  ARRAY  FUR  ANALYZING 

THE  RESIDUALS  FROM  L.UI AN  POLISHING  OF  A TWO-WAY  TABLE.  THE  USER 
SHOULD  SET  RESIUS=  1 IN  MEDPOLISH  WHEN  CODERES  IS  USED.  Tiic. 
OUTPUT  MATRIX  CONSISTS  OF  THE  SYMBOLS  -°+  DEPENDING  ON  THE  SIZE 
OF  THE  RESIDUALS. 

- : VALUE  BELOW  LOWER  QUARTILE 

° : VALUE  BETWEEN  QUARTILES 

+:  VALUE  LARGER  THAN  UPPER  QUARTILES 
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INDICATE  HISSING  /ALU ES 


C OMP  ARE  HO*/ 

SYNTAX:  COMPARE  MATRIX 

PARAMETERS: 

(1)  DEPTH - VERTICAL  HEIGHT  OF  DISPLAY  {DEFAULT-20 
LINES) . 

(2)  MISS-  NUMBER  USED  TO  INDICATE  MISSING  VALUES 
( DEFAULT  =~99999) 

(3)  NGAP-  NUMBER  OF  HORIZONTAL  SPACES  BETWEEN  THE 
BOXPLOT  DISPLAYS  ( DEF AU LT = 3) . 

GROUP:  COMPARISONS 

SUBPROGRAM:  FILL 

DESCRIPTION : COMPARE  OPERATES  ON  AN  N BY  K MATRIX  TO  PRODUCE  X 

VERTICAL  BOX  PLOTS  PLACED  NEXT  TO  EACH  OTHER  TO  ALLOW  VISUAL 
COMPARISON  OF  THE  CENTERS , SPREADS , AND  OUTLIERS  OF  THE  BATCHES 
( COLUMNS ) OF  THE  MATRIX.  THE  USER  MUST  FILL  UP  THE  MATRIX  SO 
THAT  THERE  ARE  N OBSERV ATIONS  FOR  EACH  BATCH.  THE  PARAMETER, 

MISS,  SHOULD  BE  USED  TO  FILL  IN  MISSING  VALUES.  THE  USER  SHOULD 
ASSURE  THAT  MISS  IS  DIFFERENT  FROM  ALL  VALID  ENTRIES  IN  THE  DATA 
ARRAY . 


CONDENSEHOW 

SYNTAX:  CONDENSE  MATRIX  OR  R+C ON  DENSE  MATRIX 

PARAMETERS : 

(1)  NUM-  CONTROLS  WHAT  STATISTICS  ARE  INCLUDED  IN  THE 
SUMMARY.  NU  Me [1,2, 5, 7]  ( DEFAULT  = 2 ) . 

NUM= 1 GIVES  MEDIANS  OF  EACH  COLUMN  OF  THE  MATRIX 
NUM-2  GIVES  MEDIANS  AND  INTERQUARTILE  RANGES 
NUM - 5 GIVES  MIN , QU ART  I LES , MAX , AND  SIZE 
NUM-1  GIVES  MIN , EIGHTHS , MAX  AND  SIZE 

(2)  MISS-  NUMBER  USED  TO  CODE  MISSING  VALUES 
( DEF  AU  LT- ~ 9 9 9 9 9 ) . 

GROUP:  COMPARISONS 

SUBPROGRAMS:  UTCOND,  FMT 

DESCRIPTION : CONDENSE  GENERATES  SUMMARY  STATISTICS  FOR  EACH 

GROUP  OF  DATA  REPRESENTED  BY  THE  COLUMNS  OF  THE  ARGUMENT  MATRIX. 
SINCE  THE  ARGUMENT  IS  A MATRIX,  THE  USER  MUST  FILL  UP  THE  MISSING 
VALUES  IN  THE  MATRIX  WHENEVER  THE  GROUPS  ( COLUMNS ) HAVE  DIFFERENT 
NUMBERS  OF  OBSERVATIONS . THE  PARAMETER,  MISS,  SHOULD  BE  USED  TO 
FILL  THE  MATRIX  SINCE  THE  FUNCTION  WILL  RECOGNIZE  THOSE  VALUES  AS 
MISSING  DATA  AND  WILL  IGNORE  THEM  IN  ALL  CALCULATIONS.  THE 

RESULT  FROM  CONDENSE  CAN  BE  USED  AS  THE  ARGUMENT  OF  OTHER 
FUNCTIONS  SUCH  AS  SCAT  OR  LINE  IF  AND  ONLY  IF  NUM= 2. 
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CONTI NGENCY HOW 
SYNTAX:  CONTINGENCY  MATRIX 

GROUP:  TWOWAY 

DESCRIPTION : CONTINGENCY  TAKES  A TWO-WAY  TABLE  AND  PERFORMS  A 

TEST  OP  INDEPENDENCE  OF  ROWS  AND  COLUMNS.  THE  CHI-SQUARE 
STATISTIC  AND  ITS  DEGREES  OF  FREEDOM  ARE  OUTPUT. 


KSHOW 

SYNTAX:  F KS  G 

GROUP:  GFIT 

DESCRIPTION : KS  DOES  A KOLMOGOROV -SMI RNOV  GOODNESS  OF  FIT  TEST 

OF  THE  EQUALITY  OF  TWO  DISCRETE  EMPIRICAL  PROBABILITY 
DISTRIBUTIONS . THE  ARGUMENTS  F AND  G ARE  EACH  VECTORS  OF 
OBSERV  AT  IONS . THEY  NEED  NOT  BE  SORTED , NOR  MUST  THEY  BE  THE  SAME 
SIZE.  OUTPUT  OF  KS  ARE  THE  VALUES: 

MAX  F(X)-GU)  AND  MIN  F(.X)-G(X) 


LIHEHOW 

SYNTAX:  LINE  ARRAY  OR  INLINE  ARRAY 

PARAMETER: 

RSELECT - SELECTS  THE  OUTPUT  MATRIX  OF  RESIDUALS. 

RS ELECT = 1 GIVES  ABSCISSA  VS.  RESIDUALS  ( DEFAULT- 1 ) . 
RSELECT*!  GIVES  FITTED  VALUES  VS.  RESIDUALS 
GROUP:  RELATIONS 

DESCRIPTION : LINE  FITS  A STRAIGHT  LINE  TO  A SET  OF  (X,Y)  POINTS 

BY  DIVIDING  THE  POINTS  INTO  3 REGIONS  AND  USING  THE  MEDIANS  OF 
THE  X AND  Y VALUES  IN  THE  OUTER  REGIONS  TO  DETERMINE  THE  SLOPE. 
THE  INTERCEPT  IS  THE  MEDIAN  OF  THE  DIFFERENCES  Y-SLOPE*X.  LI'ir. 
ALSO  COMPUTES  THE  hESIDUALS  AND  GIVES  AN  N BY  2 MATRIX  CONTAINING 
EITHER  THE  X -VALUES  VS.  RESIDUALS  OR  THE  FITTED  VALUES  VS. 
RESIDUALS  DEPENDING  ON  THE  PARAMETER,  RSELECT.  IF  THE  USER  DOBS 
NOT  WANT  RESIDUALS  TYPED  OUT,  HE  SHOULD  USE  THE  SYNTAX: 

Z-^LINE  ARRAY 

THE  ARGUMENT  ARRAY  CAN  BE  EITHER  A VECTOR  OR  A 2-COLUMN  MATRIX. 

IF  THE  ARGUMENT  TS  A VECTOR  OF  SIZE  N,  THE  X -V  ART  ABLE  IS 
CONSTRUCTED  TO  BE  THE  FIRST  N POSITIVE  INTEGERS.  IF  THE  ARGUMENT 
IS  A 2-COLUMN  MATRIX,  THE  FIRST  COLUMN  TS  TAKEN  TO  BE  THE  SET  OF 
X-VALUES  AND  THE  SECOND  COLUMN  THE  SET  OF  Y-VALUES.  LINE  CAN  BE 
USED  IN  CONJUNCTION  WITH  THE  FUNCTION  SCAT  TO  PRODUCE  THE  SLOPE , 
THE  INTERCEPT,  AND  A PLOT  OF  THE  RESIDUALS . TO  DO  THIS,  ENTER: 

SCAT  LINE  ARRAY 
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LSLI NEHOW 

SYNTAX : LSLI HE  ARRAY  OR  Z-LSLINE  ARRAY 

GROUP : RELATIONS 

DESCRIPTION : LSLINE  DETERMINES  THE  LEAST  SQUARES  SOLUTION , 

Y=YI+SLOPE*X . THE  ARGUMENT  ARRAY  CAN  BE  A VECTOR  OR  A 2-COLUMN 
MATRIX.  IF  THE  ARGUMENT  IS  A VECTOR,  THE  X -VALUES  ARE  TAKEN  TO 
BE  THE  INTEGERS  1 TO  N INHERE  N=p ARRAY.  IF  THE  ARGUMENT  IS  A 
2-COLUMN  MATRIX,  COLUMN  1 CONSISTS  OF  THE  X-VALUES  AND  COLUMN  2 
THE  Y -VALUES.  RESIDUALS  ARE  DETERMINED  AND  STORED  IN  A MATRIX 
WHOSE  FIRST  COLUMN  IS  THE  SET  OF  X-VALUES  AND  WHOSE  SECOND  COLUMN 
IS  THE  SET  OF  RESIDUALS.  IF  THE  USER  WANTS  TO  SUPPRESS  PRINTOUT 
OF  THE  RESIDUALS,  HE  MUST  ASSIGN  THE  RESULTS  OF  LSLINE  TO  A 
VARIABLE-,  I.E.,  TYPE  Z+LSLINE  ARRAY. 


MEDPOLISH HOW 

SYNTAX : MEDPOLISH  MATRIX  OR  Z-MEDPOLISH  MATRIX 

PARAMETERS: 

(1)  RESIDE-  CONTROLS  OUTPUT  OF  RESIDUALS.  IF  1, 

RESIDUALS  ARE  PRODUCED  IN  A TWO-WAY  TABLE:  IF  2, 

COMPARISON  VALUES,  ( RE ° . *C E) iTV  AND  RESIDUALS  ARE 
PRODUCED  IN  A TWO  COLUMN  MATRIX  {DEFAULTED  . 

(2)  EPSILON-  PROPORTION  BY  WHICH  THE  SUM  OF  THE 
ABSOLUTS  VALUES  OF  THE  RESIDUALS  MUST  BE  REDUCED  AT 
EACH  ITERATION  TO  CONTINUE  POLISHING  ( DEFAULT-0 .01). 

(3)  NORMEFFECTS - DETERMINES  WHETHER  OR  NOT  NORMALIZED 
EFFECTS  ARE  OUTPUT.  THE  DEFAULT  OF  0 SUPPRESSES  OUTPUT 
OF  NORMALIZED  VALUES.  N 0 RM E F F ECT S = 1 CAUSES  NORMALIZED 
VALUES  TO  BE  PRINTED. 

GROUP : TWOWAY 

DESCRIPTION : MEDPOLISH  ITERATIVELY  SWEEPS  OUT  MEDIANS  FROM  THE 

ROWS  AND  COLUMNS  OF  A TWO-WAY  TABLE  TO  YIELD  THE  MODEL: 

OBS  = MEDIAN  + ROW  EFFECT  + COLUMN  EFFECT  + RESIDUAL 
IT  YIELDS  RESIDUALS  FOR  TESTING  THE  ADEQUACY  OF  THE  MODEL  ( SEE 
SHOWRES  AND  CODERES).  THE  FUNCTION  CONTINUES  POLISHING  UNTIL  THE 
STOPPING  RULE  CONTROLLED  BY  EPSILON  IS  ACTIVATED.  THE  OVERALL 
MEDIAN,  THE  ROW-  AND  CO LU MN -E F F ECT S , AND  THE  SUM  OF  ABSOLUTE 
VALUES  OF  THE  RESIDUALS  ARE  GIVEN.  THE  USER  CAN  SUPPRESS 
RESIDUAL  PRINTOUT  BY  ASSIGNING  THE  FUNCTION  TO  A VARIABLE,  I.E., 
BY  WRITING: 

Z*-  MEDPOLISH  MATRIX 


CCRRELATICNHOW 
SYNTAX  « R+-CCRRELAT  ION  W 
GROUP x STATS 
SUBPROGRAM i FMT 

DESCRIPTION t CORRELATION  DETERMINES  THE  SIMPLE  PEARSON  PRODUCT 
MOMENT  CORRELATIONS  BETWEEN  EACH  PAIR  OF  VARIABLES  REPRESENTED  BY 
THE  C COLUMNS  OF  W.  THE  OUTPUT  IS  A C BY  C CORRELATION  MATRIX. 


UST  AT  SHOW 
SYNTAX  i M ST  ATS  W 
GROUP i STATS 
SUBPROGRAM t FMT 

DESCRIPTION x MSTATS  DETERMINES  CORRELATIONS , MEANS,  STANDARD 
DEVIATIONS , LOWER  QUART I LEE,  MEDIANS , AND  UPPER  QUART ILES  FOR  THE 
COLUMNS  OF  THE  ARGUMENT  MATRIX  W.  THE  CORRELATIONS  ARE  DISPLAYED 
IN  AN  N BY  N MATRIX,  WHERE  N IS  THE  NUMBER  OF  COLUMNS  OF  W, 


PARTIALHCW 
SYNTAX x PARTIAL  W 
GROUP x STATS 

SUBPROGRAMS t CORRELATION , FMT 

DESCRIPTION i PARTIAL  DETERMINES  PARTIAL  CORRELATIONS  OF  ANY 
SPECIFIED  ORDER  FOR  THE  COLUMNS  OF  THE  MATRIX  W.  THE  USER  WILL  BE 
ASKED  FOR  THE  COLUMNS  TO  BE  PARTIALLED  OUT.  THE  USER'S  RESPONSES 
MUST  BE  INTEGERS  BETWEEN  1 AND  C,  WHERE  C IS  THE  NUMBER  OF 
COLUMNS  OF  W. 
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NUMSUMHOW 

SYNTAX : NUMSUM  VECTOR  OR  NUMSUM  ARRAY 

PARAMETER: 

MISS - NUMBER  USED  TO  INDICATE  MISSING  VALUE  IN  THE  DATA 
ARRAY  (DEFAULT=~99999) . 

GROUP : COMPARISONS 

SUBPROGRAM : FMT 

DESCRIPTION : NUMSUM  OPERATES  ON  EITHER  A VECTOR  OR  AN  ARRAY  BUT , 

IN  ALL  CASES , THE  DATA  ARE  TREATED  AS  A SINGLE  BATCH.  NUMSUM 
PRODUCES  A NUMERICAL  SUMMARY  WHICH  GIVES  THE  SAMPLE  SIZE  ( AFTER 
DELETION  OF  MISSING  VALUES),  THE  EIGHTHS,  THE  EXTREMES,  THE 
SPREADS,  AND  THE  MIDPOINTS  IN  TABULAR  FORM.  THIS  SUMMARY  IS 
USEFUL  FOR  TESTING  THE  SYMMETRY  OF  A DATA  SET,  AND  TO  EVALUATE 
THE  EFFECTIVENESS  OF  DATA  TRANSFORMATIONS  IN  PRODUCING  SYMMETRY. 


ON EH  HOW 

SYNTAX:  ON  EH  VECTOR 

GROUP:  SMOOTH 

DESCRIPTION:  ON  EH  DOES  A SINGLE  HANNING  OF  A SEQUENCE  OF  DATA. 

IT  IS  USED  IN  CONJUNCTION  WITH  THE  OTHER  SMOOTHING  FUNCTIONS  AS  A 
FINAL  TOUCH-UP  TO  A SEQUENCE  OF  SMOOTHED  DATA.  THE  I -T  H RESPONSE 
IS  REPLACED  BY  0 . 2 5*Y[I -1 ] + 0.5xy[j]  + 0 . 2 5 * Y [ 1+ 1 ] . SEE  A2R  A'vO 
A 3 RSR . 


ONE  3 HOW 

SYNTAX:  ONE 3 VECTOR 

GROUP:  SMOOTH 

DESCRIPTION : ONE 3 SMOOTHS  A VECTOR  OF  DATA  USING  ONE  PASS  OF 

RUNNING  MEDIANS  OF  3 WITH  TU KEY'S  END-POINT  RULE.  WHEN  SMOOTHING 
(.X , Y ) PAIRS,  THE  USER  SHOULD  USE  AS  THE  ARGUMENT  VECTOR  THE  Y 
VALUES  ORDERED  ACCORDING  TO  THE  MAGNITUDE  OF  THE  X VALUES: 

VECTORS CiY] 

IT  IS  ASSUMED  THAT  THE  X VALUES  ARE  EQUISPACED.  REPEATED 

SMOOTHING  BY  MEDIANS  OF  3 CAN  BE  ACCOMPLISHED  WITH  A3R.  OUTPUT 
CONSISTS  OF  THE  SMOOTHED  SEQUENCE.  A PLOT  OF  THE  SMOOTHED 
SEQUENCE  IS  GIVEN  BY  TYPING: 

SCAT  0NE3  VECTOR 
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REGRESSHCW 

SYNTAX  I Z-Y  REGRESS  X 
PARAMETER x 

A INTERCEPT - DETERMINES  WHETHER  CR  NOT  AN  INTERCEPT  TERM 
IS  TO  HE  INCLUDED.  ^INTERCEPT- 1 INTERCEPT 

TERM,  AND  EINTERCEPT-0  GIVES  NO  INTERCEPT.  ( DEFAULT  IS 

1.) 

GROUP x RELATIONS 
SUBPROGRAMS  i FM?  5C4J7 

DESCRIPTION  t REGRESS  DOES  A MULTIPLE  REGRESSION  ANALYSIS 
RELATING  THE  DEPENDENT  VARIABLE  Y TO  A SET  OF  CARRIERS  X.  THE 
LEFT  ARGUMENT  Y IS  A VECTOR  OF  SIZE  N.  THE  RIGHT  ARGUMENT  X IS 
AN  N BY  K MATRIX  CONSISTING  OF  N OBSERVATIONS  ON  EACH  OF  K 
VARIABLES  OR  A VECTOR  OF  SIZE  N IF  K= 1.  OUTPUT  CONSISTS  OF  Ay 
A NOV  A TABLE,  RmSQU  ARE , STD.  ERROR,  REGRESSION  COEFFICIENTS  ( THE 
FIRST  COEFFICIENT  IS  THE  CONSTANT  TERM  IF  MNTERCEPT-1 . ) , 
TnSTATISTICS,  VARIANCE mCCVARI AN CE  MATRIX,  DURBIN "WATSON 

STATISTIC.  AND  A VECTOR  OF  PREDICTED  Y VALUES  AND  RESIDUALS. 
THERE  IS  AN  OPTION  THAT  ALLOWS  THE  USER  TO  INPUT  A VECTOR  OF  X 
VALUES  AND  USE  THE  REGRESSION  EQUATION  TO  FORECAST  Y VALUES.  THE. 
USER  CAN  ALSO  OBTAIN  A SCATTER  PLOT  OF  THE  RESIDUALS.  WHEN 
EXECUTION  TERMINATES , THE  PREDICTED  Y VALUES  AND  THE  RESIDUALS 
RESIDE  IN  THE  N BY  2 MATRIX  Z. 


STATISTICSHCW 

SYNTAX x STATISTICS  VECTOR 
GROUP x STATS 
SUBPROGRAM x FMT 

DESCRIPTION i STATISTICS  DETERMINES  THE  MEAN,  VARIANCE,  STANDARD 
DEVIATION,  COEFFICIENT  OF  VARIATION,  LOWER  AND  UPPER  QUARTILFS , 
MEDIAN,  TRIMEAN,  MIDMEAN,  MIDRANGE,  RANGE,  MEAN  ABSOLUTE 
DEVIATION,  INTERQUARTILE  RANGE,  SKEWNESS,  AND  KURTOSIS  FOR  THE. 
DATA  IN  VECTOR. 


ARRAY 


SCATHOH 
SYNTAX i SCAT 
PARAMETERS i 


(1)  HID ~ CONTROLS  THE  HORIZONTAL  SIZE 
( DEFAULT= 3 0 CHARACTERS ) . 

OF 

THE 

DISPLAY 

(2)  DEPf  CONTROLS  THE  VERTICAL  SIZE 
( DEF AULT=1 5 LINES). 

OF 

THE 

DISPLAY 

(3)  NDIVX,  NDIVY r NUMBER  OF  UNITS  ON 
RESPECTIVELY  (.DEFAULT^  4,4)  . 

X - 

AND 

Y-AXES, 

GROUP t RELATIONS , SMOOTH , COMPARISONS 

DESCRIPTION  i 5f/ir  PRODUCES  A SCATTER  PLOT  OE  THE  DATA  CONTAINED 
IN  ARRAY.  THE  ARGUMENT  ARRAY  CAN  BE  A VECTOR  OR  A MATRIX  H IT H AS 
MANY  AS  9 COLUMNS.  IF  A VECTOR  OF  SIZE  N,  THOSE  VALUES  ARE 
PLOTTED  VS.  THE  INTEGERS  1 TO  N-,  IF  A MATRIX,  THE  SECOND,  THIRL, 
ETC.  COLUMNS  ARE  PLOTTED  VS.  COLUMN  1 ON  THE  SAME  AXES. 
DIVISIONS  ON  THE  AXES  OF  THE  PLOT  ARE  HOT  EXPLICITLY  PRINTED, 
EXCEPT  AT  THE  EXTREMES  OF  THE  PLOT.  THE  USER  CAN  CONTROL  THE 
RESOLUTION  OF  THE  PLOT  BY  MODIFYING  THE  PARAMETERS  HID,  DEP, 
NDIVX  AND  NDIVY.  PRINTING  TIME  INCREASES  DRAMATICALLY  WITH 
RESOLUTION ; THEREFORE,  OTHER  PLOT  PROGRAMS  SHOULD  BE  USED  IF  HIGH 
RESOLUTION  IS  DESIRED.  FOR  A SINGLE  GROUP  OF  DATA,  THE  NUMBER  N 
(2S/V<9)  DILL  BE  PRINTED  IF  N POINTS  LIE  CLOSE  TOGETHER  ON  THE 
DISPLAY.  FOR  MULTIPLE  PLOTS  ON  THE  SAME  DISPLAY,  THE  LETTER  A 
REPRESENTS  GROUP  1,  B REPRESENTS  GROUP  2,  ETC.  FOR  2 POINTS 
CLOSE  TOGETHER  THE  LETTER  HILL  BE  PRINTED  H IT H AN  UNDERSCORE. 
THE  DISPLAY  CANNOT  HANDLE  3 OR  MORE  POINTS  CLOSE  TOGETHER  IF 
THERE  ARE  MULTIPLE  PLOTS. 


SHOWRESHOW 

SYNTAX:  SHOW  RES  MATRIX  OR  SHOWRES  MEDPOLISH  MATRIX 

SUBPROGRAM:  AND 

GROUP:  TWOWAY 

DESCRIPTION : SHOWRES  AIDS  IN  THE  ANALYSIS  OF  RESIDUALS  FROM 

MEDIAN  POLISHING  A TWO-WAY  TABLE.  THE  USER  SHOULD  SET  RESIDS= 1 
IN  MEDPOLISH  IF  SHOWRES  IS  TO  BE  USED.  THE  OUTPUT  OF  SHOWRES  IS 
A MATRIX  CONSISTING  OF  THE  SYMBOLS  0O°xX  DEPENDING  ON  THE  SIZES 
OF  THE  RESIDUALS.  LET  R REPRESENT  THE  RESIDUAL,  M THE  MIDSPREAD 
OF  ALL  RESIDUALS,  L THE  LOWER  QUARTILE,  AND  U THE  UPPER  QUARTILE. 
THE  CODES  ARE  PRINTED  AS  FOLLOWS: 

0:  R<L- 1.5/7 
O:  L-l.SM<.R<L-M 

U+M<R<U+1.3M 
X:  R>U+1.SM 


SPLITHOW 

SYNTAX:  SPLIT  A3R  VECTOR 

GROUP:  SMOOTH 

DESCRIPTION:  SPLIT  DOES  ONE  PASS  AT  DIVIDING  MESAS  {PAIRS  OF 

ADJACENT  POINTS  WITH  A COMMON  VALUE  WHICH  IS  A LOCAL  MAX  OR  MIN) 
USING  T'J KEY'S  END-POINT  RULE.  IT  IS  USED  IN  CONJUNCTION  WITH 
A3R. 
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STEMLEAFHOW 

SYNTAX:  ST EM LEAF  VECTOR 

GROUP:  DISPLAY 

PARAMETERS : 

(1)  WIDTH’  CONTROLS  THE  WIDTH  (.CHARACTERS  PER  LIVE)  OF 
THE  DISPLAY  (DEFAULT* TO). 

(2)  SCALE-  VARIES  THE  DEPTH  (STEM  INTERVAL)  OF  THE 
DISPLAY  IN  UNITS  OF  1,  2,  OR  0.5  TIMES  A POWER  OF  10. 
SELECT  AN  INTEGER  FROM  1 TO  3 (DEFAULT*!) . 

SUBPROGRAM:  LIT 

DESCRIPTION : STEMLEAF  GENERATES  A STEM  AND  LEAF  DISPLAY  OF  A 

VECTOR  OF  OBSERVATIONS . THE  FUNCTION  AUTOMATICALLY  SCALES  THE 
DATA  USING  A SCALING  ROUTINE  BASED  ON  THE  RANGE  AND  SIZE  OF  THE 
DATA  VECTOR.  BY  VARYING  THE  SCALE  PARAMETER  THE  USER  CAN  CHANGE 
THE  SIZE  OF  THE  STEM , BUT  IT  WILL  NOT  NECESSARILY  BE  SCALED  BY 
THE  AMOUNT  SPECIFIED.  TRY  VARIOUS  VALUES  LIKE  1,2,  AND  3 TO  SEE 
WHICH  CHOICE  YIELDS  THE  BEST  DISPLAY.  IF  A LEAF  CONTAINS  TOO 
MANY  CHARACTERS  TO  BE  PRINTED  ON  A LINE , THE  LINE  WILL  BE 
TRUNCATED  AND  THE  NUMBER  OF  TRUNCATED  CHARACTERS  WILL  BE  WRITTEN 
AT  THE  END  OF  THE  LINE.  IF  THE  STEM  INTERVAL  IS  TWO  TIMES  A 

POWER  OF  TEN , THE  LEAVES  MAY  CONTAIN  THE  CHARACTERS  A,B,C -J 

IN  ADDITION  TO  THE  DIGITS  0,1,...,  9.  IN  THIS  CASE , THE 
CHARACTERS  A TO  J REPRESENT  THE  NUMBERS  10  TO  19. 


TNICEHOW 

SYNTAX  : TWICE  VECTOR 
GROUP:  SMOOTH 

SUBPROGRAMS:  ONE  3,  SPLIT,  A3RSR 

DESCRIPTION:  TWICE  SMOOTHS  A SEQUENCE  OF  DATA,  THEN  SMOOTHS  THE 

RESIDUALS  AND  ADDS  THE  SMOOTHED  RESIDUALS  BACK  TO  THE  SMOOTHS 
DATA  TO  OBTAIN  THE  FINAL  SMOOTH.  THE  SMOOTHING  IS  DONE  USIN 
A3RSR. 
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nut: FAULT  VALUES  OF  TdE  PARAMETERS 


DSP 

1 5 

DEPTH 

20 

EPSILO'V 

0.01 

LEdIGTH 

50 

MISS 

*99999 

tvUTV  X 

4 

DDIVY 

4 

;1GAP 

3 

■WRRiE  FF  LCT 

0 

MUM 

2 

OPT T 0U 

1 

REST US 

1 

SCALE 

1 

RS ELECT 

1 

■VI 0 

30 

VI  uT  d 

70 

LI  liT  EhCEPT 

1 

QA 

91 


>6  f\J 


VAHD CU]V 
7 L+A  MD  BiCiD 

[1]  •>(  ( (2<pp4)v3<ppg)  , 0*  p pfi  )/  17  3 

[2]  B*~%B 

[3]  -»■(  ( (3  = ppti)Al*lpp£)  ,2  = ppi4)/  l7  7 

[4]  A + ,A 

[5]  -*•(  a/(  (p/4)  *1,0)  ,l*0~lp~24>p&)/16 

[ 6 ] /!<-(  ( ( D*p  A ) L D[  pA)  , 1)  pA 

C 7 ] -*-(l*pp£>)/9 

[8]  B*-{  ( ( pb  ) T ( 1 =pB  ) *lppA  ) , 1 )pB 

C 9 ] -*•(  (a/£*1,  lpp/$)  ♦ l*£>-lp“24>Ofi)/  16  11 

[10]  5-*-(((3  = pp5)pl),(lpp/4),lp<f>p5)p5 

[11]  -*-<3=pp6)/14 

[12]  (C+lptypA  ) p 0 ) , ( lpfipB ) p 1 ) \B 

[13]  -*0xppL[  ; if]^/4 

[14]  Z>(  1 , ( (C*-lp<t>pA  ) p 0 ) , (_l  + lp<t»p^)pl)\B 

[15]  ->0*p  p£[  ; ; 1 + if  ]<->3 

[16]  -*0  = pi>'  ARGUMENTS  OF  AMD  ARE  MOT  CONFORM  ABLE . ' 
L17]  'AN  ARGUMENT  OF  AMD  IS  OF  IMPROPER  RANK..' 

7 


V AMUi/  4[G]7 

V AtWVAiUKiBTSSi  MB  R ; R',MAX ; BBSS ; MSBL ; MSB ; V 755  ;£DF  ;TDF  iTSS -,Y  bAI 
TiBLOCK 

[1]  INPUT 

[2]  SWASU 

V 


V^3/t[G]V 

7 Z-*-/l 3 A X;yi;Y2;Xl;^2;iV;C 

[1]  Z-5^53  X 

[2]  +1*  ( pppX-*-Z  ) » 0*  + / | Z’X 

V 


7j43/?5S[j]V 

V Z*/13A'5tf  X 

[1]  a:--/; 3*  X 

[2]  Z--A2R  SPLIT  X 

[3]  + 2*  ( pppA'-*-Z  ) « 0*  + / | Z-X 

V 
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vboxplotlqw 

7 Z-tiOXPLOT  X\L\U\ti 

[1]  X-(  HE~20+U-L)*(  (y*^r  / tX)-LENGTtl*-L*-U  ,X  ) + tX*LEUGTH~\ 

[2]  0. 1XL0. 5 + lOxL;  (orifi/VCfrtf-ejp  ' * ;0.1*L0.5  + 10xy 

[33  -*-9  * i 0=i/  •»£ 

[4]  'Z+-LENGT  ti  FILL  X 

[53  *«-•  "I  r I l‘  * C 1 + Z 3 

[63  B*-b,'  *|*|*  »x  ooooooooo  •••••••••'[ 1+Z 3 

[73  £}«-£,»  _ | L | | _ 23456789  2 3 4 5 6 7 8 9 ' [ 1 +Z 3 

L 8 3 *0  x p p Z*~  ( 3 « L ERGTH  ) P fl 

[93  'RANGE  OF  DISTRIBUTION  IS  ZERO' 


7 Ctil  SQUARE  1 03  7 
V T CHI SQUARE  XiDFlC 
[ 1 3 -+3*  i ( ( pX  ) -pT ) 

[23  -+0xpl J-'LENGTH  ERROR'.  X AND  T MUST  BE  THE  SAMP:  SIZE.' 

[33  0 + /(  (X-T)*2)iT<-T*  + /X 

[43  'CHI-SQUARE  VALUE  - ';C;'  DF  = ';(pX)-l 
7 


7coDi,/?fi's[a3v 

7 Z*-C  0 DE  RES  X ; Y ;.V  ; S ; Lt7  ;£/{/ 

ci 3 tf+-py+-yUY«-.x3 

[23  3-0. 5xy[L0. 25x3  +/V*i33+y  [To.  25x1  +//x\33 

[33  Lt7— 3 [ 1 3 -D+-Q,  [ 3 3 ■»  3 [ 1 3 
[43  £/3-3[  33+k 

[53  y-u>L3)  + u>3[i3  ) + u>3L3]  )+x>uq 

[ 63  2-*- ' -* » ° ++ 1 [ i +y 3 

[73  Z*-((pZ[;l]),2xpz[l;3)p(,Z)  MiMpZ)p'  » 

7 
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VCOMPAREl D]7 

7 Z+COMPARE  X W iL;Y  -.Mir  iJ  iP 

Cll  X+-(  HE~2  0+l'*a-L)*(  ([>£'«■  T /Y)-DFPTH*L<-1/Y  + (MISS*,X)  / ,X)*X*DEPTH-1 
[2]  l4-iilE~20*-L'nL)*U'*DEPTH*L)+MISS*DEPTH-l 
[31  Z*-(DEPTH,6*C-pXlli  Dp'  * 

[4]  J+0 

[ 5 ] J J + 1 

[61  P+DEPTH  FILL  Pi  hP+(  P*M ) lP+X  [ ;«/  1 1 

[7]  Z[  ;‘2^6*«/l^»  | ”i_T  | '[1-Pl 

[ 8 1 Z[  ;“l  + 6*c7l^'  *T  = 1 = |*  OOOOOOOOO  • ••••••••'[  1-rP] 

[9]  Z[;6k^1^'  | ”l_T | 23456789  2 34567 89 ' [ 1 *P 1 

[10]  >5 *\J<C 
[111  Z-«Z 

7 


VCCNDENSEIU 17 

V R-C ON  DENSE  X ;I  ;/Y  -tNR  -,U ; fi ; Z 
[11  P<"\0 

[2]  UTCOND  X 
[31  -*-0*iU£'/4*2) 

[4]  F«-(JV/?.2)pZ 

7 


7rCA/ri/VG£7i/rY[Dl7 
7 CONTINGENCY  F ;P  ;T  {S’  ;£  ;f  RI2  ;OF 
[11  FWF 

[2]  T++/C-+/F 

[3]  E-Co.xRtT 

[41  CMI2*-*/.  ( (Fn£')*2)t£ 

[61  DF-*/( pF)”l 

[6]  ' CHI-SQUARE  VALUE = ' ;ffll 2;'  £>F  = * ;0F 

7 
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1 cot DENSE C037 

7 R-*-C  ON  DENSE  X ;I  ;N  -,NR ; l' ; Q ; Z 

[11 

[ 2 1 UTCOND  X 

[3]  +0*\(NUM*2) 

[4]  R-(NR, 2)pZ 

7 


VCCAri/VGfWmDlV 
V CONTINGENCY  F iR iC ;T iE iCUI2 ; DF 
[1]  F-+rfF 
[21  T+-+/C++/F 

[31  E+C°.*RiT 
[ 4 1 CHI2++/, ( (F»F)*2) *F 

[51  £>F«-*/(pF)-l 

[61  ' CHImSQUARE  VALUE-  'iCHl 2;’  DF  = ' ;0F 

7 


7C0flflffLil2T0ff[[]]7 

7 R+CORRELATICN  N iZ  iC -,S ‘.CgiME  ANS ‘,V  AR 
[ 1 1 C - ( $Z  ) + . * Z+V  « ( p(7 ) p (ME  ANS+-  + 4W  )*ltp(~2  + l,l.pf/)  pW  + pCH+- ' * 

[21  R+'BFB.2'  FWY,(0,\p5),[ll(\pS),rtFo.x5-H(/i4F5-»-  + /Z*2)* 

0.5 
7 
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VF«r[0]7 

V 0L-£  FMT  RiSiUiLiGiXiTiKiJ ;M;G ;P ;D ;W ;0 ;L ; B ; / ;Cfl 

[1]  tf<-£«-14M«-pP*-(  ir_2+pP)pP 

[2]  G/X  ( 1 = l+Af ) 4-  1 0 *M+-M*2*U+-l<pCti-QJL%  )pA-’  0123456789.  ' 

[3]  +S«i  (N+0=N)  vy<-l2ipS«-t£ 

[4]  L0  :-^~i/v(  xp^4xg  = px^p^^-'  » ) av/(  » >4  ' .Of’D' ) eS 

[5]  -*-(L0  + (7«-0=pS-*-J'+S)+iS=M[2j-«-l),£“(  + P*~  ’ >5  * eK+K . (J+-S\  ’ , ’ ) 

+5 

[6]  +E+*pS+'TEXT  DELIMITER ' 

C 7 ] ->-L3-3xx(pG-K=K-<-(Ke~14A)  /K)  Lf/«-pX«-(  pK«-U iG ) ) + ( - ( ) iO ) * X 

[8]  L:+(.D-~l+G-KeE)  /L3-2x(pK)*W*-l+0*-'  XA'  eK+(.~Xe  ' ,')/K 

[9]  -*-L3x  i (fl*  + /G)*x«[  2] -<-101  | 1 - A i ( 1 1-G  \ 0 ) +/( 

[10]  -+L3-t>0%  -(.L<-'  EFI'  eK)  /xtf«-10l  | 1-A\(  | 1 -G*-B\  ' . ' ) + B-*-(  1- ( 4>G  ) x 0 ) iK 
Cli]  4-^(1 + P^(  (ir P4)l(w[i]-fl),y)t4)4>4 

[12]  L3:+(HD*\H*~'X'  eK)  ,E-pX*--W  ,D*-OpP*-(  ( tf-ff , 0 ) x i , y ) pJif 

[13]  -*-Z,4-i~H'Z/,<2<-1  + P/?<-(0  1 xpP^p[;iW[2]«-aL'if[2]r3x7AD])4fi 

[14]  P«-P*10*L«-[  10*|  P+0  = P 

[15]  -+-L3  x \Q-.J  «-+/  B«-(  ,B,cK),*0=P-‘-(L0. 5+.Vx,P)  tW«-1  0*Z)-«-l  Oi  | 1-AiG  + P 

[16]  L4:->-(pl4pL)/P-pp^-<-(l  0 xpG<-J  pT\'~  ' )pJ*-J  ,04-v/T*-0>P+B/P 

[17]  +( x/xor  w«-*z » eX)[ . x~r«-(  j-+g^i+[  io*ir  |p)>o+L«-i/-z>+o+~2+z,)/L/p,p, 

I 

[18]  -+E+*pS+'  FIELD  VI DTE' 

[19]  -*-L4  + l + x(  U[2]«-£v.<0)+G«-l  + 10r  . <|  Z,- ( B/ , Z, ) +7’«-l  0 = | P ) >V -0  + 0 + 

3 

[20]  T* — J +-P[T  / \1+J  J^L^plpA’^’P*  , * +0"  ' [-7p2-xp],  A[  l + «?(OplO)T  |L] 

[21]  F:+(J  v2>D*~'T'  eK)  / 1 , N+pX+M  11, 1 +«S?  ( Op  1 0 ) T T N*l  \ |P],X 

[2  2]  0«-,  (Wm(Z?r.x$X[;2+P]*l+A)o.<0-M0-l 

[23]  X+NpX,XLD/ \pX+,X]+'  ' 

[24]  I:-+(J+-jvo  = + /0+-or  Z/-0 ) /I  +pD+pP+-G  ,A[1  + $(Z/P10)tL  | p] 

[25]  P+D  p ( ,0*-0<|>0)\(  ,0«-0  « . <(  -G  )<t>  i£+G«-l  + pG  ) / , P 

[2  6]  -*HD~\J  vL+-~ ' L ’ eK , PlT  / \ D*-l  iX+pP+P  ,£;]•*-'*  ' 

[27]  P+Xp(,4>0)\(  ,0+~Xt~0)/,P 

[2  8]  + ( ~R  ) /£  -B+-1 , D-*~QpP-*-B  V ( 0*  X*-V  *1-2*L) iP 

[2  9]  HD:CH*-(  pK-+Cl+0+-O,  (M[  2 ] [ p0 ) p 0«-(  ' , ' =CR)  / \pCH)pCH)$Cd 

[30]  0+-,(M[  2],*)t  0 "l  +(W[2],e)p(,4>D®.>xfi^r/P^l+D-~l4>D)\K 

[31]  ->(LO-/Ax3),pOp^-GL,  ( (l  = l+#)+tfxl,tf)pZJ,  ,p 

[32]  E:K*-'NO  VALID  E,  I,  OR  F PHRASE' 

[33]  * FMT  PROBLEM  ?,K;(l,pS)pS 

V 
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VFIILtCnv 

7 Z+-W  FILL  X;JV;0;X1;X2;X3;X4;Y;I  -,J -,S 
Cl]  W-pX-XCU'-.X] 

[2]  2^0.5xXCL0.2  5x3+A*\3]+XCr  0 . 2 5*1  +.V  * i 3 ] 

C 3 1 S+-(Z  [l]»0.5  + 1.5*Z?),(Z[l]«*0«  S+D+Z  [3]wZ[l]),Z[l] 

[4]  5-L0.5+5,Z[2],Z[3],(Z[3]+0.5+0),Z[3]+0.5+ 

1.5x0 

[5]  Xl-(X2s5[l])/X2-(XS.?[2])/XH0.5+X 

C 63  X2-(X2>5[l])/X2 

C 7 1 X4«-(X3!>5[7])/X3-(X*5[6])/X 

[8]  X 3-*-  ( X 3 <5  [ 7 ] )/X3 

C 9 1 Y«-(  (X>S[2])xX<5[6])/X 

[10]  Z^f/pO 

c ii]  zcr>2o+9L-*-,i(~i,xi)o.=i<-\(5-orsu)Ci] 

C 1 2 ] Z[I]-*-10+9L+jK~l,X2)«  . = r-*-5[  l]  + i 1+5C  2]  w5C  1] 

Cl 3]  Z[I]«-10-*-9L+jK'l,X3)  o.=I^_l+5C6]  + il  + 5C7]-.5C6] 

[14]  ZCI]^20  + 9L-ri<X4o.=J^‘i^5[7]  + \l^VI.5C7] 

C 1 5 3 8 x \ 0 = pX 

C 1 6 ] zu«-*itn<-z[i<-i+y]-9 

C 17]  Z[I-mOU«I+1]«-8 

C 18 ] ZC5C3]+\0r5C5]^fC3]+l]^6 

[19]  Z[5[3],5[5]  ]«-0 

[20]  Z[5[4]]-l 

[21]  Z[5[3].S[5]>ZC5[3],5[5]]  + 2 4 
7 


7JA/Pi,J7[G]7 
7 INPUT iJ %M l Z 

Cl]  ' ENTER  NUMBER  OF  TREATMENTS 

C 2 3 c/-*-0*X-»-Q 

[3]  » ENTER  IE.CZQ.R  OF  NUMBER  OF  CBS.  FCR  TREATMENTS  1 TO 

[4]  Y-(  (M»X*-r /»♦□)  ,X)pO 

C 5 ] +l*\(K=pN) 

[61  -mxppU*-' RE-ENTER  COUNT  VECTOR , ONE  ELEMENT  PER  TREA 
C 7 ] +0x,  U<J+J  + l) 

[8]  ' ENTER  * ; 2W C eX 3 ; * OBSERVATIONS  FCR  TREATMENT  ' ;J 

C 9 1 -►7xpppY[;€/]-«-(Mi4X,l)pZ,  [M+MAX  *pZ*-Q)  pO 

7 
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r T7~ 


7XF[J17 

[11 

7 P AF  C/;Z  ;/V  ;A';FF 

; /V  Cr 

;A IAX  ’,SZG:\i ; - 

U*-nr*-\  >*«-0  *.v«-p £•*■ 

ZU 

Z+F.G 1 

[21 

-MlM  (tf  < '-ri) 

[31 

iuF*-tiF , ( * /FiX*-Z  [ 

**  1 ) * pF 

[41 

i#'  e +-u  F , ( ▼ / 1 < X ) i pF 

L S 1 

: J G +■  ( ♦ / Cr  S A’ ) * p Cr 

[61 

SZSiV+hnF*-  ,PjF°  . - 

(V  G , 

( / .7  < X ) i p (j 

[71 

>i'dX*lt(  (,VF  = 

.-,/iy 

-r /UF+\.'ib)  /SZGm) 

L 8 1 

-►2  ,/VF*- 1 0 

L 9 1 

/lit'  ;>j  /S  A F ( X ) - 1/  ( X 

) = 

'■ML 

[101 

' ,x\v  m)-5U) 

- 1 . 

• 

l/u 

V 


ilme  caiv 

7 Z+LZNE  W -,X  iY  iN  iMiP-,1  -,J  iKiXUYUX2iY2iF  iSLiYI 
[ll  -►3  * t 2=  p p<7 

[2]  V«-$(2,pV)p(  lpV),V 

[3]  F-U<-pY«-*[P+4Y*-V[  ;ll  l)p*-5Z^0 

[4]  JCl<*-0. 5 *X[IH  0.5  *«■*■!  0.5*l*,'4«-[0.S-*-/H3l 

[5]  Jr2«-0.5**[#*l*n**[tf  + l«en 

L 6 1 Yl-0.5*Z[Il-'-<Z<-ZUZ«-M  + Y-i7[P;2l.-Fl)[en 

[71  Y2<"0.5*Z[Il-<-(Z«'ZUZ-(-/'«)+Yl)[«n 

[81  F+-(SL*-SL+(Y2-Yl)*X2-Xl)*X 

[91  X-X* 1 
[101  -*>6*\X<2 

[in  yi-o.5xY[Lo.5*A/+n*(y*y[4y^[P;2i«Fi)[ro.bxw-n 

[121  'SLOPE I ' iSLi'  Y "INTERCEPT i ' lYI 
[131  ZH7[  ; 2 1 «F*-YJ-*-5L*V[  ; 1 ] 

[141  +16x\ kSELECZ-1 

[151  + 0,  ,Z-$(2,A0pF,Z 

[161  Z-$(2,JV)prf[  ;ll,Z 

7 


[11 


7LITCG17 
R+-LZZ  A 

F*-,'  0123456789*  [!♦(  ( 1-^L  10*i4 ) pi  0)  T/J 1 


99 


•<n'  . 


VLSLINE [Q]V 

V Z+-LSLINE  XiNiB 

[1]  +Al*\2=ppX 

[2]  X+t$(2,N)pX+(\N*pX)  ,X 

C 3 ] Ahb*-X[.  ;2]|fi(fl)(  2 )p(  ,-f[  ; 1]  ) « (M+lipX)pl) 

[4]  * SLOPE : »;fl[l];'  Y ■’INTERCEPT : * ; 6 C 2 3 

[5]  Z-*(2.#)pU[sl]>,Jr[s2]*(B[2]+Jr[  ;l]*fl[l]) 

V 


Vtf£0P0LJStf[[]]V 

V Z+MEDPOLISH  Xik\C\S\RI  CI\*.  A\S  \RG  \RE  \C  E \TV 

[I]  /?— pJf C 5 ID 

[ 2 ] /V«-ft'x(’-«-pY[  1 ; ] 

[3]  X+X^TV+O.  5*Z[L0.  5x/i/+l]  + (Z«-ZUZ-,Y]  )[T  0. 5*0  + 1] 

[4]  RE-R  pO 

[5]  CE-CpO 

[ 6 ] Q--S*-*/  | ,X 

[7]  Z-(A.r)p(,X)U,Y  + ( ift^.xfp  lxffG«-(r /,*)’[/,*] 

[8]  X*-X*(RI*’+/0. 5*Zi  ;L  0.  5*f  +1]+Z[  ; T 0 . 5 *f  + 1 ] ) o . xfp  i 

[9]  RE+-RE  + RI 

[10]  Z-<5i(C,fi)p(,$X)U,$X  + (/?pl)°.x(  iC)x*G] 

[II]  ftp  l)o.xCI«-+^0. 5*Z[L0.5*A+1;  ] + 2[r  0. 5xft+l  ; ] 

[12]  CE+CE+CI 

[13]  +ex\EPSILON<l-+/ I tX*S 

[14]  Z*-X 

[15]  -+1 7 x i REST  US*2 

[16]  Z~$l2,N)pU,RE*.*CE)*TV),,Z 

[17]  +21* \NORMEFPECTS=0 

[18]  RE*-RE~RA*-(+/RE)*R 

[19]  CE*-CE*CA-(+/CE)*C 

[20]  TV+TV+RA+CA 

[21]  ' TYPICAL  VALUE : ' \TV 

[22]  'ROW  EFFECTS : ' -,RE 

[23]  'COLUMN  EFFECTS-.  ’ ; CE 

V 


vmstatsvjiv 

7 MS  TATS  tfiZiCiSiCtLiliVAiiSiVAkSiHiJ 
L 1 1 F-«-($Z)  + . x2«-.y  » ( qW  ) pi'jE  AN  S+-  (■*■/¥)  f/V-r-l  + p ( 2 + l,l,pV)pi7f  pr^ 

[2]  'BF  8.2'  FM?(0,  ip5)  ,[ll( \p5)  , F *5  o . x£U- ( y/ltf  c«-+/2  * 2 ) * 

0.5 

[3]  'flMEANS  D.6F8 . 2 ' FAi?  AtfF/l/VF 

[4]  ' 25i  0.  Z?F/  ffl,fiF8.2'  FMT  S+S  *(/INl)*0.5 

[5]  Z--\J-Q 

[6]  -t-END*  \ ( pF ) + l 

[7]  C«-Jr[U«"J/[;/].l 

[8]  -*6  , Z-i-Z  , r-*-0 ,5xF[L0.25x3+,Vx\3]  +C  [ T 0.25*l-*-/V*\3] 
t 9 1 EltDtZ+-(  (p5),3)pZ 

[10]  'DL.  LSi3,6F8 . 2 ' FAJ™  Z[;l] 

[11]  'ffl MEDIANS  i3,aF8.2'  FMT  Z[;2] 

[12]  ' 0^  . QRTLS^,bF8 .2  ' FMT  Z[;3] 

V 


7A£'MF£'A*[[]]7 
7 NU to SUM  XiN iQiTl;Cg 
[ 1 1 FF 

[2]  • NUMERICAL  SUMMARY ' 

[3]  7 0p  1 ■»  * 

[4]  /*-•-, *1 

[5]  ' SAMPLE  SIZE  = ' ;,V 

[ 6 ] 7 Op  1 n 1 

[7]  ^-*-0.5xX[[0. 125*7 +A/x\7]+^[T  0. 125xl  +/yx*7] 

[8]  JP1^Q[4]  , 0,Q[4]  ,0,0,  ( (Q[2]+6[6]  ) *2)  , Q [ 2 ] , 0 

[9]  Tl+Tl  , fi[  6]  , (Q[6]r>Q[2])  , ((Q[l]+Q[7])*2) 

[10]  n^ri,fl[l].0.fl[71,  (fl[7l*fl[l])  , ( U[1  ]+*[#])  *2) 

[11]  ri^-ri.mi.o.xm.xm-rm] 

[12]  T 1-  4 5 pTl 

[13]  £H-' MIDPTS , , LCQ/8/MIN , .MEDIAN , , UPQ/8/MAX , .SPREADS,' 

[14]  'BF10.2.C  | l5'  FM!T  1”1 

[15]  7 Op  ' it  ' 

7 


7FA'FF[G]7 
7 H+-ONEH  XiA-.B-.C 
[ 11  ^0,X,  0 

[21  B«-X,0,0 

[31  r«-o,o,x 

[4]  H + ~24(C+B  + 2*A)  m 

[51  ««■*[!]. 2*fl. Up*! 


Cl] 
C 2 ] 

[3] 

[4] 

[5] 


Cl] 

[2] 

[3] 

[4] 
C 5 ] 
C 63 


yOA/£3CG]7 

Z+0NR3  A;Y1;Y2;A1;*2;*;C 
Y1<-(3XA[2]  ) ’2  XA[ 3 ] 

Y2*-{  3 XA  [iV-»  1 ] ) •» 2 XA[  ~2+.V-*-p  A ] 

Al+-Yl,”l  + A 
A2*-l  + A,  Y2 

Z<-(Y1,A,Y2)[(  i/V)  + l + (C=Al<A)-(ASA2)=f+-Al<A2] 


VPAHTIALLiHV 

V P-PARTIAL  f!  iC  iCINlMiY  iXiBETA 
14*-  x 1 pN 

' ENTER  COLUMNS  Ob'  MATRIX  TO  BE  PARTI ALLED  OUT.' 

/ll:-M2xiO*  + /~(  C+-D ) e M 
Y+-t7[  iCIN+-(~MeC) /Ml 

-*•0  , pP+CORRELATION  R*-Y -X  + . *BET  A+Y&X+-1 , (2^(p^C  ;T3  ),  l)p*^C  ;/"3 
/32:-Ml,p  'UNACCEPTABLE  VALUES , ENTER  INTEGERS  FRO’4  1 TO  ' 

SI 


SJSHOWRES  C G3  V 

7 Z-SdOVRES  X-,Y  iN  iQiDiLWiLFiUWii/F 

[1]  W-‘-pY-*-Y[iY-*-,A  ] 

[2]  ty-'-O  . 5XYCL0. 2 5*  3 +.V x i 3 ] +Y [ C 0.25*1  + Vxi3] 

[3]  LN+Ql  l]-H+-WC3]’(j[l] 

[4]  y^[3]  + Z> 

[ 5 ] LF*-L'N  ■» Z>  * 2 

[6]  UF+UN+D* 2 

[7]  Y- U>  LF  ) + (*>  Z-y  ) + (*>£C  1]  ) + (*>£[  2]  ) + ( A’ >^[  3]  ) + ( X >UN  ) +A  >UF 

[8]  Z+- ' 0o°  » o o xX  ’ [ 1+Y  ] 

[9]  Z*-(  (pZ[;l]),2xp2[l;])p(,Z)  MD,(pZ)p'  ' 

V 


102 


* 


VREGRESSiQlV 

V Z*-Y  REGRESS  X. \H  ;K  ;C  ;X  Pv  T NV  ; XPY  ; BET  A ; RSS \T  SS -,S2  ; ESS  ;V  I D ; DEP 

Cl]  Jr«-(2t(pA),l)p* 

L 2 ] X^(0,1^EINTERCEPT)*1,X 

L 3 ] XPXINV+®1$X)  + . *X 

[4]  BETA+XPXniV+.*XPY*-(<S)X)  + . *Y 

[5]  RSS-(  (§BETA)  + .*XPY)*C+-(  ( + /Y  ) * 2 ) */V«-p  , Y 

[6]  ESS~(TSS*-(  (SiY)  + .xY)-C)^RSS 

[7]  S2*-,BSS*(.!t*l)*K4-(p,BETA)-'EINTERCEPT 

[8]  CR 

[9]  ' AtiOVA' 

[10]  ££-' SOURCE, DF, SUM  SQUARES, MEAN  SQUARE , F^RATIO' 

[11]  ' ' 

[12]  ,i!m,Gff£SSI0fflI),I4,B£16.4*  FMT  ( A ) , ( , RSS ) , ( , RSS  *£)•(•  RSS  *K)  iS2  , 

L13]  C.H *■ ' ' 

[14]  '7]  RESIDUALS, IH,BE1B.E  ' FkT  ( ( U- 1 ) -K) , ( , ESS ) , S 2 , 0 

L15]  ' I\T0T  AL  0,14, ££16. 4 ' FMT(N-l) , ( ,TSS) ,0,0 

[16]  » ' 

[17]  'R  SQUARE : ' ; ,RSS*T3S 

[18]  'STF  ERROR:  '-,,S 2*0.5 

[19]  COEFFICIENTS ,T  STATISTICS' 

[20]  ' FI  5 . 4 ' FMT$(2,p,6£T/1)p(,££-7'/3),(,££7M)*(l  1 4?  2 x A'i3  A JT  V i/  )* 

0.  5 

[21]  ' DO  YOU  V AUT  A PRINTOUT  OF  THE  V ART  ARC E-COV ARI ARC E MATRIX?' 

[22]  +A1*\ 'Y  ' *1+0 

[23]  ' VARIARCE^COV  ARIANCE  MATRIX:  ' 

[24]  '£12.4'  FMT  V 

[25]  A 1 : * OUR  SIR -V  AT  SON:  ' ; (+/(  (1  + ,C  )’  ( ~ 14  ,C ))  *2)  * + / ( ,C*-Y  "X+ ,*bSTA)  * 

2 

[26]  Z+$(  2,,V  )p  ( ,X  + . *5ETA  ) , ,C 

[27]  Bl:'UO  YOU  ,/  AUT  TO  FORECAST  A VALUE  FOR  Y ? ' 

[28]  >C1m  'Y  '*1+0 

[29]  'ENTER  X VECTOR  ( ' ; A ; * VALUES )' 

L 3 0 ] 'FORECAST  OF  Y VALUE:  ' ; (C-(  1’  El  INTERCEPT ) 4 1 ,J)  + . *RET  A 

[31]  'VARIANCE  OF  FORECAST  ERROR:  ' ;S2*1+C+. *XPXIRV+ . *QC 

[32]  +E1 

[33]  Cl’.'DO  YOU  J ART  TO  SCAT  RESIDUALS  VS.  PREDICTED  Y?  ' 

[34]  +0xi 'iV  ' =1+0 

[35]  DEP+-Q  . 5 xQlD+-l  / 7 0 , ( \ / ( ( 0 . 7 5 *.V  ) , 3 0 ) ) 

[36]  SCAT  Z 

V 


103 


VSfyUTLLnV 


V N+SCAT  ZiN  iX;Y  iCihiU  iSiLil  iJ  iKiUTirLiG:DitiiA;CiV 
[11  ■+3x\(2  = + /2  = /!/)v(x/^)>i-//|/-t-pZ 

[2]  Z^$(2,pZ)p( \pZ) ,Z^,Z 

[31  J,-Z[;lf  tr«-~l  + (pZ)[2ll 
[41  fi«-pZ-Jt-,Z[;ll 
[5]  W'-5«-  2p0 

[6  1 J*-l+Q*p{D+-NDlVX  >NDIVY)  tb+-WID,DEP 
[71  t'?-l  0*  [lO^rL-lE"  2 0-^  (.U'lJ  W /Z)-SW\-l/Z)  *DlJl 

[81  5[s/]^'rxL5[Jli['7,^['r[l  + t|f,£,nt'^(l  2 5)*t'n 

[9  1 UiJl+l’T*U,tJl*UT 

[10]  L[«n«-1+Gx[  (BlJ  1 -1 ) *G-( £'[•/]  »S[-7]  )*UT 

[111  Z^.Y 

[121  ■+7*\3>Ji-J  + l 

[131  /K(4>L)pO 

[14]  y^l  + [0.5  + (Z,[llnl)x(^,5[n)t£/[l]pf[n 

[15]  y-li-[0.5  + (i[2]-l)x(rn5[2])t£,[2l^5[2] 

[161  1-1 

[17]  -»-20x\l<r 

[18.1  i4[y[I;  11  ;y [II  KlOLi4[y [I;  1]  ;y[m  + l 

[191  — 1 8 + 6 x R<I-*-I  + 1 

[201  7-1 

[211  O-0=y-/3[Y[I;7]  ;Y[I]  1 

[221  i4[Y[I;7l  ;Y[Ill-(10xy>K+l)+(  (K-l)xtf=y)+  (K-3  5»*  2 *J  ) *D 

[23]  -+2 1 x \ if  sZ— J + 1 

[ 2 4 1 -21x  \CzJ-*-J  -Z-l 

[25]  C—  ( 4>P/0  L 1 T l■,■[0.5■*•(L"l)x5i5"^., 

[261  i4[;G[m^i4[;G[ll]  + 36xO  = i4[;G[m 

[27]  /l[G[2l;l-/l[C[2];]-35xO=/3[C[2l;l 

[281  V-*  <>2345678  SffiLkK&J  J,I  IHHGQ.F  £E  [HD^C  CB  £A£w  | ’ [ l + S/4  2 

1 

[29]  'RANGE  CF  Xt  »;S[l].£'[ll 

[30]  • RANGE  OF  Yt  »;5[2l,£'[2l 
’ 


VSPLI r[Q]V 

V Z-SPLIT  KiPiQiRiSiT  il  ;J  iC 

[1]  ?+~  3 + a: 

[21  Q-~2> 1*X 
[3]  R+-~  1 + 2VX 
[41  P+-3+JT 

[5]  I«-(  ((P<«)xfi->5)-r(P>e)x/?<5)xg  = i( 

[61  C+U+Xill,X,XU-\)ll+I+U*0)/I+I*\~3+M  + pX-\ 

[71  P-f(3xe)w2x^[J.l 
[81  P-X[I+2l 

[91  r*(3x^Jr[I  + 4l  )«2xiSf[J+5l 

[101  *[I+2]*-(fl*W)*(Px“W)+G*0=«K(r=P<G)'»(GSi*)=C-'-P<« 

[111  xIIt3’\*-(T*1=J)*(P*~1=J)+S*0=J+-(.C=P<S)t,(S<T)=C+-P<T 

[ 121  z«-i  + “i+x 

V 


[11 

[21 

[31 

[41 

[51 

[61 

[71 

[81 

[9] 

[10] 
[111 
[121 

[13] 

[14] 
[151 
[161 

[17] 

[18] 
[19] 


75Y,/4riPinir5[aiv 

STATISTICS  RiXBARiVARifii  iQiSiZ 

XBAR+(.+/R)iN+o  ,R 

VAR*-{+/(  (R’-'XBAR)*  2)  )*N~  1 

5-PUP] 

^■►0 . 5 x5  [ [ 0 . 2 5 x 3 -*-W  x \ 3 1 [ [ 0 . 2 5x  l+ffx  1 3 1 

'MEAN i 'iXBAR 
'VARIANCE  i '-,VAR 
•STD.  DRV.  t 1 ;5-»-i''i4/?*0.5 
•CCEFF.  CF  VARIATION  t » 

•LOWER  QU  ART  I LE  \ » ; <5  [ 1 ] 

•UPPER  QUARTILEi  ' ;Q[3l 
• MEDIAN i *;Q[2] 

•TRIMEAN  t ' ;0.25xQ[l]-rQ[3]-r2xg[2] 

•MIDMEAN  t • ; ( tPZ)x  + /Z^(  ( P>2 [ 1 1 ) a ( p<$ [ 3 1 ) ) /R 

•RANGE  t * ; ( T / R)  *>1  / R 

•MIDRANGE:  • ; ( ( T /R ) + L / R)  i 2 

'tffrffl  ABSOLUTE  DE.IATICN:  • ; ( + / | R»Ql  2]  ) *iV 

•INTERQUARTILE  RANGE t * ;Q[3]«<3[1] 

•CCEFF.  OF  SKEWNESS \ ’;((  + / (P*Xfl<4P)  *3)  i/V-1 ) *5*3 
•CCEFF.  OF  KURTCSISt  ’;(<(-■/ U**S4J?)  *4  ) f.tf-1 ) *5*4 ) «3 
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VSTEMLEAFlUW 

7 Z+STEMLEAF  X ;C -,R  ;S  iSI  ;I  ;J  ;F  i A ‘,W  i L iWW  lAA  iXW 

[1]  r-io*i*Lio»^i£"2o+u[p;ni-u-^u^,;n ) cn ) i scale 

[2]  5I«-(l-*-3x+/(flxr)>  25  50)W(p*)>  25  100 

[3]  X-LO. 5+^xrxiO*+/5I=  236 

[4]  F-  + /(.SI=\ 9)/  0.5  2 1 1 0.5  2 2 1 0.5 

[51  i4-'-20p  1 01234567 89ABCDEFGHIJ  1 ,Z*-'  ' 

[61  I«-fxU[ll*10xF 
[71  Xd+Z+WU+WIDTH**  . 

[8]  S+'m  ' [ W-Hll^O] 

[9  1 L+- 1 + ( | W+-  Uf^(10xJ)+t/x~i  + l0xF)/^)nl0xL  |J 

[101  A A*- A [ 1+  ( 1 0 10 ) T | III 
[111  + ( (pALLl)ZM)/Ll 

[121  Z-Z,  *(7p5,  AA,  ' | (VV  + >J[Ll ( LI?  ( p-4  [ L 1 ) ) ,2p  ’ ' 

[13]  -*L  2 

[14]  LltZ-Z.JWpF.yM,  ' I'  ./ULl.Wp'  ' 

[151  L2tI^I+Fxi-(J=0)  x^[  H<0 

[16]  ■>8xl0<pA'^(pV)4A' 

[171  Z-(  (l  (pZ)**(7)  ,Xl7)pZ 

7 


75£,/'l5Q[Dl7 

7 SUhSQiC  \NUMBER\B\BLDF 

[il  r«-+rfir 

[21  ?55«-(i-/(+/y*2)  (+/T)*2)*Nl'MBEk4-+/N 

[3]  BLSS+(*/(((B+-+/Y)*2)iK))-C 

[41  MSBL-BLSSiBLDF+MAX-1 
[51  ■+7*\(OPTICN  = 2) 

[6]  BLSS+bLDF- 0 

[71  WTSS-TSS-BLSS*BTSS<-(+/ ( (T*2)*N)  )-f 
[ 8 1 MSR+BTSSiK* 1 
[91  TDF+NUMBER-1 

[101  F+MSR  *MSE*-rfTSSi  E DF-*-T  DF  ~B  LDF  + K-1 
[111  BLCCK-(B*K) *YBAk+ (+/?)* NUMBER 

[12]  ' TABLE1 

[13]  ' SOURCE  DF  SS  MS  F ’ 

[141  ££«-" 

[151  'D  TREATMENTS  ,15  ,F13.2,f,11.2,F8.2'  FMT(K*1 ) , BTSS ,MSR  % F 
[161  -*•(  OPTION^  1 ) /L5 

[171  'B  #LCT/£50 ,18  ,fl3.2,F11.2,F8.2  ' FMT (MAX-1) ,BLSS ,MS6L , (MSBLiMSE 

) 

[181  L5i'D  f^i?(?i?[!l,I9,F13,2,fll.2  » FAl?  EDF  ,WTSS  ,MSE 

[19]  »fl  ^<?^/l£^5,I9,F13.2,  FM?  TDF ,TSS 

[201  »l!|F«5Q£'i4FF  = E.F5.3'  FMT  ( BTSS+B  LS  S)  *T  SS 

[21]  »B OVERALL  MEAN  - 0^10.2'  FMT  YBAR 

[22]  '^TREATMENT  EFFECTS  ft.FB.2'  FMT (T *N) -Y BAR 
[231  -+(CPTICN  = l)/0 

[24]  'flBLCCK  EFFECTS  B.F6.2'  FMT  BLOCK 
7 
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vry/caicnv 

V Z-TVICfc  X\Y 

Ll]  Z*-Y+A3RSh  X-Y-A3KSR  X 

V 


VJTCOUD  [J]V 

V UTCOND  X 

[ 1 ] X*-(  2p  ( pA  ) , 1 ) pX 

L 2 ] Z-'' 

[3]  r-*-0*A/F-*-(  p*  ) [ 2] 

L 4 j rt-p£W  [ kU*-(U*MISS ) /'J+X  C ; J+  1 ] ] 

[5]  y*0.5xy[L0.125x7+.Vxi7]+(/[r  0 . 1 2 5 x 1 +.V  x i 7 ] 

C 6]  -*-7  + 2xL/vy;n2 

[7]  Z«-Z,a[4] 

[8]  +4  + 11*  HR-J+-I  + 1 

L 9 ] Z-Z,Ql  4],iC6]-si[2] 

[10]  ■*•4  *N  R>  I *~I  + 1 

[11]  Z-Z.tf[l3.3[2]#S[4].3C6],i/[i»],V 

[12]  ■♦4  + 1 3*N  R-I+-T  + 1 

[13]  Z-Z,tf[l],3[l  2 4 6 7],£/[/V],W 

[14]  +4  + 1 6*// A = 7 <-J  + 1 

[15]  ££■*- ' MEDIANS  ' 

[16]  -*•  0 x p []*•  ' F 1 0 . 2 ' F&T  Z-UF,  DpZ 

[17]  Z-UF.oJpZ 

[ 1 8 J QJL*-*  MIN,  LG<u, MEDIAN  9UPu,  VAX,  SIZE  ' 

[19]  +0xpi>’F10.2  * F'4?  Z 

L 2 0 ] Z ■*■(/»'  rt  , 8 ) p Z 

[21]  ' MIN  ,LO 8 , LOW  , Mb,  D I AN , UPQ  JJP 8 , , SIZE  ' 

[22]  1>'F9.2'  b‘MT  Z 

V 
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